1

我正在 R 中进行彻底的“R2Winbugs”包中的贝叶斯分析。在 R 输出中,估计的参数被四舍五入。我的问题是如何控制这个?

这是我的命令:

bugs(data, inits=inits, model.file = "C:/Users/Gunal/Desktop/dummy/dummyw.txt",
    parameters, digits=5,
    n.chains = 3, n.iter = 1000, codaPkg = FALSE,
    bugs.directory = "D:/PROGRAMLAR/WinBUGS14/")

这是我输出的一部分

         mean  sd   2.5%    25%    50%    75%  97.5% Rhat n.eff
beta1    0.3 0.1    0.0    0.2    0.3    0.4    0.5  1.7     5

我想要的是将 beta1 视为 0.30000。我试过digits=5,但没有用。有谁知道该怎么做?

任何帮助将不胜感激

干杯

居纳尔

4

1 回答 1

1

是的,这很简单 - 使用print(out, dig = 6)(或 7, 8, ...) 打印结果:

out <- bugs(....)

print(out, dig = 6)
于 2012-12-07T15:18:50.123 回答