CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesc Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDescY Q TextViews.StdViewDescTextViews.ViewDesc TextModels.StdModelDescTextModels.ModelDesc TextModels.AttributesDesc1$Courier New $Courier Neww #------------------------------------------------------------------------------------------------------------------------------ model 4 #------------------------------------------------------------------------------------------------------------------------------ model{ C<-90.0 # model's likelihood for (i in 1:n){ time[i] ~ dnorm( mu[i], tau ) # stochastic componenent # link and linear predictor mu[i] <- beta[1] + beta[2] * cases[i] + beta[3] * distance[i] log.like[i] <- -0.5*log(2*3.14)-0.5*log(s2)-0.5*(time[i]-mu[i])*(time[i]-mu[i])/s2 } # prior distributions a<-0.001 b<-0.001 tau ~ dgamma( a,b) prior.prec <- 1.0E-3*tau 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) logs2<-log(s2) # inv.like <- exp( -sum(log.like[1:n]) -C ) for (k in 1:3){ log.g[k] <- -0.5*log(2*3.14)-log(g.beta.sd[k])-0.5*pow( (beta[k]-g.beta.mean[k])/g.beta.sd[k], 2) log.prior[k]<- -0.5*log(2*3.14)+0.5*log(prior.prec)-0.5*pow( prior.prec*beta[k], 2) } b0<-g.tau.mean/(g.tau.sd*g.tau.sd) a0<-g.tau.mean*b0 log.g[4] <- a0*log(b0)-loggam(a0)+(a0-1)*log(tau)-b0*tau log.prior[4] <- a*log(b)-loggam(a)+(a-1)*log(tau)-b*tau wlike <- exp( sum(log.g[1:4]) -sum(log.like[1:n]) - sum(log.prior[1:4]) - C) } INITS list( tau=1, beta=c(0,0,0)) DATA (LIST) list( n=25, 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) , g.beta.mean=c(2.38600,1.61300,0.01436), g.tau.mean=0.09371, g.beta.sd=c(1.150000,0.178300,0.003735), g.tau.sd=0.028280 ) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*