我正在使用锯齿进行贝叶斯分析,但不幸的是,我遇到了一个无法确定其原因的错误。我已经尝试研究一个星期的其他示例,但是问题仍然存在。这是我得到的错误: 运行时错误: 无法将节点插入mu [1]。尺寸不匹配。 下面是我的代码,这个问题主要是在我命令最后一行时出现的。请帮助我,因为我现在已经呆了一个星期。
mod.jags<-paste("
model{
for (i in 1:N){
perinatal[i]~dbern(p[i])
logit(p[i])<-mu[i]
mu[i]<-beta[1]+beta[2]*gender[i]+beta[3]*Bord[i]+beta[4]*Marital[i]+beta[5]*mweight[i]+beta[6]*residence[i]
+beta[7]*Birthint[i]+beta[8]*BMI+beta[9]*age[i]+beta[10]*Occupation[i]+beta[11]*SEWI[i]+beta[12]*age_firstbirth[i]
+beta[13]*Mateducat[i]
}
for (j in 1:13){
beta[j]~dbeta(2,2)
}
}
")
writeLines(mod.jags,"mod.jags")
##Parameters for posterior
params.jags<-c("beta")
inits1.jags<-list("beta"=rep(0,13))
inits.jags<-function(){
inits=list("beta"=0.0)
}
##Fitting the model using R2jags
set.seed(50)
fit.jags<-jags(data = dat1,inits = inits.jags,parameters.to.save = params.jags,n.chains = 4,n.iter = 2000,n.burnin = 1000,
model.file="mod.jags")