mets (version 1.3.7)

phreg: Fast Cox PH regression

Description

Fast Cox PH regression Robust variance is default variance with the summary.

Usage

phreg(formula, data, offset = NULL, weights = NULL, ...)

Arguments

formula

formula with 'Surv' outcome (see coxph)

data

data frame

offset

offsets for Cox model

weights

weights for Cox score equations

...

Additional arguments to lower level funtions

Author

Klaus K. Holst, Thomas Scheike

Details

influence functions (iid) will follow numerical order of given cluster variable so ordering after $id will give iid in order of data-set.

Examples

library(mets)
data(TRACE)
dcut(TRACE) <- ~.
out1 <- phreg(Surv(time,status==9)~vf+chf+strata(wmicat.4)+cluster(id),data=TRACE)
summary(out1)

par(mfrow=c(1,2))
plot(out1)

## computing robust variance for baseline
rob1 <- robust.phreg(out1)
plot(rob1,se=TRUE,robust=TRUE)

## iid decomposition, with scaled influence functions
## for regression parameters
head(iid(out1))
## making iid decomposition of baseline at a specific time-point
Aiiid <- iid(out1,time=30)
head(Aiiid)
## both iid decompositions
dd <- iidBaseline(out1,time=30)
head(dd$beta.iid)
head(dd$base.iid)