CorrBin (version 1.3)

SO.LRT: Likelihood-ratio test statistic

Description

SO.LRT computes the likelihood ratio test statistic for stochastic ordering against equality assuming marginal compatibility for both alternatives. Note that this statistic does not have a $\chi^2$ distribution, so the p-value computation is not straightforward. The SO.trend.test function implements a permutation-based evaluation of the p-value for the likelihood-ratio test.

Usage

SO.LRT(cbdata, control = soControl())

Arguments

cbdata

a CBData object

control

an optional list of control settings, usually a call to soControl. See there for the names of the settable control values and their effect.

Value

• The value of the likelihood ratio test statistic is returned with two attributes:
• ll0
  the log-likelihood under $H_0$ (equality)
• ll1
  the log-likelihood under $H_a$ (stochastic order)

See Also

SO.trend.test, soControl

Examples

data(shelltox)
LRT <- SO.LRT(shelltox, control=soControl(max.iter = 100, max.directions = 50))
LRT