[英]Error in jags.model RUNTIME ERROR: Compilation error on line 5. Index out of range taking subset of a
我对 Jags 很陌生,我正在尝试计算 model DIC,但返回了这样的错误:jags.model 中的错误(model.file,data = data,inits = init.values,n.chains = n.chains , : RUNTIME ERROR: Compilation error on line 5. 索引超出范围,我的代码的子集附在下面。
prior <- c(0.1, 0.15, 0.20, 0.25)
target <- 0.2
y<-c(0,0,0,4)
n<-c(0,5,5,8)
ndose=length(n)
bugs<-function(){
for (j in 1: ndose){
y[j]~dbin(p[j],n[j])
p[j]<-exp(3+a[j]*log(prior[j]/(1-prior[j])))/(1+exp(3+a[j]*log(prior[j]/(1-prior[j]))))
}
a~dgamma(1,1)
}
inits<-function(){
list('a'=0.001)
}
parameter<-c('a')
bugs.data<-list(y=y,n=n,prior=prior,ndose=ndose)
fit<-jags(data=bugs.data,inits=inits,parameter,bugs)
给定您提供的先验,您已将a
指定为标量,但在上面您已将其索引为1:ndose
。 如果a
确实是标量,则将代码更改为:
prior <- c(0.1, 0.15, 0.20, 0.25)
target <- 0.2
y<-c(0,0,0,4)
n<-c(0,5,5,8)
ndose=length(n)
bugs<-function(){
for (j in 1: ndose){
y[j]~dbin(p[j],n[j])
p[j]<-exp(3+a*log(prior[j]/(1-prior[j])))/(1+exp(3+a*log(prior[j]/(1-prior[j]))))
}
a~dgamma(1,1)
}
inits<-function(){
list('a'=0.001)
}
如果它是一个向量,那么你需要a
具有ndose
值。 仅使用上面的代码,可以这样完成:
target <- 0.2
y<-c(0,0,0,4)
n<-c(0,5,5,8)
ndose=length(n)
bugs<-function(){
for (j in 1: ndose){
y[j]~dbin(p[j],n[j])
p[j]<-exp(3+a[j]*log(prior[j]/(1-prior[j])))/(1+exp(3+a[j]*log(prior[j]/(1-prior[j]))))
}
for(j in 1:ndose){
a[j]~dgamma(1,1)
}}
inits<-function(){
list('a'=c(0.001, 1, 0.5, 0.7))
}
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.