0

R2Jags パッケージを使用しています。

最初に次のようにモデルを取得します。

jags.data = list("Y"=y, "N"=N, Y1=y[1], C=x1)
# The parameters that we are monitoring (must monitor at least 1)
jags.params=c("sd.q","sd.r","Y", "X1", "X2")
model.loc=("ss_model.txt")
mod_ss = jags(jags.data, parameters.to.save=jags.params, model.file=model.loc, n.chains = 3, 
              n.burnin=5000, n.thin=1, n.iter=10000, DIC=TRUE)  

しかし、次を使用して X1 の事後値を取得すると、次のようになります。

attach(mod_ss)
X1

私は得る:

エラー: オブジェクト 'X1' が見つかりません

私の質問: これを修正するにはどうすればよいですか?

前もって感謝します。

4

1 に答える 1

2

問題が解決しました:

R2Jags の新しいバージョンでは、呼び出す正しい関数は次のようです。

attach.jags(mod_ss, overwrite=TRUE)

[だからattach.jags代わりにattach()]

于 2015-04-11T18:51:26.897 に答える