Calculate net reproductive value (R0) from a matrix population model
Calculate net reproductive value (R0) from a matrix population model. The net reproduction value (R0) is the mean number of recruits produced during the mean life expectancy of an individual. See section 5.3.5 of Caswell (2001).
net_repro_rate(matU, matR, start = 1, method = "generation")
matU |
The survival component of a matrix population model (i.e. a square projection matrix reflecting survival-related transitions; e.g. progression, stasis, and retrogression). Optionally with named rows and columns indicating the corresponding life stage names. |
matR |
The reproductive component of a matrix population model (i.e. a square projection matrix only reflecting transitions due to reproduction; either sexual, clonal, or both). Optionally with named rows and columns indicating the corresponding life stage names. |
start |
Index (or stage name) of the first stage at which the author
considers the beginning of life. Only used if |
method |
The method used to calculate net reproductive value, either
|
The method argument controls how net reproductive rate is calculated.
If method = "generation", net reproductive value is calculated as the
per-generation population growth rate (i.e. the dominant eigenvalue of
matR %*% N, where N is the fundamental matrix). See Caswell
(2001) Section 5.3.4.
If method = "start", net reproductive value is calculated as the
expected lifetime production of offspring that start life in stage
start, by an individual also starting life in stage start (i.e.
(matR %*% N)[start,start]).
If offspring only arise in stage start, the two methods give the
same result.
Returns the net reproductive value. If matU is singular (often
indicating infinite life expectancy), returns NA.
Roberto Salguero-Gomez <rob.salguero@zoo.ox.ac.uk>
Hal Caswell <h.caswell@uva.nl>
Caswell, H. 2001. Matrix Population Models: Construction, Analysis, and Interpretation. Sinauer Associates; 2nd edition. ISBN: 978-0878930968
Other life history traits:
entropy_d(),
entropy_k(),
gen_time(),
life_expect_mean(),
longevity(),
repro_maturity,
shape_rep(),
shape_surv()
data(mpm1) net_repro_rate(mpm1$matU, mpm1$matF) # calculate R0 using the start method, specifying either the life stage index # or name net_repro_rate(mpm1$matU, mpm1$matF, method = "start", start = 2) net_repro_rate(mpm1$matU, mpm1$matF, method = "start", start = "small")
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.