求助,openBUGS在load data的时候报错expected a number or an NA,请问该怎么修改啊?机械专业的,完全不懂这个软件啊!谢谢啦
model{
C<-10000
for (i in 1: N){
zeros[i] <- 0
zeros[i] ~ dpois(phi[i])
}
for (j in 1: N) {
phi1[j]<-g*loga+d*logv+e*logf+log(p*lambda*beta*exp(-lambda*pow(time[j]-b,beta))*pow(time[j]-b,beta-1)-r*s*exp(-s*pow(time[j]-c,r))*pow(time[j]-c,r-1)*(p-1))+(pow(a,g)*pow(v,d)*pow(f,e)-1)*log(p*exp(-lambda*pow(time[j]-b,beta))-exp(-s*pow(time[j]-c,r))*(p-1)+1)
phi[j]<--phi1[j]+C
}
beta~ dunif(1, 2.5)
r~ dunif(1, 2.5)
lambda~ dgamma(0.001, 0.001)
s~ dgamma(0.001, 0.001)
b~ dgamma(0.001, 0.001)
c~ dgamma(0.001, 0.001)
d~ dgamma(0.001, 0.001)
e~ dgamma(0.001, 0.001)
g~ dgamma(0.001, 0.001)
p~ dgamma(0.001, 0.001)
}
list(N=24, time=c(119, 13, 20, 30, 36, 65, 37, 79, 79, 7, 50, 35, 4, 140, 9, 96, 24, 55, 4.5, 119, 8, 65, 31, 48),
v=c(72, 145, 145, 72, 104, 104, 145, 72, 72, 145, 104, 104, 206, 50, 104, 104, 104, 104, 206, 50, 104, 104, 104, 104),
f=c(0.12, 0.3, 0.12, 0.3, 0.2, 0.2, 0.12, 0.3, 0.12, 0.3, 0.2, 0.2, 0.2, 0.2, 0.6, 0.06, 0.2, 0.2, 0.2, 0.2, 0.6, 0.06, 0.2, 0.2),
a=c(2, 2, 2.5, 2.5, 2.25, 2.25, 2, 2, 2.5, 2.5, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 3, 1.5, 2.25, 2.25, 2.25, 2.25, 3, 1.5))
list(beta= 1.2235, lambda= 0.0079, r= 1.2235, s= 0.0079, b=0.5, c= 0.9573, d= 2.2899, e= 0.9573, g=1.3135, p= 0.4)