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