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.