PHInfiniteEstimates (version 1.6)

fixcoxph: Remove observations from a proportional hazards regression, and return the fit of the reduced model.

Description

This function implements the approximate conditional inferential approach of kz19;textualPHInfiniteEstimates to proportional hazards regression.

Usage

fixcoxph(randdat, xxx, iv, verbose = FALSE)

Arguments

randdat

A list with at least the component y, representing the Surv() object. I expect that this will be output from an initial non-convergent regression.

xxx

a design matrix for the gregression. I expect that this will be the $x component of the output from an initial non-convergent regression, run with x=TRUE .

iv

name of the variable of interest, as a character string

verbose

logical flag governing printing.

Value

Fitted survival analysis regression parameter of class coxph, fitted form data set with observations forcing infinite estimation removed.

References

kz19PHInfiniteEstimates

Examples

# NOT RUN {
data(breast) # From library coxphf
bcfit<-coxph(Surv(TIME,CENS)~ T+ N+ G+ CD,data=breast,x=TRUE)
# }
# NOT RUN {
fixcoxph(bcfit,bcfit$x,"T",Surv(TIME,CENS)~ T+ N+ G+ CD)
# }