Above are all the numbers you change for Day‑1 and Day‑2. Everything below is read‑only unless you want to explore further!
The chunk below reproduces your two‑sex, 8‑class model from
Population dynamic model.R. It is collapsed and eval=FALSE
so the default knit ignores it. Feel free to enable and explore!
# ---- Two-sex 8‑class model (excerpt) ----
# 1. parameters
ages <- c("JuvF","SubF","AdultF","OldF","JuvM","SubM","AdultM","OldM")
S <- c(0.55,0.80,0.88,0.75, 0.50,0.70,0.75,0.70)
Fec <- c(0,0,1.25,0, 0,0,0,0) # females only
# 2. build 8×8 Leslie & project...
# (paste your existing loop or %*% calculation here)
# 3. plot total females, males, harvest vector, etc.