私はJAGSで彼のモデルコードを実行しようとしています:
model{
## priors
b0 ~ dnorm(0, 5)
b1 ~ dnorm(0, 5)
b2 ~ dnorm(0, 5)
b3 ~ dnorm(0, 5)
b4 ~ dnorm(0, 5)
b5 ~ dnorm(0, 5)
b6 ~ dnorm(0, 5)
sigmayear ~ dt(0, 1, 2) I(0, )
## likelihood
for(i in 1:nind){
for(t in 2:last[i]){
y[i,t] ~ dbern(prop[i, t - 1])
prop[i, t - 1] <- y[i, t - 1] * (1 - d[i, t - 1])
logit(d[i, t - 1]) <-
b0 + b1 * temp[i, t] + b2 * sex[i] + b3 * day[t] +
b4 * strat[i] + b5 * temp[i, t] * sex[i] +
b6 * temp[i, t] * strat[i] + sigmayear * ryear[years[i]]
} # close t
} # close i
for(k in 1:nyears){
ryear[k] ~ dnorm(0, 1)
} #close k
}
私が作成したデータ マトリックスは問題ないはずです。
> str(bugs.data)
List of 9
$ y : num [1:124, 1:203] 1 1 1 1 1 1 1 1 1 1 ...
$ nind : int 124
$ last : num [1:124] 98 98 78 91 78 98 87 88 98 111 ...
$ temp : num [1:124, 1:203] 1.33 1.33 1.33 1.33 1.33 ...
$ sex : num [1:124] 1 0 1 1 1 0 0 0 0 1 ...
$ day : num [1:203] -1.72 -1.7 -1.69 -1.67 -1.65 ...
$ strat : num [1:124] 0 0 0 0 0 0 0 0 0 0 ...
$ years : num [1:124] 7 7 7 7 7 7 7 7 7 1 ...
$ nyears: int 8
これらのパラメーターを使用した実際の呼び出しは、以前は機能していました...
ni <- 5000 #number of sims
nt <- 1 #get rid of
nb <- 2000 #warm in
nc <- 3 # number of chains
しかし、何らかの理由で、私がそれを呼び出したいとき、これが起こります:
> ms <- jags(bugs.data, inits, parameters, "knownfateSurvivalfat.txt", n.chains = nc, n.thin = nt, n.iter = ni, n.burnin = nb, working.directory = getwd())
Compiling model graph
Resolving undeclared variables
Deleting model
Fehler in jags.model(model.file, data = data, inits = init.values, n.chains = n.chains, :
RUNTIME ERROR:
Compilation error on line 14.
Index out of range taking subset of y
「y」が範囲外であるべき理由を誰かが考えていますか?
「i」は常にマトリックス (bugs.data) に適合する 1 から 124 までで、「t」は 2 から 203 までで、これも問題ありません。