winbugs第二个程序有错求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

model{
    c<-10
    for(i in 1:n){zeros<-0
    zeros~dpois(zeros.mean)
    zeros.mean<--log(l)+c  
    l<-alpha*beta*(pow(y,(beta-1)))*(exp(-alpha*((pow(y,beta)-(pow(t,beta))))))}
    for(j in 2:n){log(beta[j])<-log(beta[j-1])+(eta[j])
    eta[j]~dnorm(0,lamda)}     
    s<-sqrt(1/lamda)
    lamda~dgamma(1,0.5)
    alpha~dgamma(1,0.5)
    beta[1]~dlnorm(0,0.1)
}
#Data    list(n=19,y=c(0.3945,0.7435,0.9055,1.1075,1.4631,1.8113,1.9136,2.1390,2.3514,2.7121,2.7612,2.7807,2.8761,2.9606,2.9802,3.0674,3.1310,3.1719,3.2531))
#Initial values         
list(alpha=44,mu=0.01,tau=0.01,beta=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1))
错在什么地方？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS bug Win 求助 程序 winbugs

for (j in 2:n),需要定义j=1时的node吧

按照提示的错误，for (j in 2:n),
log（beta[j]).....,index 应该是一个二维矩阵

麻烦3楼的朋友能把改好的程序传上来吗？谢谢！

你的第一循环中的beta是一个系数，怎么在第二个循环中的beta成了beta[j],他们俩个是一个参数？把第二个循环中的beta换成b，就好了，但是我不知道你的模型是什么道理
for(j in 2:n){

log(b[j])<-log(b[j-1])+(eta[j])
    eta[j]~dnorm(0,lamda)
}

model{
    c<-10
    for(i in 1:n){zeros<-0
    zeros~dpois(zeros.mean)
    zeros.mean<--log(l)+c  
    l<-alpha*beta*(pow(y,(beta-1)))*(exp(-alpha*((pow(y,beta)-(pow(t,beta))))))}
   for(j in 2:n){

log(b[j])<-log(b[j-1])+(eta[j])
    eta[j]~dnorm(0,lamda)
}
    s<-sqrt(1/lamda)
    lamda~dgamma(1,0.5)
    alpha~dgamma(1,0.5)
    beta~dlnorm(0,0.1)
}
#Data    list(n=19,y=c(0.3945,0.7435,0.9055,1.1075,1.4631,1.8113,1.9136,2.1390,2.3514,2.7121,2.7612,2.7807,2.8761,2.9606,2.9802,3.0674,3.1310,3.1719,3.2531))
#Initial values         
list(alpha=44,mu=0.01,tau=0.01,beta=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1))

现在模型check正确，但在load data时，提示以下错误：
expected a number or an NA
希望楼上的朋友看看！

data里应该给b值吧

model{
    c<-10
    for(i in 1:n){zeros<-0
    zeros~dpois(zeros.mean)
    zeros.mean<--log(l)+c  
    l<-alpha*beta*(pow(y,(beta-1)))*(exp(-alpha*((pow(y,beta)-(pow(t,beta))))))}
   for(j in 2:n){

log(b[j])<-log(b[j-1])+(eta[j])
    eta[j]~dnorm(0,lamda)
}
    s<-sqrt(1/lamda)
    lamda~dgamma(1,0.5)
    alpha~dgamma(1,0.5)
    beta~dlnorm(0,0.1)
}
#Data    list(n=19,y=c(0.3945,0.7435,0.9055,1.1075,1.4631,1.8113,1.9136,2.1390,2.3514,2.7121,2.7612,2.7807,2.8761,2.9606,2.9802,3.0674,3.1310,3.1719,3.2531),b=c(0.3212,0.7232,0.9215,1.1075,1.4231,1.2313))
#Initial values         
list(alpha=44,mu=0.01,tau=0.01,beta=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1))
在data里给b赋值后，错误提示还是那样？如何解决？谢谢！

您好，请问出现报错expected a number or an NA是什么原因啊