Я пишу сценарий JAGS
(иерархическая байесовская модель), в котором времена событий моделируются как гонка между двумя процессами.
Наблюдения: time
- это измеренное время событий.
Модель: два процесса с гауссовыми скоростями - какой процесс завершается первым, запускаетсясобытие.
Цель: Оценить скорости двух процессов.
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!
}
}
Вопрос: Как я могу указать, что отсчет времени производится по большей из двух скоростей, каждая из которых отсчитывается по разному распределению?
Пример гистограмма смоделированных time
данных с использованием mu1=3, mu2=3
с различными стандартными отклонениями для двух процессов (фиксированные на 1 и 0,1)