全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 winbugs及其他软件专版
1381 0
2014-04-13
为什么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)
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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