全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件
1502 2
2011-12-22
向高手请教
用winbugs 做死亡率修云的时候,需要假设死亡率随年龄递增
如以theta 表示死亡率,则需限制条件theta[1]<theta[2]<theta[3]<  
以下是自编代码
model{
dx[1]~dpois(lambda[1]);
lambda[1]<-ex[1]*theta[1];
theta[1]~dgamma(1.49,beta)I(0,theta[2]);
for(i in 2:29)
{
   dx~dpois(lambda);
   lambda<-ex*theta;
   theta~dgamma(1.49,beta)I(theta[i-1],theta[i+1]);
   }
dx[30]~dpois(lambda[30]);
  lambda[30]<-ex[30]*theta[30];
  theta[30]~dgamma(1.49,beta)I(theta[29],0.02);

  beta<-pow(gamma,-1);
  gamma~dgamma(3.0,115);
}
但运行时候出现
“cannot sample from interval censored gamma full conditonal theta[1]"
我感觉是I(theta[i-1],theta[i+1])的问题
请问高手这种情况怎么处理,如何用winbugs 实现这种限制条件

非常感谢
二维码

扫码加我 拉你入群

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

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

全部回复
2011-12-24 23:08:27
哈哈 楼主好牛啊
二维码

扫码加我 拉你入群

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

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

2013-2-25 14:11:19
我也发生了类似的问题,请问楼主怎么解决的
二维码

扫码加我 拉你入群

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

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

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

分享

扫码加好友,拉您进群