Alternative test of the Buckley-James estimator by Empirical Likelihood

Share:

Description

Use the empirical likelihood ratio (alternative form) and Wilks theorem to test if the regression coefficient is equal to beta, based on the estimating equations.

The log empirical likelihood been maximized is

∑_{j=1}^n \log p_j ; ~ ∑ p_j =1

where the probability p_j is for the jth martingale differences of the estimating equations.

Usage

1
bjtestII(y, d, x, beta)

Arguments

y

a vector of length N, containing the censored responses.

d

a vector of length N. Either 1's or 0's. d=1 means y is uncensored; d=0 means y is right censored.

x

a matrix of size N by q.

beta

a vector of length q. The value of the regression coefficient to be tested in the model Y_i = β x_i + ε_i

Details

The above likelihood should be understood as the likelihood of the martingale difference terms. For the definition of the Buckley-James martingale or estimating equation, please see the (2015) book in the reference list.

The estimation equations used when maximize the empirical likelihood is

0 = ∑ d_i Δ F(e_i) (x \cdot m[,i])/(n w_i)

where e_i is the residuals, other details are described in the reference book of 2015 below.

The final test is carried out by el.test. So the output is similar to the output of el.test.

Value

A list with the following components:

"-2LLR"

the -2 loglikelihood ratio; have approximate chisq distribution under H_o.

logel2

the log empirical likelihood, under estimating equation.

logel

the log empirical likelihood of the Kaplan-Meier of e's.

prob

the probabilities that max the empirical likelihood under estimating equation.

Author(s)

Mai Zhou.

References

Zhou, M. (2016) Empirical Likelihood Methods in Survival Analysis. CRC Press.

Buckley, J. and James, I. (1979). Linear regression with censored data. Biometrika, 66 429-36.

Zhou, M. and Li, G. (2008). Empirical likelihood analysis of the Buckley-James estimator. Journal of Multivariate Analysis, 99, 649–664.

Zhu, H. (2007) Smoothed Empirical Likelihood for Quantiles and Some Variations/Extension of Empirical Likelihood for Buckley-James Estimator, Ph.D. dissertation, University of Kentucky.

Examples

1
2
data(myeloma)
bjtestII(y=myeloma[,1], d=myeloma[,2], x=cbind(1, myeloma[,3]), beta=c(37, -3.4))

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.