fixcoxph: Remove observations from a proportional hazards regression,...

Description Usage Arguments Value References Examples

View source: R/fixcoxph.R

Description

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

Usage

1
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 regression. 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

\insertRef

kz19PHInfiniteEstimates

Examples

1
2
3
4
5
6
7
8
9
data(breast) # From library coxphf
bcfit<-coxph(Surv(TIME,CENS)~ T+ N+ G+ CD,data=breast,x=TRUE)

fixcoxph(bcfit,bcfit$x,"T",Surv(TIME,CENS)~ T+ N+ G+ CD)

testdat2 <- data.frame(Time=c(4,3,1,1,2,2,3),
  Cen=c(1,1,1,0,0,0,0), Primary=c(0,2,1,1,1,0,0), Sex=c(0,0,0,0,1,1,1))
(bcfit<-coxph(Surv(Time,Cen)~Primary + Sex, testdat2, x=TRUE, ties="breslow"))
fixcoxph(bcfit,bcfit$x,"Primary")

PHInfiniteEstimates documentation built on July 24, 2021, 1:07 a.m.