hyper_matrix2: Function to build a matrix, used in creating the A and B...

In Giappo/mbd: Multiple Birth Death Diversification

Description

Function to build a matrix, used in creating the A and B operators. It produces the structure q ^ (m - n) * (1 - q) ^ (k + 2 * n-m) * sum_j 2 ^ j choose(k, j) * choose(n, m - n - j)

Usage

hyper_matrix2(n_species, k, q)

Arguments

n_species	number of species.
k	the number of visible species in the phylogeny at a given time.
q	the single-lineage speciation probability at a triggered event.

Author(s)

Hanno Hildenbrandt, adapted by Giovanni Laudanno

Giappo/mbd documentation built on March 3, 2020, 3:36 a.m.