我有这段代码
out <- bugs(data, inits=inits,...)
并且其输出的一部分是
mean sd 2.5% 25% 50%
beta1 0.0 0.0 0.0 0.0 0.0
beta2 0.1 0.0 0.0 0.1 0.1
我想知道如何使用测试版。例如,我怎么能计算beta1 * 10。我试过了10*out$beta1
,但它没有用。
答案 0 :(得分:2)
函数bugs()
生成列表列表。参数的估计平均值位于名为mean
的子列表中。您可以按函数str(out)
检查输出的整个结构。
out$mean$beta1
out$mean$beta2