CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDescDocuments.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc TextViews.StdViewDescTextViews.ViewDescTextModels.StdModelDescTextModels.ModelDescTextModels.AttributesDesc'* *n#------------------------------------------------------------------------------------------------------------------------------ Version 1: model m0 #------------------------------------------------------------------------------------------------------------------------------ model { for (t in 1:YEARS){ # stochastic part of the likelihood y[t] ~ dbin( p, N[t] ) log.like[t] <- loggam( N[t]+1 ) - loggam( y[t]+1 )- loggam( N[t] - y[t]+1 ) + y[t]*log( p ) + (N[t]-y[t])*log( 1-p ) } p~dbeta( 1, 1) # computation of the inverse likelihood inv.like <- exp( - sum(log.like[1:YEARS]) ) ps <- sum( y[1:YEARS] )/sum( N[1:YEARS] ) qs <- 1-ps var <-ps*qs/sum( N[1:YEARS] ) log.g <- -0.5*log( 2*3.14 ) - 0.5*log( var ) - 0.5*(p-ps)*(p-ps)/var log.prior<-0.0 # computation of the inverse wlikelihood wlike <- exp( log.g-sum( log.like[1:YEARS] )-log.prior ) } 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=0.5) TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$ ZGo * ,[ @Documents.ControllerDesc t]s ' `h*