Use the empirical likelihood ratio and Wilks theorem to test if the regression coefficient is equal to beta, based on the rank estimator for the AFT model.
The log empirical likelihood been maximized is
∑_{d=1} \log Δ F(e_i) + ∑_{d=0} \log [1F(e_i)];
where e_i are the residuals.
1  RankRegTest(y, d, x, beta, type="Gehan")

y 
a vector of length N, containing the censored responses. 
d 
a vector (length N) of 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 
.
type 
default to Gehan type. The other option is Logrank type. 
The estimator of beta can be obtained by function
rankaft( )
in the package rankreg
. But here you may test other values of
beta. If you test the beta value that is obtained from the rankaft( )
,
then the 2LLR should be 0 and the pvalue should be 1.
The above likelihood should be understood as the likelihood of the error term, so in the regression model the error e_i should be iid.
The estimation equation used when maximize the empirical likelihood is
0 = ∑_i φ (e_i) d_i Δ F(e_i) (x_i  \bar x_i )/(n w_i)
which was described in detail in the references below.
See also the function RankRegTestH
, which is based on the hazard likelihood.
A list with the following components:
"2LLR" 
the 2 loglikelihood ratio; should have approximate chisq distribution under H_o. 
logel2 
the log empirical likelihood, under estimating equation. 
logel 
the log empirical likelihood of the KaplanMeier of e's. 
prob 
the probabilities that max the empirical likelihood under rank estimating equation constraint. 
Mai Zhou.
Kalbfleisch, J. and Prentice, R. (2002) The Statistical Analysis of Failure Time Data. 2nd Ed. Wiley, New York. (Chapter 7)
Jin, Z., Lin, D.Y., Wei, L. J. and Ying, Z. (2003). Rankbased inference for the accelerated failure time model. Biometrika, 90, 34153.
Zhou, M. (2005). Empirical likelihood analysis of the rank estimator for the censored accelerated failure time model. Biometrika, 92, 49298.
1 2 3 4 
Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
Please suggest features or report bugs with the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.