pclades: Joint probability of two clades

View source: R/pclades.R

pcladesR Documentation

Joint probability of two clades

Description

Joint probability of two clades

Usage

pclades(a, b, n, case = 1, log = FALSE)

Arguments

a

size of clade A.

b

size of clade B.

n

number of taxa.

case

see details

log

logical if TRUE (default is FALSE) returns log probability

Details

case 1: A == B;

case 2: A is contained in B;

case 3: B is contained in A;

case 4: A and B are disjoint and a + b = n;

case 5: A and B are disjoint and a + b < n;

case 6: A and B are incompatible.

Value

(log) joint probability of clades A and B given their sizes.

Examples

pclades(a = 2,  b = 3, n = 5,  case = 2)

maxbiostat/cladeCorrelation documentation built on July 28, 2023, 2:42 a.m.