best2nullLRT: Calculate likelihood ratio tests on GMYC objects

In spiderDev: Species Identity and Evolution in R--development

Description

Calculation of likelihood ratio tests on GMYC objects. Includes comparing the best model with the null, and comparing single-threshold and multiple-threshold models.

Usage

best2nullLRT(gmyc)
sing2multLRT(test1, test2)

Arguments

gmyc	A ‘gmyc’ object created with gmyc.
test1	A ‘gmyc’ object created with gmyc.
test2	A ‘gmyc’ object created with gmyc.

Details

best2null compares the best model with the null model.

sing2multLRT compares a multiple-threshold model with a single-threshold model. It is a re-implementation of compare in the splits package, to facilitate the use of data in downstream applications.

Value

A list with the following elements:

chi.sq	2 times the difference between the likelihoods
df	Degrees of freedom
probability	Probability based on a chi-square distribution

Author(s)

Samuel Brown <[e-mail]>

See Also

gmyc

Examples

data(doloGMYC)
data(doloMultGMYC)

best2nullLRT(doloGMYC)

sing2multLRT(doloGMYC, doloMultGMYC)

spiderDev documentation built on May 2, 2019, 5:23 p.m.