R/fitYP4.R

In ELYP: Empirical Likelihood Analysis for the Cox Model and Yang-Prentice (2005) Model

fitYP4 <- function(Y, d, Z, beta1=1, beta2= -1, maxiter=60){
#### Do two things: (1) for the data, and given beta1, beta2; 
#### compute the alpha and baseline that max the EL. i.e. (NPMLE|beta1, beta2)
#### (2) Given the baseline and the alpha + beta(s), compute the (log) EL value.

temp1 <- YP4(y=Y, d=d, Z=Z, b1=beta1, b2=beta2, k=maxiter)
 
ELval <- ELcomp(Haz=temp1$Hazw, Sur=temp1$Survival, gam=temp1$gam)

list(EmpLik=ELval, BaselineH=temp1$Hazw, alpha=temp1$alpha)
}