calcP_n: Calculate Z (part of equation 6.4 of Harte 2011)

In BioGeoBEARS: BioGeography with Bayesian (and Likelihood) Evolutionary Analysis in R Scripts

Description

This function is a used by get_probvals, which calculates the Maximum Entropy (Harte (2011)) discrete probability distribution of a number of ordered states (e.g., faces of a 6-sided die) given the mean of many rolls. Here, this is merely used so that a single parameter can control the probability distribution of small versus large descendant areas during cladogenesis.

Usage

  calcP_n(n, lambda1, Z)

Arguments

n	Value of the state (e.g., which of a number of faces on a die, or number of different size classes of geographic range).
lambda1	Lambda parameter (Harte2011).
Z	numeric values from calcZ_part.

Details

See also: Maximum Entropy probability distribution for discrete variable with given mean (and discrete uniform flat prior) http://en.wikipedia.org/wiki/Maximum_entropy_probability_distribution

Value

Prob_n, numeric value of the probability of state n.

Note

Go BEARS!

Author(s)

Nicholas J. Matzke [e-mail]

References

http://phylo.wikidot.com/matzke-2013-international-biogeography-society-poster http://en.wikipedia.org/wiki/Maximum_entropy_probability_distribution

Harte2011

Matzke_2012_IBS

See Also

calcZ_part, maxent, symbolic_to_relprob_matrix_sp

Examples

testval=1
n = 6
lambda1 = 0.5
Z = 1
calcP_n(n, lambda1, Z)

BioGeoBEARS documentation built on May 29, 2017, 8:36 p.m.