Winbugs中的一个简单模型,但它说“此链包含未初始化的变量”

时间:2016-07-29 18:19:50

标签: winbugs weibull

我有一些简单的时间来处理事件数据,没有协变量。我试图将Weibull分布拟合到它。所以我有以下代码。在我加载首字母之前,一切看起来都很好。它说“这个链包含未初始化的变量”。但我不明白。我认为Weibull dist只有2个参数,我已经指定了它们。您能否提一些建议?谢谢!

model
{   
    for(i in 1 : N) {                       
            t[i] ~ dweib(r, mu)I(t.cen[i],)
    }   
    mu ~ dexp(0.001)
    r ~ dexp(0.001)
}
# Data
list(
t.cen=c(0,3.91,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,21.95,23.98,33.08),
t=c(2.34,NA,5.16,5.63,6.17,6.8,7.03,8.05,8.13,8.36,8.83,10.16,
10.55,10.94,11.48,11.95,13.05,13.59,16.02,20.08,NA,NA,
NA),
N=23
)
# Initial values
list(
r=3,mu=3
)

1 个答案:

答案 0 :(得分:0)

其他未初始化的变量是t向量中的缺失(NA)值。请记住,BUGS语言不区分数据和参数,而提供值为NA的数据则相当于不将其作为数据提供。