This function applies a KL-divergence approach to a list of indicator species in phyloregions.

indicators(
  theta,
  top_indicators = 5,
  method = c("poisson", "bernoulli"),
  options = c("min", "max"),
  shared = FALSE
)

Arguments

theta

A matrix or data.frame of cluster probability distributions from a topics modeling.

top_indicators

Integer to obtain the top driving species in clusters.

method

The model assumption for KL divergence measurement. Available choices are "poisson" (default) and "bernoulli".

options

Option "min" selects species that maximize the minimum KL divergence of a phyloregion vs all other phyloregions. Option "max" selects species that maximize the maximum KL divergence of a phyloregion against all other phyloregions.

shared

Logical if TRUE, lists top species driving patterns in more than one phyloregion.

Value

A list of top indicator species and their indicator values

Examples

data(africa)
indsp <- indicators(africa$theta, top_indicators = 5,
                    options = "max", method = "poisson")