Pre- and Post-reproduction census
p0 <- 0.2
p1 <- 0.9
p2 <- 0.6
m2 <- 3.0
m3 <- 6.0
A1 <- matrix(c(0.0, m2 * p1, m3 * p2,
p0, 0.0, 0.0,
0.0,p1,0), byrow = TRUE, nrow = 3)
Convert to pre-breeding
A2 <- matrix(c(0.0, m2 * p0, m3 * p0,
p1, 0.0, 0.0,
0.0,p2,0), byrow = TRUE, nrow = 3)
Compare population growth rates
(popdemo::eigs(A1, what = "lambda"))
## [1] 1.070261
(popdemo::eigs(A2, what = "lambda"))
## [1] 1.070261