3

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