Fast additive hazards model with robust standard errors

Fast Lin-Ying additive hazards model with a possibly stratified baseline. Robust variance is default variance with the summary.

Usage

aalenMets(formula, data = data, no.baseline = FALSE, ...)

Arguments

formula

formula with 'Surv' outcome (see coxph)

data

data frame

no.baseline

to fit model without baseline hazard

...

Additional arguments to phreg

Details

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

Author

Thomas Scheike

Examples

data(bmt); bmt$time <- bmt$time+runif(408)*0.001
out <- aalenMets(Surv(time,cause==1)~tcell+platelet+age,data=bmt)
summary(out)
#> 
#>    n events
#>  408    161
#> 
#>  408 clusters
#> coeffients:
#>            Estimate       S.E.    dU^-1/2 P-value
#> tcell    -0.0129601  0.0041293  0.2303828  0.0017
#> platelet -0.0087422  0.0028056  0.1664323  0.0018
#> age       0.0066203  0.0013880  0.0789265  0.0000
#> 
#> exp(coeffients):
#>          Estimate    2.5%  97.5%
#> tcell     0.98712 0.97917 0.9951
#> platelet  0.99130 0.98586 0.9968
#> age       1.00664 1.00391 1.0094
#> 

## out2 <- timereg::aalen(Surv(time,cause==1)~const(tcell)+const(platelet)+const(age),data=bmt)
## summary(out2)