Population Dynamics Model with Age Classes for an insect Exactly the same model as population.age.model, but written as an ordinary differential equation system (ode) with deSolve package. 7 states variables E : egg stage. homogenous population (density) (number per ha) L1 : larvae1 stage. homogenous population (density) (number per ha) L2 : larvae2 stage. homogenous population (density) (number per ha) L3 : larvae3 stage. homogenous population (density) (number per ha) L4 : larvae4 stage. homogenous population (density) (number per ha) P : pupae stage. homogenous population (density) (number per ha) A : adult stage. homogenous population (density) (number per ha)
population.age.model.ode(rb = 3.5, mE = 0.017, rE = 0.172,
m1 = 0.06, r12 = 0.217, m2 = 0.032, r23 = 0.313, m3 = 0.022,
r34 = 0.222, m4 = 0.02, r4P = 0.135, mP = 0.02, rPA = 0.099,
mA = 0.027, iA = 0, duration = 100, dt = 1, method = "euler")
: eggs laid per adult per unit area (day-1)
: relative mortality rate of egg (day-1)
: eggs hatch (day-1)
: relative mortality rate of larvae L1 (day-1)
: relative rate L1->L2 (day-1)
: relative mortality rate of larvae L2 (day-1)
: relative rate L2->L3 (day-1)
: relative mortality rate of larvae L3 (day-1)
: relative rate L3->L4 (day-1)
: relative mortality rate of larvae L4 (day-1)
: relative rate L4->P (day-1)
: relative mortality rate of purpae (day-1)
: relative rate P->A (day-1)
: relative mortality rate of adult L1 (day-1)
: input rate of adult (unit.day-1)
: simulation duration
: time step for integration
: integration method (euler, rk4,...)
data.frame with values for state variables for each time step.