MCMC法与观测数据得到的参数最佳值误差如何计算
How to calculate errors of best value of parameters that obtained from MCMC method and observational data
我有一个模型和一些观测数据。我使用 MCMC 方法获得最佳自由参数,并使用一些编码绘制 1 到 3 西格玛置信水平的等高线(如您在图中看到的那样)。我想要任何置信水平下每个最佳值的 sigma 的 +/- 值,但我知道它不是对称的。所以,我没有发现平方方差的通用方程有用。还有其他方法可以计算+/-误差吗?
this is my contours
this is what I want to get
我在最终的高斯概率中使用了 np. percentile(w, [15, 85])
、np. percentile(w, [5, 95])
和 np. percentile(w, [0.5, 99.5])
,并正确地找到了每个置信度的错误
我有一个模型和一些观测数据。我使用 MCMC 方法获得最佳自由参数,并使用一些编码绘制 1 到 3 西格玛置信水平的等高线(如您在图中看到的那样)。我想要任何置信水平下每个最佳值的 sigma 的 +/- 值,但我知道它不是对称的。所以,我没有发现平方方差的通用方程有用。还有其他方法可以计算+/-误差吗? this is my contours
this is what I want to get
我在最终的高斯概率中使用了 np. percentile(w, [15, 85])
、np. percentile(w, [5, 95])
和 np. percentile(w, [0.5, 99.5])
,并正确地找到了每个置信度的错误