calculating baseline hazard and baseline cumulative hazard using the result from a coxtv or coxtp object

The baseline estimation is the baseline hazard at each observed failure time when holding all the covariates to be zero.

baseline(fit)

Arguments

fit

model from coxtv or coxtp.

Value

A list with three components:

time

the unique observed failure times.

hazard

the baseline hazard corresponding to each unique failure time point.

cumulHaz

the cumulative baseline hazard corresponding to each unique failure time point.

Examples

data(ExampleData)
z <- ExampleData$z
time  <- ExampleData$time
event <- ExampleData$event

fit   <- coxtv(event = event, z = z, time = time)
#> Iter 1: Obj fun = -3.2982771; Stopping crit = 1.0000000e+00;
#> Iter 2: Obj fun = -3.2916285; Stopping crit = 2.1424865e-02;
#> Iter 3: Obj fun = -3.2916034; Stopping crit = 8.0884492e-05;
#> Iter 4: Obj fun = -3.2916034; Stopping crit = 1.8232153e-09;
#> Algorithm converged after 4 iterations!
base.est <- baseline(fit)