CDOoDocuments.StdDocumentDescDocuments.DocumentDescContainers.ViewDescViews.ViewDescStores.StoreDesc Documents.ModelDescContainers.ModelDescModels.ModelDescStores.ElemDesc TextViews.StdViewDescTextViews.ViewDesc TextModels.StdModelDescTextModels.ModelDesc  TextModels.AttributesDesc+ArialPS '*PS 6* Arial\model{ # PRIORS alpha[1] <- 0; # zero contrast for baseline food for (k in 2 : K) { alpha[k] ~ dnorm(0, 0.00001)} # vague priors # Loop around lakes: for (k in 1 : K){ beta[1, k] <- 0 } # corner-point contrast with first lake for (i in 2 : I) { beta[i, 1] <- 0; # zero contrast for baseline food for (k in 2 : K) { beta[i, k] ~ dnorm(0, 0.00001)} # vague priors } # Loop around sizes: for (k in 1 : K){ gamma[1, k] <- 0} # corner-point contrast with first size for (j in 2 : J) { gamma[j, 1] <- 0 ; # zero contrast for baseline food for ( k in 2 : K){ gamma[j, k] ~ dnorm(0, 0.00001) } # vague priors } # LIKELIHOOD # i = loop around lakes, j = loop around sizes, k = loop around foods for (i in 1 : I) { for (j in 1 : J) { X[i,j,1:K] ~ dmulti( p[i,j,1:K] , n[i,j] ) n[i,j] <- sum(X[i,j,]) for (k in 1:K) { p[i,j,k] <- phi[i,j,k] / sum(phi[i,j,]) log(phi[i,j,k]) <- alpha[k] + beta[i,k] + gamma[j,k] } } } } list( I = 4, J = 2, K = 5, X = structure(.Data = c(23, 4, 2, 2, 8, 7, 0, 1, 3, 5, 5, 11, 1, 0, 3, 13, 8, 6, 1, 0, 5, 11, 2, 1, 5, 8, 7, 6, 3, 5, 16, 19, 1, 2, 3, 17, 1, 0, 1, 3), .Dim = c(4, 2, 5))) list(alpha = c(NA, 0, 0, 0, 0), beta = structure(.Data = c(NA, NA, NA, NA, NA, NA, 0, 0, 0, 0, NA, 0, 0, 0, 0, NA, 0, 0, 0, 0), .Dim = c(4, 5)), gamma = structure(.Data = c(NA, NA, NA, NA, NA, NA, 0, 0, 0, 0), .Dim = c(2, 5)) )TextControllers.StdCtrlDescTextControllers.ControllerDescContainers.ControllerDescControllers.ControllerDesc TextRulers.StdRulerDescTextRulers.RulerDescTextRulers.StdStyleDescTextRulers.StyleDescZTextRulers.AttributesDesc$0ZGo * <[ @Documents.ControllerDesc vz pk ~