イベントデータまでの時間は簡単で、共変量はありません。私はそれにワイブル分布を当てはめようとしていました。だから私は次のコードを持っています。イニシャルをロードするまでは、すべて問題ないようです。「このチェーンには初期化されていない変数が含まれています」と表示されます。しかし、私は理解していません。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
)