fitYP3 {ELYP} R Documentation

## Compute Baseline Hazard for the Given Data, Given Parameters: beta1, beta2, lam, and fun. Also, Given the Baseline, Compute the empirical likelihood value.

### Description

This function fits the baseline for given beta1 beta2 and lam; and then compute the empirical likelihood.

### Usage

fitYP3(Y, d, Z, beta1, beta2, lam, fun)


### Arguments

 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.

### Details

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.

### Value

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

### References

Zhou, M. (2002). Computing censored empirical likelihood ratio by EM algorithm. Tech Report, Univ. of Kentucky, Dept of Statistics

### Examples

## 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)


[Package ELYP version 0.7-5 Index]