lambda_N: Computation of lambda_epsilon

Description Usage Arguments Details Author(s) References

View source: R/mcmc_alpha.R

Description

Computes the parameter lambda_epsilon for n in 1:N

Usage

1
lambda_N(epsilon, beta, N, upper = 10000)

Arguments

epsilon

Minimum size of unsampled splits (see appendix 1)

beta

The beta parameter quantifying the balance of the tree

N

Tip number

upper

Upper bound from the integral (see lambda_n expression, default to 10000)

Details

λ_ε = 2 \int_(ε,∞) exp(-(β+n+1)x)(1-exp(-x))^β dx

Author(s)

Odile Maliet, Fanny Gascuel & Amaury Lambert

References

Maliet O., Gascuel F., Lambert A. (2018) Ranked tree shapes, non-random extinctions and the loss of phylogenetic diversity, bioRxiv 224295, doi: https://doi.org/10.1101/224295


apTreeshape documentation built on Jan. 8, 2021, 2:07 a.m.