Log-likelihood ratio test for a Dirichlet mean vector {Compositional}R Documentation

Log-likelihood ratio test for a Dirichlet mean vector

Description

Log-likelihood ratio test for a Dirichlet mean vector.

Usage

dirimean.test(x, a)

Arguments

x

A matrix with the compositional data. No zero values are allowed.

a

A compositional mean vector. The concentration parameter is estimated at first. If the elements do not sum to 1, it is assumed that the Dirichlet parameters are supplied.

Details

Log-likelihood ratio test is performed for the hypothesis the given vector of parameters "a" describes the compositional data well.

Value

If there are no zeros in the data, a list including:

param

A matrix with the estimated parameters under the null and the alternative hypothesis.

loglik

The log-likelihood under the alternative and the null hypothesis.

info

The value of the test statistic and its relevant p-value.

Author(s)

Michail Tsagris

R implementation and documentation: Michail Tsagris mtsagris@uoc.gr and Giorgos Athineou <gioathineou@gmail.com>

References

Ng Kai Wang, Guo-Liang Tian and Man-Lai Tang (2011). Dirichlet and related distributions: Theory, methods and applications. John Wiley & Sons.

See Also

sym.test, diri.nr, diri.est, rdiri, ddiri

Examples

x <- rdiri( 100, c(1, 2, 3) )
dirimean.test(x, c(1, 2, 3) )
dirimean.test( x, c(1, 2, 3)/6 )

[Package Compositional version 5.2 Index]