Dealing with missing data and custom distribution

That doesn’t sound like an issue on the random side anymore but on the observed data. You still see the problem even without the missing value no?

What do you see if you write model.check_test_point()?