Winbugs程序问题

ym2009

2350

收藏 2013-10-23

悬赏 10 个论坛币已解决

model{
d<-100000
for(i in 1:n){
zeros<-0
zeros~dpois(zeros.mean)
zeros.mean<--log(l)+d
l<-theta*theta*(1+x)*exp(-theta*x)/(1+theta)
}
theta~dgamma(0.5,5.0E-04)
}
list(n=30,x=c(0.8,2.9,4.3,5.0,6.7,8.2,9.7,11.9,14.1,19.9,0.8,3.1,4.3,5.3,6.9,8.6,9.8,12.4,15.4,20.6,1.3,3.2,4.2,5.5,7.1,8.6,10.7,12.5,15.4,21.3))
list(theta=1)
帮我看看这程序，怎么运行不了，我做的是Lindley分布的贝叶斯估计

最佳答案

TimeT 查看完整内容

我改了一下，现可以运行了。WINBUGS运行结果如下： node mean sd MC error 2.5% median 97.5% start sample theta 0.212 0.02732 1.422E-4 0.1618 0.2109 0.2684 10001 40000 我改过的程序如下： model{ d

扫码加我拉你入群

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

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

全部回复

TimeT

2013-10-23 12:04:04

我改了一下，现可以运行了。WINBUGS运行结果如下：
node mean sd MC error 2.5% median 97.5% start sample
theta 0.212 0.02732 1.422E-4 0.1618 0.2109 0.2684 10001 40000

我改过的程序如下：
model{
d<-100000
for( k in 1 : N ) {
zeros[k]<-0
phi[k]<- d-log(L[k])
zeros[k]~dpois(phi[k])
L[k]<- theta*theta/(1+theta)*(1+x[k])*exp(-theta*x[k])
}
theta~dgamma(0.5,5.0E-04)
}

list中的数据就不抄一遍了。

你运行一下试试。