Page 1 of 1

Getdist: values in the .margestats file

Posted: March 09 2005
by Ariel Sanchez
Hi,

I am using the COSMOMC code to obtain constraints on cosmological parameters and GETDIST to analise the output files. I have been comparing my results with the outputs from other code and the results for simple six parameter models agree quite well in general, but with some subtleties. I have a doubt about what do the values quoted in the .margestats file really mean. I have interpreted those values as the mean value, the standard deviation, and the limits of the allowed regions for the parameters correspondent to the values given by the parameters "contour1" and "contour2" in the distparams.ini file (by default 68% and 95%). Is that correct?
I will also like to know the precise definition of the limits, that is: how are they calculated?
I hope you can help me. Thanks in advance

Ariel Sanchez

Re: Getdist: values in the .margestats file

Posted: March 09 2005
by Antony Lewis
Yes, that is correct.

Unless you have set one of the limitsxxx parameters to impose a prior on one of the parameters (e.g. that [tex]m_\nu > 0[/tex]), the limits are two-tail limits. Otherwise they are one-tail.

For the 95% limits this means that 2.5% of the samples are larger than the upper limit, 2.5% are lower than the lower limit. For a one tailed test, it would mean 5% of the samples larger than the upper limit.

So the limits are integrated probabilities. In particular, the marginalized probability value at the two limiting points will not in general be the same (unless the posterior is symmetrical). An alternative protocol would be to choose limits so that they have the same marginalized probability, with a given total area outside them (so the integrals at the upper and lower ends would not be the same). GetDist does not to this, largely because it is much harder as it requires an estimate of the marginalized probability density at a point (rather than just counting samples).

Getdist: values in the .margestats file

Posted: September 15 2015
by Cassio Pigozzo
Dear Antony,

One question about this topic: if I have a parameter with a limit and which PDF doesn't peak at that limit (and isn't a tail in fact), how would be the best way to proceed? Both one or two tails procedures you described don't seem to fit well to this case. Wouldn't it be possible to evaluate the cumulative probability density from the peak?

Thank you.

Cassio

Re: Getdist: values in the .margestats file

Posted: September 15 2015
by Antony Lewis
The latest GetDist follows the recipe outlined in the Planck 2013 parameter paper for whether to quote 1, 2 or no-tail limits for bounded parameters, depending on the marginalized density at the boundary.