Winbugs程序问题

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分布的贝叶斯估计

我改了一下，现可以运行了。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中的数据就不抄一遍了。

你运行一下试试。

d<-100000
    for(i in 1:n){

这里的n是否要改成d？

TimeT 发表于 2013-10-23 12:04
我改了一下，现可以运行了。WINBUGS运行结果如下：
node         mean         sd         MC error        2.5%        median        97.5%        start        sa ...

可以运行了，产生一个新的分布我的程序中出问题了。

我想用winbugs解决一个简单的Bayes求解过程，已经知道先验分布和样本值，怎么编这段程序啊？我是初学者期待您的回复，谢谢

TimeT 发表于 2013-10-23 12:04
我改了一下，现可以运行了。WINBUGS运行结果如下：
node         mean         sd         MC error        2.5%        median        97.5%        start        sa ...

Could you give some reason for this change? Thanks

Lisrelchen 发表于 2013-12-7 22:21
Could you give some reason for this change? Thanks

WINBUGS的manual中有一节“Specifying a new sampling distribution”，那里介绍了“ZEROS TRICK”，看了那里，我才这么改的。

nuomin 发表于 2013-10-23 12:38
d

d与N不同，不应改。你可以看看WINBUGS的manual中有一节“Specifying a new sampling distribution”，那里介绍了“ZEROS TRICK”，那里C就是这里的d，那里的N就是这里的N。