Я хочу совместно смоделировать распределение отсчетов, соответствующее отказу лечебных единиц в RJAGS.Время отказов для лечебного учреждения распределено по Вейбуллу.Общий период времени, в течение которого субъект остается на лечении (общее время), также распределяется по Вейбуллу.Два Вейбулла связаны друг с другом.
Я знаю параметры обоих распределений Вейбулла в параметризации R dweib.Я преобразовал это в параметризацию RJAGS в данных.
Я знаю о различиях в параметризации R и RJAGS Вейбулла.https://stats.stackexchange.com/questions/18550/how-do-i-parameterize-a-weibull-distribution-in-jags-bugs.
Я думал о следующей стратегии:
# Survival of individual treatment units
survivaltime ~ dweib(shape, scale)
# Survival of total treatment time
totaltime ~ dweib(shape1,scale1)
# MODEL FOR COUNTS. Not sure about this!!!
counts ~ dpois(cumulativehazard)
# Cumulative hazard of the Weibull is time^a
cumulativehazard<-pow(totaltime/shape,a)
# I know that the Weibull scale a in R is different from RJAGS parameterization
a<-pow(scale,-1/shape)
Мне было интересно, является ли теоретически правильной основная идея выражения распределения счета как распределения Пуассона накопленной опасности?