JAGSイベントの時間が 2 つのプロセス間の競争としてモデル化されるスクリプト (階層ベイジアン モデル)を書いています。
観察: timeイベントの測定時間です。
モデル:ガウス レートの 2 つのプロセス - どちらのプロセスが先に終了しても、イベントがトリガーされます。
目標: 2 つのプロセスの速度を推定します。
model{
# Priors
mu1 ~ dnorm( 0,1 ) # rate of one process
mu2 ~ dnorm( 0,1 ) # rate of other process
sigma1 <- 1 # variability in rate
sigma2 <- 0.1 # variability in rate
# Observations
for (i in 1:N)
rate1[i] ~ dnorm( mu1, sigma1 ) # Sample the two
rate2[i] ~ dnorm( mu2, sigma2 ) # racing processes.
rmax[i] <- max( rate1[i], rate2[i] ) # which was faster?
time[i] ~ 1/rmax[i] #### This is wrong!
}
}
質問:それぞれが異なる分布からサンプリングされた 2 つのレートのうち大きい方から時間がサンプリングされていることをどのように示すことができますか?
2 つのプロセスの異なる標準偏差 (1 と 0.1 に固定)timeを使用した、シミュレートされたデータのヒストグラムの例mu1=3, mu2=3
