Rel.homog: Relative Homogeneity

In vegetarian: Jost Diversity Measures for Community Data

Description

Provides a comparison between communities with unequal weights based on the numbers equivalent of the Shannon beta diversity (Jost 2007). When community weights are not equal, the Shannon measures of diversity where q=1, are still valid and can be used to calculate a relative measure of homogeneity derived from MacArthur's homogeneity measure (Equation 22 from Jost 2007).

Usage

Rel.homog(abundances, abundances2 = NULL, wts = FALSE, boot = FALSE, boot.arg = list(s.sizes = NULL, num.iter = 100))

Arguments

abundances	Community data as a matrix where columns are individual species and rows are sites or a vector of different species within a site. Matrix and vector elements are abundance data (e.g. counts, percent cover estimates).
abundances2	Community data, a vector of different species within a site. Vector elements are abundance data (e.g. counts, percent cover estimates). If abundances is given a matrix, then abundances2 defaults to NULL.
wts	A vector of community weights whose length equals the number of communities being compared. See details below for examples of when community weights might be treated as equal or unequal. Defaults to wts=FALSE where all communities are equally weighted.
boot	Logical indicating whether to use bootstrapping to estimate uncertainty.
boot.arg	(optional) List of arguments to pass bootstrapping function: list(s.sizes=number you specify, num.iter=number you specify)

Details

This measure assumes the Shannon case where the order q=1.

Value

REL.HOMOG	A scalar between zero (all communities distinct) and one (all communities identical).
StdErr	(optional) Standard error of value estimated through bootstrapping.

Author(s)

Noah Charney, Sydne Record

References

Jost, L. 2007. Partitioning diversity into independent alpha and beta components. Ecology 88(10): 2427-2439.

See Also

M.homog, similarity, turnover.

Examples

data(simesants)
Rel.homog(simesants[1:2,-1])
hemlock<-subset(simesants,Habitat=="Hemlock")[,-1]
hardwood<-subset(simesants,Habitat=="Hardwood")[,-1]
Rel.homog(abundances=hemlock,abundances2=hardwood)
Rel.homog(simesants[1:2,-1], boot=TRUE)

Example output

[1] 0.8464179
[1] 0.8464179
$REL.HOMOG
[1] 0.8464179

$StdErr
[1] 0.02597221

vegetarian documentation built on May 1, 2019, 8:01 p.m.