cosmomc : getdist marginalised probability
-
- Posts: 13
- Joined: May 11 2006
- Affiliation: Ohio State University
cosmomc : getdist marginalised probability
Hi,
I was under the impression that:
a) filename.margestats showed statistics from the marginalised posterior pdf.
b) the solid line plotted from filename_px.dat plotted the marginalised posterior pdf.
However, I am finding that there is a conflict between these. The value of param for which column 2 (which is the normalised likelihood) of *.dat =1 is not the mean value in file.margestats
I was under the impression that:
a) filename.margestats showed statistics from the marginalised posterior pdf.
b) the solid line plotted from filename_px.dat plotted the marginalised posterior pdf.
However, I am finding that there is a conflict between these. The value of param for which column 2 (which is the normalised likelihood) of *.dat =1 is not the mean value in file.margestats
-
- Posts: 1943
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: cosmomc : getdist marginalised probability
The maximum of the posterior is not the same as the mean unless the posterior happens to be exactly symmetric.
-
- Posts: 13
- Joined: May 11 2006
- Affiliation: Ohio State University
cosmomc : getdist marginalised probability
Thanks for that clarification.
1) I understand that the limits produced by margestats are based on having an equal number of samples on either side.
2) Is it possible to get limits based on the area enclosed at a particular likelihood from getdist or would we have to compute it from the chain(s) ourselves ?
1) I understand that the limits produced by margestats are based on having an equal number of samples on either side.
2) Is it possible to get limits based on the area enclosed at a particular likelihood from getdist or would we have to compute it from the chain(s) ourselves ?
-
- Posts: 1943
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: cosmomc : getdist marginalised probability
There's no option to calculate the other choice of limit definition at the moment.