Learn R Programming

meteR (version 1.2)

sipd: Generic method to obtain the species-level individual power distribution (SIPD)

Description

Extract species level individual power distribution from ESF object and return object inheriting from meteDist. This distribution (Theta) describes the distribution of metabolic rates across the individuals of a species with n individuls

Usage

sipd(x, ...)
"sipd"(x, sppID, n, ...)

Arguments

x
An object of class meteESF (i.e. the fitted distribution $R(n,e)$)
...
arguments to be passed to methods
sppID
the name or index of the species of interest as listed in the spp argument passed to meteESF
n
integer. Alternatively can extract METE prediction by indicating number of individuals in the species

Value

An object of class meteDist. The object contains a list with the following elements.
data
The data used to construct the prediction
d
density funciton
p
cumulative density function
q
quantile funtion
r
random number generator
La
Vector of Lagrange multipliers
state.var
State variables used to constrain entropy maximization
type
Specifies the type of distribution is 'sad'

Details

If n is provided then only the theoretical prediction is returned (because data from multiple species could map to the same n). Thus if data and prediction are desired use sppID.

References

Harte, J. 2011. Maximum entropy and ecology: a theory of abundance, distribution, and energetics. Oxford University Press.

See Also

sad.meteESF, ipd.meteESF, metePsi

Other Theta: meteTheta

Examples

Run this code
data(arth)
esf1 <- meteESF(spp=arth$spp,
                abund=arth$count,
                power=arth$mass^(.75),
                minE=min(arth$mass^(.75)))
sipd1 <- sipd(esf1, sppID=5)
sipd1

Run the code above in your browser using DataLab