急 请教一个winbugs 问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

向高手请教
用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 实现这种限制条件

非常感谢

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS bug Win 死亡率

哈哈 楼主好牛啊

我也发生了类似的问题，请问楼主怎么解决的