I dind’t have too much luck with batch_size or obj_n_mc.
However, I noticed that while the pm.Deterministic() calculation of “cov” is wrong, if I simply calculate the cov matrix with the point estimates of the parameters after fitting, the corr matrix shows much less deviation from the true:
Left is corr matrix calculated post fitting, right is estimate calculated during fitting through Deterministic
