Hi JAGS community, I am getting familiar with the Bayesian reasoning and I am stuck, after having searched for tips/examples/codes throughout the web, at implementing a Posterior Predictive Check for this Dirichlet-Multinomial model. Any hel/tip/indication on how to accomplish this task would be extremely appreciated. Stef model{ #loop observations for (i in 1:numObs){ choice[i]~dcat(theta[catid[i],1:numSubstrato]) theta[i,1:numSubstrato] ~ ddirch(vcounts[1:numSubstrato]+0.05) } alphasbstr <- rep(1,numSubstrato)...