为什么winbugs运行结果总是出现educational version can not do this model;
程序如下:model
{c<-10000
for(i in 1:N)
{zeros[i ]<-0
zeros[i ]~dpois(zeros.p[i ])
zeros.p[i ] <--l[i ]+c
l[i ]<--log(alpha)-(1+1/alpha)*log(1+alpha*x[ i ]/beta)
}
theta~dgamma(0.1,1.0)
beta<-1/theta
alpha~dpar(0.1,0.1)I(0.0001,10)
}
list(x=c(0.10,0.10,0.10,0.10,0.10,0.10,0.10,0.10,0.10,0.10,0.10,0.10,0.10,0.10,
0.10,0.10,0.10,0.10,0.10,2.10,8.10,10.10,10.10,10.10,10.10,10.10,10.10,
10.20,20.10,20.10,20.10,20.10,20.10,20.10,20.10,20.10,20.10,20.10,20.10
,20.10,20.10,25.10,25.10,25.10,30.10,30.10,30.10,30.90,40.10,40.10,40.1
0,50.10,50.10,50.10,50.10,50.10,50.10,50.10,50.10,50.10,50.10,50.10,55.
10,60.10,70.10,70.10,70.10,70.10,78.10,78.10,78.10,80.10,80.10,80.10,80
.10,80.10,80.10,80.10,88.10,90.10,95.10,100.10,100.10,100.10,110.10,110
.10,120.10,120.10,140.10,200.10,250.10,250.10,300.10,400.10,480.10,500.
10,800.10),N=96)
list(alpha=0.1,theta=1)