Learn R Programming

SurvRegCensCov (version 1.7)

ConvertWeibull: Transformation of survreg output for the Weibull distribution

Description

Transforms output from survreg using the Weibull distribution to a more natural parameterization. See details and the vignette for more information.

Usage

ConvertWeibull(model, conf.level = 0.95)

Value

vars

A matrix containing the values of the transformed parameters and their standard errors

HR

A matrix containing the hazard ratios for the covariates, and \(1 - \)conf.level\( / 2\) confidence intervals.

ETR

A matrix containing the event time ratios for the covariates, and \(1 - \)conf.level\( / 2\) confidence intervals.

Arguments

model

A survreg model, with dist = "weibull".

conf.level

Confidence level used to produce two-sided \(1-\alpha/2\) confidence intervals for the hazard and event time ratios.

Author

Sarah R. Haile, Epidemiology, Biostatistics and Prevention Institute (EBPI), University of Zurich, sarah.haile@uzh.ch

Details

The survreg function fits a Weibull accelerated failure time model of the form

$$\log t = \mu + \alpha^T Z + \sigma W,$$

where \(Z\) is a matrix of covariates, and \(W\) has the extreme value distribution, \(\mu\) is the intercept, \(\alpha\) is a vector of parameters for each of the covariates, and \(\sigma\) is the scale. The usual parameterization of the model, however, is defined by hazard function

$$h(t|Z) = \gamma \lambda t^{\gamma - 1} \exp(\beta^T Z).$$

The transformation is as follows: \(\gamma = 1/\sigma\), \(\lambda = \exp(-\mu/\sigma)\), and \(\beta=-\alpha/\sigma\), and estimates of the standard errors can be found using the delta method.

The Weibull distribution has the advantage of having two separate interpretations. The first, via proportional hazards, leads to a hazard ratio, defined by \(\exp \beta\). The second, of accelerated failure times, leads to an event time ratio (also known as an acceleration factor), defined by \(\exp (-\beta/\gamma)\).

Further details regarding the transformations of the parameters and their standard errors can be found in Klein and Moeschberger (2003, Chapter 12). An explanation of event time ratios for the accelerated failure time interpretation of the model can be found in Carroll (2003). A general overview can be found in the vignette("weibull") of this package.

References

Carroll, K. (2003). On the use and utility of the Weibull model in the analysis of survival data. Controlled Clinical Trials, 24, 682--701.

Klein, J. and Moeschberger, M. (2003). Survival analysis: techniques for censored and truncated data. 2nd edition, Springer.

See Also

This function is used by WeibullReg.

Examples

Run this code
data(larynx)
ConvertWeibull(survreg(Surv(time, death) ~ stage + age, larynx), conf.level = 0.95)

Run the code above in your browser using DataLab