AICstats_2models: Calculate all the AIC and LRT stats between two models
In BioGeoBEARS: BioGeography with Bayesian (and Likelihood) Evolutionary Analysis in R Scripts

Description Usage Arguments Details Value Note Author(s) References See Also Examples

View source: R/BioGeoBEARS_generics_v1.R

The Likelihood Ratio Test (LRT) is a standard method for testing whether or not the data likelihood conferred by a more complex is significantly better than the data likelihood conferred by the simpler model. See lrttest and lrttest_on_summary_table for more discussion.

1	AICstats_2models(LnL_1, LnL_2, numparams1, numparams2)

`LnL_1`	Log-likelihood of more complex model.
`LnL_2`	Log-likelihood of simpler complex model.
`numparams1`	Number of free parameters of the more complex model.
`numparams2`	Number of free parameters of the less complex model.

See Burnham et al. (2002) and http://www.brianomeara.info/tutorials/aic for discussion of AIC and its uses.

This function assumes that LnL_1 and numparams1 refer to the more complex model, and that LnL_2 and numparams2 refer to the simpler model nested within the more complex one.

LRT_AIC_results A table of LRT and AIC results.

Go BEARS!

Nicholas J. Matzke matzke@berkeley.edu

http://phylo.wikidot.com/matzke-2013-international-biogeography-society-poster

Burnham_Anderson_2002

Matzke_2012_IBS

lrttest, lrttest_on_summary_table

test=1

Loading required package: rexpokit
Loading required package: SparseM

Attaching package: 'SparseM'

The following object is masked from 'package:base':

    backsolve

Loading required package: Rcpp
Loading required package: cladoRcpp
Loading required package: ape
Loading required package: phylobase

Attaching package: 'phylobase'

The following object is masked from 'package:ape':

    edges