CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesc Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc#  TextViews.StdViewDescTextViews.ViewDesc TextModels.StdModelDescTextModels.ModelDescTextModels.AttributesDesc1$Courier New $Courier New,M   #------------------------------------------------------------------------------------------------------------------------------ all models depending on g #------------------------------------------------------------------------------------------------------------------------------ model{ C<-90.0 pi<-3.141593 # model's likelihood for (i in 1:n){ time[i] ~ dnorm( mu[i], tau ) # stochastic componenent # link and linear predictor mu[i] <- beta[1] + g[1] * beta[2] * cases[i] + g[2] * beta[3] * distance[i] log.like[i] <- -0.5*log(2*pi)-0.5*log(s2)-0.5*(time[i]-mu[i])*(time[i]-mu[i])/s2 like[i] <- exp( log.like[i] ) } dm <- 2 + g[1] + g[2] Deviance <- -2*sum(log.like[1:n]) AIC <- Deviance + dm*2 BIC <- Deviance + dm*log(n) L <- prod( like[1:n] ) # prior distributions a<-0.001 b<-0.001 tau ~ dgamma(a,b) prior.prec <- tau/c2 beta[1] ~ dnorm( 0.0, prior.prec ) beta[2] ~ dnorm( 0.0, prior.prec ) beta[3] ~ dnorm( 0.0, prior.prec ) # definition of sigma s2<-1/tau s <-sqrt(s2) } INITS list( tau=1, beta=c(0,0,0)) # use g to define the fitted model # m1 g=c(0,0) # m2 g=c(1,0) # m3 g=c(0,1) # m4 g=c(1,1) DATA (LIST) list( n=25, g=c(1,1), c2=1000, time = c(16.68, 11.5, 12.03, 14.88, 13.75, 18.11, 8, 17.83, 79.24, 21.5, 40.33, 21, 13.5, 19.75, 24, 29, 15.35, 19, 9.5, 35.1, 17.9, 52.32, 18.75, 19.83, 10.75), distance = c(560, 220, 340, 80, 150, 330, 110, 210, 1460, 605, 688, 215, 255, 462, 448, 776, 200, 132, 36, 770, 140, 810, 450, 635, 150), cases = c( 7, 3, 3, 4, 6, 7, 2, 7, 30, 5, 16, 10, 4, 6, 9, 10, 6, 7, 3, 17, 10, 26, 9, 8, 4) ) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*