CosmoMC: Run stopped before specified convergence achieved

Use of Cobaya. camb, CLASS, cosmomc, compilers, etc.
Post Reply
Agne Semenaite
Posts: 19
Joined: June 08 2020
Affiliation: Max Planck Institute for Extraterrestrial Physics

CosmoMC: Run stopped before specified convergence achieved

Post by Agne Semenaite » 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!

Antony Lewis
Posts: 1936
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: CosmoMC: Run stopped before specified convergence achieved

Post by Antony Lewis » May 21 2021

The run didn't time out on the server?

Agne Semenaite
Posts: 19
Joined: June 08 2020
Affiliation: Max Planck Institute for Extraterrestrial Physics

Re: CosmoMC: Run stopped before specified convergence achieved

Post by Agne Semenaite » May 21 2021

No, that was what I also suspected initially...

Antony Lewis
Posts: 1936
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: CosmoMC: Run stopped before specified convergence achieved

Post by Antony Lewis » May 22 2021

Have a look at the chains' .converge_stat file and see if it is consistent, and check the .inputparams actually has what you expect.

Agne Semenaite
Posts: 19
Joined: June 08 2020
Affiliation: Max Planck Institute for Extraterrestrial Physics

Re: CosmoMC: Run stopped before specified convergence achieved

Post by Agne Semenaite » May 25 2021

The .converge_stat file shows a value that is ~an order of magnitude higher than what is specified in MPI_Converge_Stop in .inputparams...

I might just stop the chains manually rather than relying on MPI_Converge_Stop if there is no obvious solution to this behaviour. But I really appreciate you took the time to reply Antony!

Post Reply