Nothing
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)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.