CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesck+Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc* *TextViews.StdViewDescTextViews.ViewDesc[*TextModels.StdModelDescTextModels.ModelDesce(](TextModels.AttributesDesc'* *L CD 6  J # # models 1-4 (NORMAL MODEL) model{ for (i in 1:N){ # -------------------- # model for oed # -------------------- oed[i]~dnorm( mu.oed[i], tau.oed[1]) mu.oed[i] <- b.oed[1] + b.oed[2] *(period[i]-1)+ b.oed[3] *(treatment[i]-1) + g[1]*a.oed[ patient[i] ] + (1-g[1])*a.oed.fixed[ patient[i] ] + g[2]*b.oed[4]*oedbase[i] # -------------------- # model for dbp # -------------------- dbp[i]~dnorm( mu.dbp[i], tau.dbp[1]) mu.dbp[i] <- b.dbp[1] + b.dbp[2] *(period[i]-1)+ b.dbp[3] *(treatment[i]-1) + g[1]*a.dbp[ patient[i] ]+(1-g[1])*a.dbp.fixed[ patient[i] ] + g[2]*b.dbp[4]*dbpbase[i] } for (i in 1:n){ # # Hyper priors for individual/patients random effects a.oed[i]~dnorm( 0.0, tau.oed[2]) a.dbp[i]~dnorm( 0.0, tau.dbp[2]) # # Non-informative priors for individual/patients fixed effects a.oed.fixed[i]~dnorm( 0.0, 0.001) a.dbp.fixed[i]~dnorm( 0.0, 0.001) } for (i in 1:p){ b.oed[i]~dnorm( 0.0, 0.001) b.dbp[i]~dnorm( 0.0, 0.001) } tau.oed[1]~dgamma( 0.001,0.001) tau.oed[2]~dgamma( a, a) tau.dbp[1]~dgamma( 0.001,0.001) tau.dbp[2]~dgamma( a, a) a <- g[1]*0.001 + (1-g[1])*1 # # s2[1]<-1/tau.oed[1] s2[2]<-1/tau.oed[2] s2[3]<-1/tau.dbp[1] s2[4]<-1/tau.dbp[2] for( i in 1:N ){ res1[i] <- oed[i] - mu.oed[i] res2[i] <- dbp[i] - mu.dbp[i] } R[1] <- 1 - pow( sd(res1[1:N])/sd(oed[1:N]), 2) R[2] <- 1 - pow( sd(res2[1:N])/sd(dbp[1:N]), 2) } INITS list( b.oed=c(0,0,0,0), b.dbp=c(0,0,0,0), tau.oed=c(1,1), tau.dbp=c(1,1)) DATA list(N=188, n=94, p=4, g=c(1,1), patient = c(1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8, 8, 9, 9, 10, 10, 11, 11, 12, 12, 13, 13, 14, 14, 15, 15, 16, 16, 17, 17, 18, 18, 19, 19, 20, 20, 21, 21, 22, 22, 23, 23, 24, 24, 25, 25, 26, 26, 27, 27, 28, 28, 29, 29, 30, 30, 31, 31, 32, 32, 33, 33, 34, 34, 35, 35, 36, 36, 37, 37, 38, 38, 39, 39, 40, 40, 41, 41, 42, 42, 43, 43, 44, 44, 45, 45, 46, 46, 47, 47, 48, 48, 49, 49, 50, 50, 51, 51, 52, 52, 53, 53, 54, 54, 55, 55, 56, 56, 57, 57, 58, 58, 59, 59, 60, 60, 61, 61, 62, 62, 63, 63, 64, 64, 65, 65, 66, 66, 67, 67, 68, 68, 69, 69, 70, 70, 71, 71, 72, 72, 73, 73, 74, 74, 75, 75, 76, 76, 77, 77, 78, 78, 79, 79, 80, 80, 81, 81, 82, 82, 83, 83, 84, 84, 85, 85, 86, 86, 87, 87, 88, 88, 89, 89, 90, 90, 91, 91, 92, 92, 93, 93, 94, 94), treatment = c(2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 1, 2, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 1, 2, 1, 2, 2, 1, 2, 1, 2, 1, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 1, 2, 2, 1, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 1, 2, 2, 1, 2, 1, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 2, 1, 2, 1), period = c(1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2), oedbase = c(45, 45, 51, 51, 53, 53, 49, 49, 46, 46, 61, 61, 55, 55, 52, 52, 54, 54, 50, 50, 58, 58, 52, 52, 62, 62, 54, 54, 56, 56, 62, 62, 59, 59, 50, 50, 53, 53, 49, 49, 53, 53, 59, 59, 47, 47, 39, 39, 53, 53, 68, 68, 64, 64, 59, 59, 68, 68, 64, 64, 68, 68, 70, 70, 68, 68, 68, 68, 66, 66, 63, 63, 56.5, 56.5, 59, 59, 55, 55, 57, 57, 66, 66, 47, 47, 64, 64, 76, 76, 89, 89, 69, 69, 72, 72, 43, 43, 52, 52, 49, 49, 54, 54, 64, 64, 56, 56, 58, 58, 54, 54, 49, 49, 49, 49, 45, 45, 53, 53, 48, 48, 50.5, 50.5, 50, 50, 50, 50, 56, 56, 48, 48, 45, 45, 53, 53, 55, 55, 51, 51, 47, 47, 53, 53, 53, 53, 61, 61, 57, 57, 59, 59, 55, 55, 53, 53, 52, 52, 54, 54, 59, 59, 53, 53, 57, 57, 60, 60, 58, 58, 56, 56, 49, 49, 52, 52, 60, 60, 45, 45, 57, 57, 63, 63, 60, 60, 54, 54, 73, 73), dbpbase = c(60, 60, 50, 50, 70, 70, 68.3333333333333, 68.3333333333333, 65, 65, 95.3333333333333, 95.3333333333333, 84.6666666666667, 84.6666666666667, 85.3333333333333, 85.3333333333333, 94.6666666666667, 94.6666666666667, 89.3333333333333, 89.3333333333333, 96.6666666666667, 96.6666666666667, 88.6666666666667, 88.6666666666667, 94.6666666666667, 94.6666666666667, 96, 96, 91.3333333333333, 91.3333333333333, 99.3333333333333, 99.3333333333333, 98.6666666666667, 98.6666666666667, 70.6666666666667, 70.6666666666667, 89.3333333333333, 89.3333333333333, 80, 80, 70.6666666666667, 70.6666666666667, 86, 86, 66.6666666666667, 66.6666666666667, 69.3333333333333, 69.3333333333333, 86.6666666666667, 86.6666666666667, 84.6666666666667, 84.6666666666667, 89.3333333333333, 89.3333333333333, 78, 78, 90, 90, 90, 90, 88, 88, 92, 92, 83.3333333333333, 83.3333333333333, 89.3333333333333, 89.3333333333333, 87.3333333333333, 87.3333333333333, 82, 82, 102, 102, 101.333333333333, 101.333333333333, 90.6666666666667, 90.6666666666667, 90, 90, 90, 90, 101.333333333333, 101.333333333333, 100.666666666667, 100.666666666667, 82, 82, 102.666666666667, 102.666666666667, 92.6666666666667, 92.6666666666667, 80, 80, 81.6666666666667, 81.6666666666667, 85, 85, 83.3333333333333, 83.3333333333333, 78.3333333333333, 78.3333333333333, 80, 80, 80, 80, 68.3333333333333, 68.3333333333333, 78.3333333333333, 78.3333333333333, 100, 100, 90, 90, 90, 90, 83.3333333333333, 83.3333333333333, 67.3333333333333, 67.3333333333333, 80, 80, 60.6666666666667, 60.6666666666667, 93.3333333333333, 93.3333333333333, 78.3333333333333, 78.3333333333333, 79.3333333333333, 79.3333333333333, 100.666666666667, 100.666666666667, 103.666666666667, 103.666666666667, 92.6666666666667, 92.6666666666667, 95.6666666666667, 95.6666666666667, 95.6666666666667, 95.6666666666667, 86.6666666666667, 86.6666666666667, 90.6666666666667, 90.6666666666667, 94, 94, 94.6666666666667, 94.6666666666667, 95.3333333333333, 95.3333333333333, 84.6666666666667, 84.6666666666667, 98.6666666666667, 98.6666666666667, 81.3333333333333, 81.3333333333333, 89.3333333333333, 89.3333333333333, 90.6666666666667, 90.6666666666667, 90.6666666666667, 90.6666666666667, 97.3333333333333, 97.3333333333333, 89.3333333333333, 89.3333333333333, 94.6666666666667, 94.6666666666667, 90.6666666666667, 90.6666666666667, 82, 82, 80.6666666666667, 80.6666666666667, 88.6666666666667, 88.6666666666667, 84.6666666666667, 84.6666666666667, 89.3333333333333, 89.3333333333333, 86, 86, 80.6666666666667, 80.6666666666667, 81.3333333333333, 81.3333333333333, 88, 88), oed = c(45, 45, 48, 48, 50, 52, 47, 47, 45, NA, 60, 59, 54, 54, 51, 52, 53, 54, 49, NA, 56, 56, 52, 51, 60.5, 60, 53, 53, 56, NA, 62, 61, 58, 58, 50, 48, 47, 48, 46, NA, 44.5, 42, 56, 58, 48, 48, 38, 39, 50, NA, 66, 66, 62, 62, 58, 58, 68, 66, 62, NA, 66, 66, 68, 68, 66, 66, 68, 66, 66, NA, 63, 62, 49.5, 49, 59, 57, 59, 60, 58, NA, 65, 66, 48, 49, 63, 63, 81, 80, 83, NA, 68, 68, 78, 77, 42, NA, 50, 49, 47, NA, 54, NA, 61, 61, 54, 51, 55, 55, 52, NA, 49, 49, 46, 44, 45, 43, 51, 52, 48, NA, 51, 51, 51.5, 51, 48, 50, 54, 56, 44, NA, 45, 45, 53, 53, 54, 55, 51, 50, 47, NA, 52, 52, 53, 53, 60, 60, 55, 55, 58, NA, 55, 54, 53, 53, 52, 52, 54, 53, 58, NA, 52, 52, 56, 56, 59, 59, 58, 58, 55, NA, 48, 48, 52, 50, 59, 57, 45, 45, 56, NA, 60, 59, 57, 57, 54, 54, 72, 70), dbp = c(55, 60, 60, 65, 70, 80, 60, 60, 60, NA, 94.6666666666667, 96.6666666666667, 86.6666666666667, 86.6666666666667, 81.3333333333333, 82, 94.6666666666667, 94.6666666666667, 87.3333333333333, NA, 92, 91.3333333333333, 89.3333333333333, 86.6666666666667, 93.3333333333333, 91.3333333333333, 94, 94.6666666666667, 88.6666666666667, NA, 97.3333333333333, 94, 94, 91.3333333333333, 65.3333333333333, 64.6666666666667, 69.3333333333333, 80.6666666666667, 71.3333333333333, NA, 90, 70, 60, 68.6666666666667, 80.6666666666667, 70, 67.3333333333333, 70, 84.6666666666667, NA, 80.6666666666667, 82, 89.3333333333333, 88, 74.6666666666667, 77.3333333333333, 90, 56.6666666666667, 87.3333333333333, NA, 90, 86.6666666666667, 90.6666666666667, 90, 90.6666666666667, 90.6666666666667, 89.3333333333333, 90.6666666666667, 86, NA, 84.6666666666667, 80.6666666666667, 81.3333333333333, 87.3333333333333, 87.3333333333333, 82.6666666666667, 100.666666666667, 94.3333333333333, 91.6666666666667, NA, 91.6666666666667, 95, 101.333333333333, 102, 90.6666666666667, 90.6666666666667, 84.6666666666667, 85, 92, NA, 90.6666666666667, 90, 82, 82.6666666666667, 85, NA, 70, 81.6666666666667, 66.6666666666667, NA, 76.6666666666667, NA, 83.3333333333333, 71.6666666666667, 80, 71.6666666666667, 80, 80, 76.6666666666667, NA, 81.6666666666667, 68.3333333333333, 85, 91.6666666666667, 90, 77.6666666666667, 83.3333333333333, 71.6666666666667, 62, NA, 73.3333333333333, 74, 77.3333333333333, 87, 86.6666666666667, 83.3333333333333, 73.3333333333333, 80, 92, NA, 93.6666666666667, 90.6666666666667, 98.6666666666667, 96, 93.3333333333333, 88.6666666666667, 99.3333333333333, 97.3333333333333, 96.6666666666667, NA, 85.3333333333333, 85.3333333333333, 85.3333333333333, 79.3333333333333, 85.3333333333333, 82.6666666666667, 88, 84.6666666666667, 90, NA, 82.3333333333333, 81.3333333333333, 94.3333333333333, 96, 80.6666666666667, 79.3333333333333, 83.3333333333333, 90, 90.6666666666667, NA, 88, 87.3333333333333, 96, 96.6666666666667, 90.6666666666667, 90.6666666666667, 94.6666666666667, 91.3333333333333, 92.6666666666667, NA, 79.3333333333333, 80, 80, 80, 86.6666666666667, 84.6666666666667, 84, 82.6666666666667, 85.3333333333333, NA, 80.6666666666667, 80.6666666666667, 76.6666666666667, 77.3333333333333, 80, 80.6666666666667, 85.3333333333333, 84.6666666666667) ) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ Zo * ,[ @Documents.ControllerDesc Ws,! [h$