在coda包

时间:2016-03-23 16:22:33

标签: summary mcmc coda

我编写了一个函数,在一些处理后提取已经稀疏和后期预热的MCMC样本,输出是一个mcmc对象。但是,当我运行summary()函数时,它自然不会识别输出已经处理,因此在标题中显示不正确的信息,即迭代等。

我的问题是,当我调用summary(mcmc对象)时,有没有办法不显示标题信息?



> summary(output)

Iterations = 1:10000
Thinning interval = 1 
Number of chains = 1 
Sample size per chain = 10000

1. Empirical mean and standard deviation for each variable,
   plus standard error of the mean:

           Mean        SD  Naive SE Time-series SE
beta  1.052e+00 3.189e-02 8.914e-04      9.185e-04
df    3.849e+00 2.916e-01 8.150e-03      1.516e-02
sigma 1.056e-02 2.504e-04 6.998e-06      1.000e-05
tau   8.990e+03 4.273e+02 1.194e+01      1.710e+01

2. Quantiles for each variable:

           2.5%       25%       50%       75%     97.5%
beta  9.891e-01 1.032e+00 1.052e+00 1.073e+00 1.113e+00
df    3.304e+00 3.650e+00 3.836e+00 4.042e+00 4.450e+00
sigma 1.004e-02 1.039e-02 1.055e-02 1.072e-02 1.105e-02
tau   8.197e+03 8.700e+03 8.977e+03 9.263e+03 9.917e+03




0 个答案:

没有答案