fitYP3 {ELYP} | R Documentation |
This function fits the baseline for given beta1 beta2 and lam; and then compute the empirical likelihood.
fitYP3(Y, d, Z, beta1, beta2, lam, fun)
Y |
a vector containing the observed survival times. |
d |
a vector containing the censoring indicators, 1-uncensored; 0-right censored. |
Z |
a matrix of covariates (Xi and Zi)... |
beta1 |
a vector = (alpha, beta1). |
beta2 |
a vector = (alpha, beta2). |
lam |
a scalar. tilting parameter for the baseline. When lam=0, then there is no tilting. |
fun |
a function. It determine what feature of the baseline lam tilts. |
This function computes the log empirical likelihood. The parameters are given: beta1, beta2 and lam.
What baseline feature the lam corresponds to is determined by the fun(t), as in int f(t) dH(t). This integral value for the lam is also in the output as Mulam.
A list with the following components:
LogEmpLik |
this is actually the log empirical likelihood value. |
Mulam |
The value of int f(t) d H(t) for corresponding lam. This is also called a baseline feature. |
BaseHazw |
The baseline hazard jumps. |
Mai Zhou
Zhou, M. (2002). Computing censored empirical likelihood ratio by EM algorithm. Tech Report, Univ. of Kentucky, Dept of Statistics
## censored regression with one right censored observation.
## we check the estimation equation, with the MLE inside myfun7.
y <- c(3, 5.3, 6.4, 9.1, 14.1, 15.4, 18.1, 15.3, 14, 5.8, 7.3, 14.4)
x <- c(1, 1.5, 2, 3, 4, 5, 6, 5, 4, 1, 2, 4.5)
d <- c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0)