SO.LRT: Likelihood-ratio test statistic

In CorrBin: Nonparametrics with Clustered Binary and Multinomial Data

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-squared 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 H0 (equality)
ll1	the log-likelihood under Ha (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 May 2, 2019, 4:46 p.m.