Learn R Programming

CorrBin (version 1.6.2)

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())

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)

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.

Author

Aniko Szabo

See Also

SO.trend.test, soControl

Examples

Run this code

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

Run the code above in your browser using DataLab