Skip to content

Instantly share code, notes, and snippets.

@wfolta

wfolta/Stan1checking.R

Created October 14, 2013 22:05
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save wfolta/6983012 to your computer and use it in GitHub Desktop.
Save wfolta/6983012 to your computer and use it in GitHub Desktop.
Download ZIP
Stan1 model checking
Raw
Stan1checking.R
library (coda)
gelman.diag (mcmc.list (mcmc (extract (stan.model, pars="alpha", permuted=F)[,1,1]),
mcmc (extract (stan.model, pars="alpha", permuted=F)[,2,1]),
mcmc (extract (stan.model, pars="alpha", permuted=F)[,3,1]),
mcmc (extract (stan.model, pars="alpha", permuted=F)[,4,1])))
heidel.diag (mcmc.list (mcmc (extract (stan.model, pars="alpha", permuted=F)[,1,1]),
mcmc (extract (stan.model, pars="alpha", permuted=F)[,2,1]),
mcmc (extract (stan.model, pars="alpha", permuted=F)[,3,1]),
mcmc (extract (stan.model, pars="alpha", permuted=F)[,4,1])))
geweke.diag (mcmc.list (mcmc (extract (stan.model, pars="alpha", permuted=F)[,1,1]),
mcmc (extract (stan.model, pars="alpha", permuted=F)[,2,1]),
mcmc (extract (stan.model, pars="alpha", permuted=F)[,3,1]),
mcmc (extract (stan.model, pars="alpha", permuted=F)[,4,1])))
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment