全部版块 我的主页
论坛 提问 悬赏 求职 新闻 读书 功能一区 经管百科 爱问频道
1261 3
2016-11-24
求助,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)

二维码

扫码加我 拉你入群

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

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

全部回复
2016-11-24 16:11:52
谢谢啦
二维码

扫码加我 拉你入群

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

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

2016-11-24 20:00:39
求大神指导啊,谢谢
二维码

扫码加我 拉你入群

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

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

2016-11-25 15:33:19
求帮助啊
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群