3

我有一些简单的时间来事件数据,没有协变量。我试图将 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
)
4

1 回答 1

2

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

于 2016-07-30T06:22:01.357 回答