CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDescxDocuments.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc TextViews.StdViewDescTextViews.ViewDeschTextModels.StdModelDescTextModels.ModelDescr j QTextModels.AttributesDesc* # HERE IS THE WEIBULL REGRESSION MODEL, CARLIN AND HODGES (1999) DATA: model { for(i in 1:nsub[1]) { t1[i] ~ dweib(r[1], mu1[i])I(t1.cen[i],) X1[i] <- 2*drug1[i] - 3 mu1[i] <- exp(beta0+beta1*X1[i]+W[1]) } for(i in 1:nsub[2]) { t2[i] ~ dweib(r[2], mu2[i])I(t2.cen[i],) X2[i] <- 2*drug2[i] - 3 mu2[i] <- exp(beta0+beta1*X2[i]+W[2]) } for(i in 1:nsub[3]) { t3[i] ~ dweib(r[3], mu3[i])I(t3.cen[i],) X3[i] <- 2*drug3[i] - 3 mu3[i] <- exp(beta0+beta1*X3[i]+W[3]) } for(i in 1:nsub[4]) { t4[i] ~ dweib(r[4], mu4[i])I(t4.cen[i],) X4[i] <- 2*drug4[i] - 3 mu4[i] <- exp(beta0+beta1*X4[i]+W[4]) } for(i in 1:nsub[5]) { t5[i] ~ dweib(r[5], mu5[i])I(t5.cen[i],) X5[i] <- 2*drug5[i] - 3 mu5[i] <- exp(beta0+beta1*X5[i]+W[5]) } for(i in 1:nsub[6]) { t6[i] ~ dweib(r[6], mu6[i])I(t6.cen[i],) X6[i] <- 2*drug6[i] - 3 mu6[i] <- exp(beta0+beta1*X6[i]+W[6]) } for(i in 1:nsub[7]) { t7[i] ~ dweib(r[7], mu7[i])I(t7.cen[i],) X7[i] <- 2*drug7[i] - 3 mu7[i] <- exp(beta0+beta1*X7[i]+W[7]) } for(i in 1:nsub[8]) { t8[i] ~ dweib(r[8], mu8[i])I(t8.cen[i],) X8[i] <- 2*drug8[i] - 3 mu8[i] <- exp(beta0+beta1*X8[i]+W[8]) } for(i in 1:nsub[9]) { t9[i] ~ dweib(r[9], mu9[i])I(t9.cen[i],) X9[i] <- 2*drug9[i] - 3 mu9[i] <- exp(beta0+beta1*X9[i]+W[9]) } for(i in 1:nsub[10]) { t10[i] ~ dweib(r[10], mu10[i])I(t10.cen[i],) X10[i] <- 2*drug10[i] - 3 mu10[i] <- exp(beta0+beta1*X10[i]+W[10]) } for(i in 1:nsub[11]) { t11[i] ~ dweib(r[11], mu11[i])I(t11.cen[i],) X11[i] <- 2*drug11[i] - 3 mu11[i] <- exp(beta0+beta1*X11[i]+W[11]) } for (i in 1:NSITES){ b0star[i] <- beta0 + W[i] r[i] ~ dgamma(alpha, alpha) W[i] ~ dnorm(0, tau) } alpha <- 10 tau ~ dgamma(1, 1); sigmasqu <- 1 / tau; beta0 ~ dnorm(0.0, 0.0001) beta1 ~ dnorm(0.0, 0.0001) RR <- exp(2*beta1) } # END OF C-H FRAILTY MODEL # Initial values: list(r=c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), tau=1.0, beta1=1.0, beta0=-10.0, W=c(0,0,0,0,0,0,0,0,0,0,0)) # Here are the survival data: list(NSITES=11, nsub=c(4,2,1,10,8,12,6,6,13,4,3), t1=c(NA,248,NA,344), t1.cen=c(74,0,272,0), t2=c(NA,NA), t2.cen=c(4,156), t3=c(NA), t3.cen=c(100), t4=c(NA,64,88,NA,NA,NA,NA,NA,NA,NA), t4.cen=c(20,0,0,148,162,184,188,198,382,436), t5=c(NA,NA,82,NA,NA,214,NA,262), t5.cen=c(50,64,0,186,214,0,228,0), t6=c(6,NA,76,80,202,NA,NA,NA,NA,NA,NA,NA), t6.cen=c(0,16,0,0,0,258,268,368,380,424,428,436), t7=c(NA,NA,102,NA,NA,NA), t7.cen=c(32,64,0,162,182,364), t8=c(NA,NA,NA,NA,NA,162), t8.cen=c(22,22,74,88,148,0), t9=c(8,NA,40,NA,NA,NA,NA,276,NA,366,NA,NA,NA), t9.cen=c(0,16,0,120,168,174,268,0,286,0,396,466,468), t10=c(NA,NA,NA,254), t10.cen=c(18,36,160,0), t11=c(NA,NA,NA), t11.cen=c(28,70,106), drug1=c(1,2,1,2), drug2=c(2,1), drug3=c(2), drug4=c(2,2,2,2,1,1,1,1,1,1), drug5=c(1,2,2,1,1,1,2,2), drug6=c(1,2,1,2,2,1,1,2,1,1,2,2), drug7=c(2,1,1,2,2,1), drug8=c(2,1,1,1,1,2), drug9=c(2,2,2,1,1,2,1,2,1,1,2,2,1), drug10=c(1,1,2,2), drug11=c(1,1,2)) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$0ZGo * <[ @Documents.ControllerDesc v x k{