R/myffitYP3.R

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

Defines functions myffitYP3

Documented in myffitYP3

myffitYP3 <- function(x, myY, myd, myZ) {
x1 <- x[1]
x2 <- x[2]
x3 <- x[3]
tempT <- fitYP3(Y=myY,d=myd,Z=myZ, beta1=c(x3, x1), beta2= c(x3, x2), lam=0, fun=function(t){as.numeric(t <= 0.175)} )
return(- tempT$LogEmpLik)
}

Try the ELYP package in your browser