SizeGamma Size Overlap Randomization Algorithm

Usage

size_gamma(speciesData = rnorm(50, mean = 100, sd = 1))

Arguments

speciesData
a vector of body sizes or other trait measurements of species. All values must be positive real numbers.

Value

Returns a vector of simulated body sizes as the same length as speciesData.

Description

Function to generate a random distribution of body sizes by drawing from a gamma distribution. Shape and rate parameters of the gamma are estimated from the empirical data.

Note

The shape and rate parameters are estimated from the real data using the maximum likelihood estimators generated from the fitdr function in the MASS library. The flexible gamma distribution can be fit to a variety of normal, log-normal, and exponential distributions that are typical for trait data measured on a continuous non-negative scale.

Examples

obsOverlap <- size_gamma(speciesData=rnorm(50,mean=100,sd=1))

See also

fitdistr in the MASS library.