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