inst/examples/Jags-Ydich-Xnom1subj-MbernBetaModelCompPseudoPrior.R

In variani/dbda: R code from Doing Bayesian Data Analysis (DBDA)

# Jags-Ydich-Xnom1subj-MbernBetaModelComp.R
# Accompanies the book:
#   Kruschke, J. K. (2014). Doing Bayesian Data Analysis: 
#   A Tutorial with R, JAGS, and Stan. 2nd Edition. Academic Press / Elsevier.


JagsYdichXnom1subjMbernBetaModelCompPseudoPrior <- function()
{
  out <- list()
  
  oldClass(out) <- "JagsYdichXnom1subjMbernBetaModelCompPseudoPrior"
  return(out)
}

mod <- JagsYdichXnom1subjMbernBetaModelCompPseudoPrior()

#------------------------------------------------------------------------------- 
# Generate the MCMC chain:
codaSamples <- genMCMC(mod)
  
#------------------------------------------------------------------------------- 
# Display diagnostics of chain:

parameterNames = varnames(codaSamples) # get all parameter names
for ( parName in parameterNames ) {
  diagMCMC( codaSamples , parName=parName)
}

#------------------------------------------------------------------------------
# EXAMINE THE RESULTS.

plotMCMC(mod, codaSamples)