Learn R Programming

phytools (version 0.2-40)

threshDIC: Deviance Information Criterion from the threshold model

Description

This function computes the Deviance Information Criterion from the MCMC object returned by ancThresh.

Usage

threshDIC(tree, x, mcmc, burnin=NULL, sequence=NULL, method="pD")

Arguments

tree
phylogenetic tree.
x
a named vector containing discrete character states; or a matrix containing the tip species, in rows, and probabilities of being in each state, in columns.
mcmc
list object returned by ancThresh.
burnin
number of generations (not samples) to exclude as burn in; if not supplied then 20 percent of generations are excluded.
sequence
assumed ordering of the discrete character state. If not supplied and x is a vector then numerical-alphabetical order is assumed; if not supplied and x is a matrix, then the column order of x is used.
method
method for computing the effective number of parameters (options are "pD" and "pV").

Value

  • a vector containing the mean deviance and deviance for the parameter means, the effective number of parameters, and the DIC.

See Also

ancThresh