hba.MC {MethComp}R Documentation

A MethComp object from the hba1c data

Description

This object is included for illustrative purposes. It is a result of a 5-hour run using MethComp, with n.iter=100000.

Usage

data(hba.MC)

Format

The format is a MethComp object.

Details

The data are the venous measurements from the hba1c dataset, using the day of analysis as replicate. Measurements are taken to be linked within replicate (=day of analysis).

Examples

data(hba.MC)
attr(hba.MC,"mcmc.par")
print.MethComp(hba.MC)
# One of the chains is really fishy (it's the first one)
trace.MethComp(hba.MC)
trace.MethComp(hba.MC,"beta")
# Try to have a look, excluding the first chain
hba.MCsub <- subset.MethComp(hba.MC,chains=-1)
trace.MethComp(hba.MCsub)
trace.MethComp(hba.MCsub,"beta")
# A MethComp object also has class mcmc.list, so we can use the
# coda functions for covergence diagnostics:
acfplot( subset.MethComp(hba.MC, subset="sigma"))

[Package MethComp version 0.3.0 Index]