openBUGS求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求助，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)

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：openbugs BUGS Open bug Pen expected number 机械 软件 专业

谢谢啦

求大神指导啊，谢谢

求帮助啊