diag.mcmc: MCMC diagnostics plot.

In lixiangchun/lxctk: Li Xiangchun's tool-kit (lxctk)

Description

MCMC diagnostics plot.

Usage

diag.mcmc( codaObject, parName=varnames(codaObject)[1], 
          DBDAplColors = c("skyblue","black","royalblue","steelblue"),
          figName=NULL)
##diagStanFit(stanFit,parName,saveName=NULL,saveType="jpg") ## Does not work yet

Arguments

codaObject	A coda::mcmc.list object or rstan::stanfit object.
stanFit	A rstan::stanfit object.
parName	The parameter name.
DBDAplColors	Colors to be used for plotting.
figName	A output figure name, must be ended with .png.

See Also

best.robust_t_test,best.student_t

Examples

library(lxctk)
library(coda)
library(rstan)

x <- rt(1000, df=1)
fit <- best.student_t(x, stanDso)
codaObject <- stanfit2mcmc.list(fit)
diag.mcmc(codaObject)


data('x')
##fit <- best.student_t(x, stanDso) ## The above error will be generated
fit <- best.student_t(x, stanDso, unifLo=0, unifHi=.Machine$double.xmax) ## No errors
codaObject <- stanfit2mcmc.list(fit)
diag.mcmc(codaObject)

lixiangchun/lxctk documentation built on May 21, 2019, 6:44 a.m.