| hba.MC {MethComp} | R Documentation |
This object is included for illustrative purposes. It is a result of a
5-hour run using MethComp, with n.iter=100000.
data(hba.MC)
The format is a MethComp object.
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).
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"))