我正在MacOS上使用rjags包运行贝叶斯模型。
我试图运行返回的jags.model函数: “ jags.model(“ vinho.model.txt”中的错误,数据=列表(z = log.vinhos),n.chains = 3 ,: 运行时错误: 第4行出现编译错误。 z“的子集表达式中的维数不匹配
model{
for(t in 1:145){
a[t] <- mu[t] + gamma[t]
z[t] ~ dnorm(a[t], tau.y)
}
for(t in 143:145){
y[t] <- exp(z[t])
}
mu[1] ~ dnorm(0.0,0.001)
beta[1] ~ dnorm(0.0,0.001)
for(t in 2:145){
b[t] <- mu[t-1] + beta[t-1]
mu[t] ~ dnorm(b[t], tau.w1 )
beta[t] ~ dnorm(beta[t-1], tau.w2 )
}
for(t in 1:11){
gamma[t] ~ dnorm(0.0, 0.01)
}
for(t in 12:145){
c[t] <- -sum(gamma[(t-11):(t-1)])
gamma[t] ~ dnorm(c[t], tau.w3 )
}
tau.y ~ dgamma(0.01, 0.01)
tau.w1 ~ dgamma(0.01, 0.01)
tau.w2 ~ dgamma(0.01, 0.01)
tau.w3 ~ dgamma(0.01, 0.01)
V <- 1/tau.y
sig1.2 <- 1/tau.w1
sig2.2 <- 1/tau.w2
sig3.2 <- 1/tau.w3
}
因为我是JAGS的新手,所以我们将不胜感激!