我想对RJAGS中处理单元故障所对应的计数分布进行共建模。处理单元的故障时间是威布尔分布的。受试者停留治疗的总时间段(总时间)也是威布尔分布的。两个威布尔互相耦合。
我知道R dweib参数化中两个Weibull分布的参数。我已将其转换为数据中的RJAGS参数化。
我知道Weibull的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)
我想知道将计数分布表示为累积危险的泊松分布的基本思想在理论上是否正确?