summary.MRH: Summarizes the results of an MRH object

In MRH: Multi-Resolution Estimation of the Hazard Rate

Description

This function summarizes the MRH object (either the fitted object or the MCMC chains) and returns estimates of the covariate effects, the hazard rate, the survival curve, the cumulative hazard, and the log ratio for non-propoortional hazards models. If desired, the user can specify the alpha-level for the credible intervals.

Usage

## S3 method for class 'MRH'
summary(object, alpha.level, maxStudyTime, ...)

Arguments

object	object is an MRHobject: Either 1) an MRH fitted object or 2) MCMC chains produced by the estimateMRH routine, converted to an MRH object using MRH() or as.MRH().
alpha.level	The width of the credible intervals for the parameters, with the lower bound calculated as the alpha/2 percentile, and the upper bound calculated as the 1-alpha/2 percentile.
maxStudyTime	The maximum study period (or the censoring time) used in the analysis. This is only required if the MCMC chains are used instead of the fitted MRH model.
...	Arguments to be passed to methods, such as graphical parameters (see par).

Value

d	The hazard increments within each bin (i.e. the cumulative hazard within each bin)
hazardRate	The hazard rate across the bins
SurvivalCurve	The survival curve (i.e. S(t)) over the course of the study.
CumulativeHazard	The cumulative hazard (i.e. H(t)) over the course of the study.
beta	The estimated covariate effects. If the model contains both PH and NPH parameters, the covariate effects will be denoted with 'betaPH' and 'betaNPH'.
H	The parameter H, which is the cumulative hazard at the end of the study.
Rmp	The split parameters.
maxStudyTime	The maximum study time (i.e. the censoring time).

Author(s)

Yolanda Hagar <[e-mail]> and Vanja Dukic

Examples

# These MRH fit portion of the examples are from the 
# estimateMRH() help page.  
# They do not need to be re-run if the objects 
# are already in the active workspace.
data(cancer)
cancer$censorvar = cancer$status - 1
## Not run: 
fit.lung = estimateMRH(formula = Surv(time, censorvar) ~ 
	age + as.factor(sex) + ph.karno, data = cancer,
	M = 3, maxStudyTime = 960, burnIn = 200, maxIter = 1000, 
	thin = 1, outfolder = 'MRH_lung')
## End(Not run)

data(tongue)
## Not run: 
fit.tongue = estimateMRH(formula = Surv(time, delta) ~ 
	nph(type), data = tongue, M = 4, 
	burnIn = 200, maxIter = 2000, thin = 1, outfolder = 'MRH_tongue_nph')
## End(Not run)
	
# Summarize the models
## Not run: 
summary(fit.lung)
summary(fit.lung, alpha.level = .01)

# Read in from the saved output file, converting to an MRH object,
# then summarize.
mcmc.lung = as.MRH(read.table('MRH_lung/MCMCchains.txt', header = TRUE))
summary(mcmc.lung, maxStudyTime = 960)
summary(mcmc.lung, maxStudyTime = 960, alpha.level = .01)

## End(Not run)

MRH documentation built on May 2, 2019, 11:10 a.m.