get_gamma_error_fun: Get an error function that uses the difference in gamma...

View source: R/get_gamma_error_fun.R

get_gamma_error_funR Documentation

Get an error function that uses the difference in gamma statistic.

Description

Get an error function that uses the difference in gamma statistic.

Usage

get_gamma_error_fun()

Value

an error function

Note

the gamma statistic can be found in Pybus and Harvey, 2000, 'Testing macro-evolutionary models using incomplete molecular phylogenies.', Proc. R. Soc. Lond. B, 267, 2267-2272.

Author(s)

Richèl J.C. Bilderbeek

Examples

error_fun <- get_gamma_error_fun()

phylogeny <- ape::read.tree(text = "((A:1.5, B:1.5):1.5, C:3.0);")

tree_1 <- ape::read.tree(text = "((A:1.0, B:1.0):2.0, C:3.0);")
tree_2 <- ape::read.tree(text = "((A:2.0, B:2.0):1.0, C:3.0);")
trees <- c(tree_1, tree_2)

error_fun(phylogeny, c(phylogeny))
error_fun(phylogeny, c(tree_1))
error_fun(phylogeny, c(tree_2))

richelbilderbeek/pirouette documentation built on Oct. 18, 2023, 3:09 p.m.