CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesc Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDescP H TextViews.StdViewDescTextViews.ViewDesc TextModels.StdModelDescTextModels.ModelDescTextModels.AttributesDesc1$Courier New $Courier New3* #------------------------------------------------------------------------------------------------------------------------------ Version 1 (simple) #------------------------------------------------------------------------------------------------------------------------------ model{ # model's likelihood for (i in 1:n){ y[i] ~ dnorm( mu[i], tau ) mu[i] <- beta0[drug[i]] + beta1[drug[i]]*dose[i]/ranked(dose[],1) } # # prior distributions #normal priors beta0[1] ~ dnorm( 0.0, 0.001) # constant for standard treatment beta0[2] ~ dnorm( 0.0, 0.001) # constant for standard treatment beta1[1] ~ dnorm( 0.0, 0.001) # constant for test treatment beta1[2] ~ dnorm( 0.0, 0.001) # slope ##improper flat priors #beta0[1] ~ dflat() # constant for standard treatment #beta0[2] ~ dflat() # constant for standard treatment #beta1[1] ~ dflat() # constant for test treatment #beta1[2] ~ dflat() # slope tau ~ dgamma( 0.001, 0.001) # precision of regression model s <- 1/sqrt(tau) # standard error of regression # intercept.difference<-beta0[2]-beta0[1] p2 <- step(intercept.difference) # # r2 measures sy2<-sd(y[])*sd(y[]) R2B <- 1 - 1/(tau*sy2) for (i in 1:n){ res[i]<-y[i]-mu[i] } ss2 <- inprod( res[],res[] )/(n-4) R2B2 <- 1 - ss2/sy2 } INITS list( beta1=c(0,0), beta0=c(0,0), tau=1 ) DATA (LIST) list( n=24, y=c(68.8, 67.6, 68.1, 67.6, 69.0, 67.9, 68.6, 68.3, 61.4, 59.8, 62.3, 60.6, 60.9, 60.3, 61.6, 61.8, 53.5, 51.9, 53.6, 52.2, 53.8, 54.9, 54.1, 54.2), dose=c(0.025, 0.025, 0.025, 0.025, 0.025, 0.025, 0.025, 0.025, 0.050, 0.050, 0.050, 0.050, 0.050, 0.050, 0.050, 0.050, 0.100, 0.100, 0.100, 0.100, 0.100, 0.100, 0.100, 0.100), drug=c(1,1,1,1,2,2,2,2, 1,1,1,1,2,2,2,2, 1,1,1,1,2,2,2,2) ) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*