全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
3590 8
2011-06-20
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))
错在什么地方?
附件列表

666.rar

大小:870 Bytes

 马上下载

本附件包括:

  • 666.odc

二维码

扫码加我 拉你入群

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

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

全部回复
2011-6-22 08:10:51
for (j in 2:n),需要定义j=1时的node吧
二维码

扫码加我 拉你入群

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

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

2011-6-22 08:36:48
按照提示的错误,for (j in 2:n),
log(beta[j]).....,index 应该是一个二维矩阵
二维码

扫码加我 拉你入群

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

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

2011-6-22 08:50:15
麻烦3楼的朋友能把改好的程序传上来吗?谢谢!
二维码

扫码加我 拉你入群

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

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

2011-6-22 09:20:42
你的第一循环中的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)
}
二维码

扫码加我 拉你入群

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

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

2011-6-22 15:38:05
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
希望楼上的朋友看看!
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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