winBUGS 泄漏频率模型 - 脚本错误

问题描述

我正在对两个组件中的泄漏频率分布进行建模。该模型假设任何组件的平均泄漏频率 (mu) 与泄漏的分数大小的对数线性相关。假设线性关系的系数 (alpha1,alpha2) 服从正态分布。

假设泄漏频率本身与参数 (mu,tau) 呈对数正态分布，tau 是由 gamma 描述的精度,和 mu 假设泊松过程在定义的时间内记录。 winBUGS 中的模型和加载数据是成功的，但是在初始化 alpha1,alpha2,tau 时它返回“此链包含未初始化的变量”。 如果我遗漏了脚本中的某些部分或无法成功初始化的原因，任何人都可以提供帮助吗？

model 
{
  time.yr <- time/8760                       #Convert hours to years  
  for(i in 1:nitems) {                       #Components
      for(j in 1:nleaks) {
          
          k[i,j] ~ dpois(mean.poisson[i,j])          #Poisson distr. for 'k' events in time
          mean.poisson[i,j] <- lambda[i,j]*time.yr   #Poisson parameter

          lambda[i,j] ~ dlnorm(mu[i,j],tau[i])      #Lognormal prior for lambda
 
        mu[i,j] <- pow(10,alpha1[i])*pow(fr.leak[i,alpha2[i])   #model of leak freq. 

     }
     alpha1[i] ~ dnorm(0,1.0E-3)      #Diffuse priors on alpha1,alpha2 and tau
     alpha2[i] ~ dnorm(0,1.0E-3) 
     tau[i] ~ dgamma(1,1)                         
   }
}


data
list(
      nleaks = 4,nitems = 2,time = 591300,#constants
      fr.leak = structure(                                                  
                      .Data = c(0.01E-2,0.1E-2,1.0E-2,10.0E-2,#recorded leak size 
                                0.01E-2,10.0E-2),.Dim = c(2,4)),#input structure of (nitems x nleak)=8 observations                      
              k = structure(                                                                                     
                             .Data = c(8,8,1,#comp.1 leak records
                                       2,2,0),#comp.2 leak records                          
                             .Dim = c(2,4))
)

inits
list(alpha1=c(0.1,0.1),alpha2=c(0.01,0.01),tau=c(1.0E-3,1.0E-3))
list(alpha1=c(0.1,alpha2=c(-0.01,-0.01),1.0E-3))

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）