R/summary.mcmc.list.dc.R

In datacloning/dclone: Data Cloning and MCMC Tools for Maximum Likelihood Methods

summary.mcmc.list.dc <-
function(object, quantiles = c(0.025, 0.25, 0.5, 0.75, 0.975), ...)
{
    class(object) <- "mcmc.list"
    out <- summary(object, quantiles = quantiles, ...)
    out$nclones <- nclones(object)
    if (nchain(object) > 1) {
        abin <- getOption("dcoptions")$autoburnin
        rhat <- gelman.diag(object, autoburnin = abin)$psrf[,1]
        if (inherits(rhat, "try-error"))
            rhat <- NA
    } else rhat <- NA
    if (nvar(object) > 1) {
        out$statistics <- cbind(out$statistics[,1:2],
            "DC SD"=unname(dcsd(object)),
            out$statistics[,3:4],
            "R hat" = unname(rhat))
    } else {
        tmp <- c(out$statistics[1:2],
            "DC SD"=unname(dcsd(object)),
            out$statistics[3:4],
            "R hat" = unname(rhat))
        out$statistics <- matrix(tmp, nrow=1)
        dimnames(out$statistics) <- list(varnames(object), names(tmp))
    }
    class(out) <- "summary.mcmc.dc"
    out
}