如何解决锯齿中的尺寸问题

时间:2019-07-15 01:15:15

标签: r jags rjags

我正在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的新手,所以我们将不胜感激!

0 个答案:

没有答案