1

我正在使用R 中bayesm的包运行 Hierarchical Lin Regr 模型。我有一个包含一个依赖变量和 6 个预测变量的数据集。有 207 个独特的受访者，每个受访者有 35 个观察值。

我开始使用

print(out$betadraw)

然后我阅读了有关将sink其输出out$betadraw到文件的功能。我认为该sink功能会捕获所有的平局。相反，在一定数量的平局之后平局被截断。

我需要捕捉所有的平局。此外，是否可以将对象从bayesm一个包传递到另一个包以coda进行收敛诊断？任何帮助将不胜感激。

4

1 回答 1

0

如果没有可重现的示例，很难确定发生了什么。

您应该能够使用?file将open参数设置为write. 然后，您可以使用?write捕获输出并将其写入文件，并将append参数设置为TRUE. 以下在我的机器上运行良好： 

> zz <- file(description="some name.txt", open="w")
> isOpen(zz)
[1] TRUE
> for(i in 1:100000){
+   x <- rbeta(1000, shape1=10, shape2=10)
+   write(x, file=zz, append=TRUE)
+ }
> close(zz)

（注意，我不会尝试运行它；它花了将近半个小时，创建了一个 962 MB 的文件，只能用 EditPad 打开。）

于 2013-10-26T04:03:44.493 回答