我有一些简单的时间事件数据,没有协变量。我试图将 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
)
最佳答案
其他未初始化的变量是 t 向量中缺失的 (NA) 值。请记住,BUGS 语言在数据和参数之间没有区别,并且将值 NA 的某物作为数据提供等同于不将其作为数据提供。
关于winbugs - Winbugs 中的一个简单模型,但显示为 "This chain contains uninitialized variables",我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/38665210/