Description Usage Arguments Value References See Also Examples
supELtest
provides a maximal deviation type statistics that
is better adapted at detecting local differences:
\sup_{t\in U}\{-2\log R(t)\},
where R(t) is an empirical likelihood (EL) ratio that compares two survival functions at each time point t in the set of observed uncensored lifetimes, U.
1 2 3 |
data |
a data frame/matrix with 3 columns. The first column is
the survival time. The second is the censoring indicator. The last is
the grouping variable. An example as the input to |
g1 |
the group with longer survival in one-sided testing with the default value of 1. |
t1 |
pre-specified t_1 based on domain knowledge with the default value of 0 |
t2 |
pre-specified t_2 based on domain knowledge with the default value of ∞ |
sided |
2 if two-sided test, and 1 if one-sided test. It assumes the default value of 2. |
nboot |
number of bootstrap replications in calculating critical values with the defualt value of 1000. |
alpha |
pre-specified significance level of the test with the default value of 0.05 |
compo |
FALSE if taking the standardized square of the difference as the local statisic for two-sided testing, and TRUE if constructing for one-sided testing, but only the positive part of the difference included. It assumes the default value of FALSE. |
seed |
the parameter with the default value of 1011 to |
nlimit |
the splitting unit with the default value of 200. To deal with large data problems, the bootstrap algorithm is
to split the number of bootstrap replicates into |
supELtest
returns a list with three elements:
teststat
the resulting integrated test statistic
critval
the critical value
pvalue
the p-value based on the integrated statistic
H.-w. Chang and I. W. McKeague, "Empirical likelihood based tests for stochastic ordering under right censorship," Electronic Journal of Statistics, Vol. 10, No. 2, pp. 2511-2536 (2016).
hazardcross
, intELtest
, ptwiseELtest
1 2 3 4 5 6 7 8 9 10 11 12 | library(EL2Surv)
supELtest(hazardcross)
## OUTPUT:
## $teststat
## [1] 8.945539
##
## $critval
## [1] 8.738189
##
## $pvalue
## [1] 0.045
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.