SO.LRT: Likelihood-ratio test statistic

In CorrBin: Nonparametrics with Clustered Binary and Multinomial Data

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)

Author(s)

Aniko Szabo

See Also

SO.trend.test, soControl

Examples

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

CorrBin documentation built on Sept. 29, 2023, 9:07 a.m.