CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDescoDocuments.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc TextViews.StdViewDescTextViews.ViewDesc_TextModels.StdModelDescTextModels.ModelDesciaTextModels.AttributesDesc'* *n5#------------------------------------------------------------------------------------------------------------------------------ Version 1: model m1 #------------------------------------------------------------------------------------------------------------------------------ model { for (t in 1:YEARS){ # stochastic part of the likelihood y[t] ~ dbin( p[t], N[t] ) p[t]~dbeta( 1, 1) log.like[t] <- loggam( N[t]+1 ) - loggam( y[t]+1 )- loggam( N[t] - y[t]+1 ) + y[t]*log( p[t] ) + (N[t]-y[t])*log( 1-p[t] ) } # computation of the inverse likelihood inv.like <- exp( - sum(log.like[1:YEARS]) ) for(t in 1:YEARS){ ps[t] <- y[t]/N[t] qs[t] <- 1-ps[t] var[t] <-ps[t]*qs[t]/N[t] log.g[t] <- -0.5*log( 2*3.14 ) - 0.5*log( var[t] ) - 0.5*(p[t]-ps[t])*(p[t]-ps[t])/var[t] # the prior here is one log.prior[t]<-0.0 wl[t] <- log.g[t] - log.like[t] - log.prior[t] } # computation of the inverse wlikelihood wlike <- exp( sum( wl[1:YEARS] ) ) } DATA list(YEARS=8) y[] N[] 554 1183 701 1510 749 1597 868 1924 516 1178 573 1324 978 2173 399 845 END INITS list( p=c(0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*