CosmoMC: Run stopped before specified convergence achieved
Posted: May 20 2021
Hello,
I am running some chains on CosmoMC and I would like them to have the convergence of R-1~0.005, which I specify at MPI_Converge_Stop. However, for the cases when I am combining Planck with multiple datasets, the final chains seem to have finished running before this convergence is achieved (at ~0.02)
I am aware that getdist measures convergence differently than it is calculated during the run, but even if only use the second half of the chains with no derived parameters, the Gelman-Rubin parameter of final chains is larger than what I specified (~0.04)
Has anyone faced this issue before and have any idea how I could work around this?
Thanks a lot for any insight!
I am running some chains on CosmoMC and I would like them to have the convergence of R-1~0.005, which I specify at MPI_Converge_Stop. However, for the cases when I am combining Planck with multiple datasets, the final chains seem to have finished running before this convergence is achieved (at ~0.02)
I am aware that getdist measures convergence differently than it is calculated during the run, but even if only use the second half of the chains with no derived parameters, the Gelman-Rubin parameter of final chains is larger than what I specified (~0.04)
Has anyone faced this issue before and have any idea how I could work around this?
Thanks a lot for any insight!