sipd {meteR}

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, ...)

## S3 method for class 'meteESF'
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

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.

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'

Author(s)

Andy Rominger <[e-mail]>, Cory Merow

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

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