Sorry I don’t think the random has anything to do with it. The sampling of missing values will still be done with MCMC.
Is it possible that proposed values ever evaluate to -inf / nan in your custom logp?
It could be a failure in the DensityDist as well. I don’t know if it was ever tested with missing values.