私は最近、JAGS の使用を開始し、R 内でそれを呼び出しました。最終的に、コードを使用して R にリンクされたジャグを取得しました。
install.packages("rjags")
library(rjags)
そして出力を得ました
Linked to JAGS 3.4.0
Loaded modules: basemod,bugs
JAGSモデルのデータも別ファイルにBUG形式で保存しました(教えてもらった通り)。
ただし、データを実行しようとすると、エラーメッセージが表示され続けます。
Error in file(modfile, "rt") : cannot open the connection
In addition: Warning message:
In file(modfile, "rt") :
cannot open file 'age_problem.bug': No such file or directory
Error in jags.model("age_problem.bug", data = list(X = X, N = length(X)), :
Cannot open model file "age_problem.bug"
と
Error in update(jags, 1000) : object 'jags' not found
私が見逃している重要なステップはありますか?
編集:問題の例のコード
N <- 1000
x <- rnorm(N, 0, 5)
write.table(x,
file = 'example1.data',
row.names = FALSE,
col.names = FALSE)
library('rjags')
jags <- jags.model('example1.bug',
data = list('x' = x,
'N' = N),
n.chains = 4,
n.adapt = 100)
update(jags, 1000)
jags.samples(jags,
c('mu', 'tau'),
1000)
JAGS モデル:
model {for (i in 1:N) {
x[i] ~ dnorm(mu, tau)}
mu ~ dnorm(0, .0001)
tau <- pow(sigma, -2)
sigma ~ dunif(0, 100)}