Learn R Programming

entropart (version 1.1.3)

MetaCommunity: Metacommunity class

Description

Methods for objects of type "MetaCommunity".

Usage

MetaCommunity(Abundances, Weights = rep(1, ncol(Abundances)))
is.MetaCommunity(x)
## S3 method for class 'MetaCommunity':
summary(object, \dots)
## S3 method for class 'MetaCommunity':
plot(x, \dots)

Arguments

Abundances
A dataframe containing the number of observations (lines are species, columns are communities). The first column of the dataframe may contain the species names.
Weights
A vector of positive numbers equal to community weights or a dataframe containing a vector named Weights. It does not have to be normalized. Weights are equal by default.
x
An object to be tested or plotted.
object
A MetaCommunity object to be summarized.
...
Additional arguments to be passed to the generic methods.

Value

  • An object of class MetaCommunity is a list:
  • NsiA matrix containing abundance data, species in line, communities in column.
  • NsA vector containing the number of individuals of each species.
  • NiA vector containing the number of individuals of each community.
  • NThe total number of individuals.
  • PsiA matrix whose columns are the probability vectors of communities (each of them sums to 1).
  • WiA vector containing the normalized community weights (sum to 1).
  • PsA vector containing the probability vector of the metacommunity.
  • NspeciesThe number of species.
  • NcommunitiesThe number of communities.
  • SampleCoverageThe sample coverage of the metacommunity.
  • SampleCoverage.communitiesA vector containing the sample coverages of each community.
  • is.MetaCommunity returns TRUE if the object is of class MetaCommunity. summary.MetaCommunity returns a summary of the object's value. plot.MetaCommunity plots it.

Details

In the entropart package, individuals of different "species" are counted in several "communities" which are agregated to define a "metacommunity". This is a naming convention, which may correspond to plots in a forest inventory or any data organized the same way. Alpha and beta entropies of communities are summed according to Weights and the probability to find a species in the metacommunity is the weighted average of probabilities in communities. The simplest way to import data is to organize it into two text files. The first file should contain abundance data: the first column named Species for species names, and a column for each community. lcc{ Species NameOfCommunity1 NameOfCommunity2 NameOfSpecies1 1 5 NameOfSpecies2 4 2 ... ... ... } The second file should contain the community weights and be organized as follows: lc{ Communities Weights NameOfCommunity1 3 NameOfCommunity2 1 } Files can be read and data imported by code such as: Abundances <- read.csv(file="Abundances.csv") Weights <- read.csv(file="Weights.csv") MC <- MetaCommunity(Abundances, Weights)

Examples

Run this code
# Use BCI data from vegan package
  if (require(vegan, quietly = TRUE)) {
    # Load BCI data (number of trees per species in each 1-ha plot of a tropical forest)
    data(BCI)
    # BCI dataframe must be transposed (its lines are plots, not species)
    BCI.df <- as.data.frame(t(BCI))
    # Create a metacommunity object from a matrix of abundances and a vector of weights
    # (here, all plots have a weight equal to 1)
    MC <- MetaCommunity(BCI.df, rep(1,50))
  }

Run the code above in your browser using DataLab