Calculates the phylogenetic diversity of order *q* of a probability vector.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | ```
PhyloDiversity(NorP, q = 1, Tree, Normalize = TRUE, Correction = "Best",
CheckArguments = TRUE, Ps = NULL, Ns = NULL)
bcPhyloDiversity(Ns, q = 1, Tree, Normalize = TRUE, Correction = "Best",
CheckArguments = TRUE)
## S3 method for class 'ProbaVector'
PhyloDiversity(NorP, q = 1, Tree, Normalize = TRUE, Correction = "Best",
CheckArguments = TRUE, Ps = NULL, Ns = NULL)
## S3 method for class 'AbdVector'
PhyloDiversity(NorP, q = 1, Tree, Normalize = TRUE, Correction = "Best",
CheckArguments = TRUE, Ps = NULL, Ns = NULL)
## S3 method for class 'integer'
PhyloDiversity(NorP, q = 1, Tree, Normalize = TRUE, Correction = "Best",
CheckArguments = TRUE, Ps = NULL, Ns = NULL)
## S3 method for class 'numeric'
PhyloDiversity(NorP, q = 1, Tree, Normalize = TRUE, Correction = "Best",
CheckArguments = TRUE, Ps = NULL, Ns = NULL)
is.PhyloDiversity(x)
## S3 method for class 'PhyloDiversity'
summary(object, ...)
``` |

`Ps` |
A probability vector, summing to 1. |

`Ns` |
A numeric vector containing species abundances. |

`NorP` |
A numeric vector, an integer vector, an abundance vector ( |

`q` |
A number: the order of diversity. Default is 1. |

`Tree` |
An object of class |

`Normalize` |
If |

`Correction` |
A string containing one of the possible corrections: |

`CheckArguments` |
Logical; if |

`x` |
An object to be tested or plotted |

`object` |
A |

`...` |
Additional arguments to be passed to the generic methods. |

The phylogenetic entropy is its generalization of HCDT entropy to unequal species distances (Pavoine et al., 2009).

Diversity is obtained by transforming generalized entropy.

Bias correction requires the number of individuals to estimate sample `Coverage`

. Use `bcPhyloDiversity`

and choose the `Correction`

.

The functions are designed to be used as simply as possible. `PhyloDiversity`

is a generic method. If its first argument is an abundance vector, an integer vector or a numeric vector which does not sum to 1, the bias corrected function `bcPhyloDiversity`

is called. Explicit calls to `bcPhyloDiversity`

(with bias correction) or to `PhyloDiversity.ProbaVector`

(without correction) are possible to avoid ambiguity. The `.integer`

and `.numeric`

methods accept `Ps`

or `Ns`

arguments instead of `NorP`

for backward compatibility.

An object of class `PhyloDiversity`

is a list:

`Distribution` |
The distribution used to calculate diversity |

`Function` |
The function used to calculate diversity |

`Tree` |
The functional or phylogenetic tree used to calculate diversity |

`Normalized` |
Logical. Indicates whether phylodiversity is normalized or proportional to the height of the tree. |

`Type` |
The type of diversity ("alpha", "beta" or "gamma"). |

`Order` |
The order of diversity |

`Cuts` |
A named vector containing values of neutral diversity along the tree. Names are cut ends, |

`Total` |
A value equal the total diversity (obtained by transforming the total normalized entropy), multiplied by the tree height if |

`is.PhyloDiversity`

returns `TRUE`

if the object is of class `PhyloDiversity`

.

`summary.PhyloDiversity`

returns a summary of the object's value.

`PhyloDiversity`

objects can be plotted by `plot.PhyloValue`

because `PhyloDiversity`

objects are also of class `PhyloValue`

.

The tree must contain all species of the probability vector. If it contains extra species, computation time will just be increased.

Eric Marcon <Eric.Marcon@ecofog.gf>

Chao, A., Chiu, C.-H. and Jost, L. (2010). Phylogenetic diversity measures based on Hill numbers. *Philosophical Transactions of the Royal Society B* 365(1558): 3599-609.

Marcon, E., Herault, B. (2015). Decomposing Phylodiversity. *Methods in Ecology and Evolution* 6(3): 333-339.

Pavoine, S., Love, M. S. and Bonsall, M. B. (2009). Hierarchical partitioning of evolutionary and ecological patterns in the organization of phylogenetically-structured species assemblages: Application to rockfish (genus: Sebastes) in the Southern California Bight. *Ecology Letters* 12(9): 898-908.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 | ```
# Load Paracou data (number of trees per species in two 1-ha plot of a tropical forest
# and their taxonomy)
data(Paracou618)
# Ps is the vector of probabilities
Ps <- as.ProbaVector(Paracou618.MC$Ps)
# Calculate the phylogenetic Shannon diversity of the plot
summary(PhyloDiversity(Ps, 1, Paracou618.Taxonomy) -> d)
plot(d)
# Ns is the vector of abundances of the metacommunity
Ns <- as.AbdVector(Paracou618.MC$Ns)
# Calculate the phylogenetic Shannon diversity of the plot
summary(bcPhyloDiversity(Ns, 1, Paracou618.Taxonomy, Correction = "Best") -> d)
plot(d)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.