CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesc-Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc{ sTextViews.StdViewDescTextViews.ViewDescTextModels.StdModelDescTextModels.ModelDesc'TextModels.AttributesDesc'*: *3 model{ for (k in 1:K1){ logor[k] <- log(or[k]) selogor[k] <- log(U[k]/L[k])/(2*1.96) precision.logor[k] <- 1/pow( selogor[k], 2) logor[k] ~ dnorm( theta[k], precision.logor[k] ) theta[k]~dnorm( mu.theta, tau.theta ) OR[k] <- exp(theta[k]) } for (k in 1:K2){ for (i in 1:2){ N[i,k] <- Y[i,1,k]+Y[i,2,k] Y[i,1,k] ~ dbin( p[i,k], N[i,k] ) logit(p[i,k]) <- a[k] + theta[K1+k] * equals(i,1) } theta[K1+k] ~ dnorm( mu.theta, tau.theta ) OR[K1+k] <- exp(theta[K1+k]) } # priors for( k in 1:3) { a[k] ~ dnorm( 0.0, 0.001) } mu.theta ~ dnorm( 0.0, 0.001) tau.theta ~ dgamma( 0.001, 0.001) s2.theta <- 1/tau.theta s.theta <- sqrt( s2.theta ) OR.total <- exp(mu.theta) } INITS list( mu.theta=3.5, tau.theta=1, a=c(0,0,0) ) DATA list( K1=7, K2=3, or=c(3.89, 3.97, 3.88, 17.47, 5.35, 9.1, 3.41), L=c(0.92, 2.2, 2.47, 14.24, 2.44, 5.57, 2.94), U=c(16.3, 7.16, 6.08, 21.43, 11.74, 14.86, 3.96), Y=structure(.Data=c(49, 12, 29, 29958, 0, 4, 33, 89, 171, 70186, 118, 81), .Dim = c(2, 2, 3)) ) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ Zo * ,[ @Documents.ControllerDesc Ws,! [h$