我有一些简单的时间来处理事件数据,没有协变量。我试图将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
)
答案 0 :(得分:0)
其他未初始化的变量是t向量中的缺失(NA)值。请记住,BUGS语言不区分数据和参数,而提供值为NA的数据则相当于不将其作为数据提供。