GetDist Convergence or Gelman Rubin test
-
- Posts: 3
- Joined: March 19 2020
- Affiliation: UNAM
- Contact:
GetDist Convergence or Gelman Rubin test
Hello,
I want to do a convergence test given two text files resulting from two mcmc samples (in cosmoMC format). How can I use the getConvergeTests function or getGelmanRubin function?
I suppose my main difficulty is that I can't create a mcsamples object containing the samples from both text files. Can you help me?
Thank you.
I want to do a convergence test given two text files resulting from two mcmc samples (in cosmoMC format). How can I use the getConvergeTests function or getGelmanRubin function?
I suppose my main difficulty is that I can't create a mcsamples object containing the samples from both text files. Can you help me?
Thank you.
-
- Posts: 1943
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: GetDist Convergence or Gelman Rubin test
This function will load a set of cosmomc chain files:
https://getdist.readthedocs.io/en/latest/mcsamples.html#getdist.mcsamples.loadMCSamples
https://getdist.readthedocs.io/en/latest/mcsamples.html#getdist.mcsamples.loadMCSamples
-
- Posts: 3
- Joined: March 19 2020
- Affiliation: UNAM
- Contact:
Re: GetDist Convergence or Gelman Rubin test
Hello,
Thanks a lot for your answer. I tried to use this function, but when I want to involve multiple files of chains, do I need to use a tuple, a list, a set or the comma separated file names? I tried all these four options without success.
If I use mcsamples.loadMCSamples with only one file, I don't have any problem, but with more files the error is
'set' can be 'list', 'tuple' or 'nonetype'.
Thank you.
Thanks a lot for your answer. I tried to use this function, but when I want to involve multiple files of chains, do I need to use a tuple, a list, a set or the comma separated file names? I tried all these four options without success.
If I use mcsamples.loadMCSamples with only one file, I don't have any problem, but with more files the error is
Code: Select all
AttributeError: 'set' object has no attribute 'endswith'
Thank you.
-
- Posts: 1943
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: GetDist Convergence or Gelman Rubin test
just give the base name, 'chain' will load files 'chain_1.txt','chain_2.txt' etc.
-
- Posts: 3
- Joined: March 19 2020
- Affiliation: UNAM
- Contact:
Re: GetDist Convergence or Gelman Rubin test
Thank you!