1

RでJags/Winbugsを実行していますが、これは正常に機能しますが、出力をテキストまたはcsvファイルに保存できません(どちらも機能します)

TEST.sim<- jags(data=jags.data, 
                parameters.to.save=jags.params, 
                n.iter=200000,jags.seed=123, 
                model.file="~/Documents/ARRModel.txt")

必要なのは、次のように入力したときに表示される出力を保存する方法です。

print(TEST.sim)

write.tableを試してみると、次のエラーが発生します。クラス'"bugs"'をdata.frameに強制できません。

4

1 に答える 1

2

ライブラリR2jagsと関数を使用している場合jags()、関数print()はリスト要素BUGSoutputとサブリストに格納されている統計のテーブルを作成しますsummary。これらのデータに直接アクセスして他のオブジェクトとして保存し(または関数で直接使用してwrite.table())、テキストファイルに書き込むことができます。

jag.sum<-TEST.sim$BUGSoutput$summary
jag.sum

                mean         sd        2.5%         25%         50%          75%        97.5%     Rhat n.eff
alpha     19.6399640 1.22982919  17.1888046  18.8401629  19.6391880  20.45108889  22.03860486 1.001312  4100
beta1     -0.1063905 0.01755548  -0.1409861  -0.1180097  -0.1064871  -0.09469226  -0.07198614 1.001028 18000
beta2     -0.1857631 0.03667742  -0.2573464  -0.2103279  -0.1860142  -0.16155874  -0.11304929 1.001138  7700
deviance 192.6742580 3.02780455 188.9038177 190.4371154 192.0046398 194.13875526 200.28772929 1.001252  4800

write.table(x=jag.sum,file="out.txt",sep="\t")
于 2013-03-04T12:22:29.893 に答える