bestbeta {PHInfiniteEstimates}

Newton Raphson Fitter for partial likelihood

Description

This function implements the approximate conditional inferential approach of Kolassa and Zhang (2019) to proportional hazards regression.

Usage

bestbeta(fit, exclude = NULL, start = NULL, touse = NA, usecc = FALSE)

Arguments

fit

Output from a Cox PH regression, with x=TRUE and y=TRUE

exclude

data set with stratum and patient number to exclude.

start

Starting value

touse

columns of the design matrix to use.

usecc

Logical variable indicating whether to use a continuity correction, or nuerical variable representing teh continuity correction.

Value

Fitted survival analysis regression parameter of class coxph

References

Kolassa JE, Zhang J (2019). https://higherlogicdownload.s3.amazonaws.com/AMSTAT/fa4dd52c-8429-41d0-abdf-0011047bfa19/UploadedImages/NCB_Conference/Presentations/2019/kolassa_toxslides.pdf. Accessed: 2019-07-14.

Examples

bfit<-coxph(Surv(TIME,CENS)~T+N+CD,data=breast,x=TRUE)
noccfit<-bestbeta(bfit)
bestbeta(bfit,usecc=TRUE,start=noccfit$start)