SeedProduction {comsimitv}R Documentation

Calculating number of produced seeds

Description

Number of seeds calculated following the formula used by Botta-Dukat & Czucz (2016). This built-in function can be replaced by a user-defined one.

Usage

SeedProduction(compet, b0 = 1, K = 200, seed.distrib = c("pois", "binom"), ...)

Arguments

compet

Matrix of pairwise competition coefficients

b0

Probability of producing seed, if no competition

K

Critical level of competition (See Details)

seed.distrib

Distribution of seed numbers (See Details)

...

any additional parameters

Details

Expected value of produced seeds is a decreasing sigmoid function of strength of competition (sum of abundances weighted by competition coefficients). If strength of competition is higher than parameter K, probability is set to zero. See vignette("competition") for formulas

In simulation of Botta-Dukat & Czucz (2016) each individual produces one seed or does not produce seed at all. In this case number of seeds follows binomial distribution (i.e. distrib="binom"). A more realistic alternative is using Poisson distribution (distrib="pois").

Value

Matrix of produced seeds


[Package comsimitv version 0.1.5 Index]