Transitioning from pm.sample() to pm.fit()

Same dataset.

I assumed it’s a feature of NaN handling in our VI code vs MCMC. MCMC will just give the sample -inf probability if it’s nan, whereas VI fails outright (as reported).