CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesc Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc TextViews.StdViewDescTextViews.ViewDesc TextModels.StdModelDescTextModels.ModelDescTextModels.AttributesDesc1$Courier New3 $Courier New #------------------------------------------------------------------------------------------------------------------------------ Version 3 (using CR/STZ design matrix ) [rescaled] #------------------------------------------------------------------------------------------------------------------------------ model{ for (i in 1:n){ # # creating the design matrix X[i,1]<-1.0 # beta1=constant term X[i,2]<-dose[i]/ranked(dose[],1) # beta2=rescaled dose for standard treatment X[i,3]<- equals( drug[i], 2 )*X[i,2] # beta3=CR dummy for interaction (slope difference) #X[i,3]<- (equals( drug[i], 2 )-equals( drug[i], 1 ))*X[i,2] # beta3=STZ dummy for interaction (slope difference) # model likelihood y[i] ~ dnorm( mu[i], tau ) mu[i] <- inprod( beta[], X[i,]) } # rho <- (beta[3]+beta[2])/beta[2] # relative potency in CR #rho <-(beta[2]+beta[3])/(beta[2]-beta[3]) # relative potency in stz # potency estimate potency <- rho * 1.2 # # prior distributions for (j in 1:3){ beta[j]~dnorm( 0.0, 0.0001) } #for (j in 1:3){ beta[j]~dflat() } tau ~ dgamma( 0.001, 0.001) # precision of regression model s <- 1/sqrt(tau) # standard error of regression } INITS list( beta=c(0,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*