1

JAGS を使用する場合、次の形式のモデルから出力を受け取る方法は次のとおりです。

Inference for Bugs model at "model.txt", fit using jags,
3 chains, each with 10000 iterations (first 5000 discarded)
n.sims = 15000 iterations saved
mu.vect sd.vect 2.5% 25% 50% 75% 97.5% Rhat n.eff
mu 9.950 0.288 9.390 9.755 9.951 10.146 10.505 1.001 11000
sd.obs 3.545 0.228 3.170 3.401 3.534 3.675 3.978 1.001 13000
deviance 820.611 3.460 818.595 819.132 819.961 821.366 825.871 1.001 15000

バグと同様に、モデルが完了すると表示されると思いましたが、次の形式でしか取得できません。

Compiling model graph
   Resolving undeclared variables
   Allocating nodes
Graph information:
   Observed stochastic nodes: 1785
   Unobserved stochastic nodes: 1843
   Total graph size: 61542

Initializing model

  |++++++++++++++++++++++++++++++++++++++++++++++++++| 100%

基本的な質問で申し訳ありません。誰かが有用な JAGS 入門資料を提供できれば、それもまた有用です。

敬具。

4

1 に答える 1

1

「プラス」記号しか表示されない場合は、モデルを初期化しただけであることを意味します。ジャグが実際に実行されると、通常、後に「***」記号が生成されます。したがって、ここに行がありません (コードを確認できればよかったのに)。たとえば、r2jags を使用する場合は、次のように記述します。

out <- jags(data = data, parameters.to.save = params, n.chains = 3, n.iter = 90000,n.burnin = 5000,
                      model.file = modFile)
out.upd <- update(abundance.out.mod, n.iter=10000)
于 2016-10-13T14:05:57.227 に答える