phreg: Parametric Proportional Hazards Regression

Description Usage Arguments Value Note Author(s) Examples

Description

Proportional hazards model with parametric baseline hazard(s). Allows for stratification with different scale and shape in each stratum, and left truncated and right censored data.

Usage

1
2
3
4
phreg(formula = formula(data), data = parent.frame(),
  na.action = getOption("na.action"), dist = "weibull", cuts = NULL, init,
  shape = 0, control = list(eps = 1e-08, maxiter = 20, trace = FALSE),
  singular.ok = TRUE, model = FALSE, x = FALSE, y = TRUE)

Arguments

formula

a formula object, with the response on the left of a ~ operator, and the terms on the right. The response must be a survival object as returned by the Surv function.

data

a data.frame in which to interpret the variables named in the formula.

na.action

a missing-data filter function, applied to the model.frame, after any subset argument has been used. Default is options()$na.action.

dist

Which distribution? Default is "weibull", with the alternatives "ev" (Extreme value), and "pch" (piecewise constant hazards function). A special case like the exponential can be obtained by choosing "weibull" in combination with shape = 1, or "pch" without cuts.

cuts

Only used with dist = "pch". Specifies the points in time where the hazard function jumps. If omitted, an exponential model is fitted.

init

vector of initial values of the iteration. Default initial value is zero for all variables.

shape

If positive, the shape parameter is fixed at that value (in each stratum). If zero or negative, the shape parameter is estimated. If more than one stratum is present in data, each stratum gets its own estimate.

control

a list with components eps (convergence criterion), maxiter (maximum number of iterations), and silent (logical, controlling amount of output). You can change any component without mention the other(s).

singular.ok

Not used.

model

Not used.

x

Return the design matrix in the model object? Default is FALSE

y

Return the response in the model object? Default is TRUE.

Value

A list of class "phreg"

Note

The gompertz distribution is not included here.

Author(s)

Göran Broström

Examples

1
2
fit <- phreg(Surv(enter, exit, event) ~ ses.50, data = oldmort)
fit

goranbrostrom/eha2 documentation built on May 17, 2019, 7:59 a.m.