Simulate a population size state during \([1, t_{lim}]\) times. The state at time \(t+1\) is dependent on state at time \(t\) and the projection matrix (IPM). The simulation will run until \(t_{lim}\) and if the equilibrium is not reached, it will continue. Only simulation in \([1, t_{lim}]\), and equilibrium state are returned.
Usage
sim_deter_forest(
Forest,
tlim = 3000,
equil_dist = 250,
equil_diff = 1,
equil_time = 10000,
harvest = c("default", "Uneven", "Even", "Favoured_Uneven"),
targetBA = 20,
targetRDI = 0.9,
targetKg = 0.9,
final_harv = 100,
climate = NULL,
disturbance = NULL,
correction = "none",
SurfEch = 0.03,
verbose = FALSE
)Arguments
- Forest
Group of species that each contains IPM for deterministic transition for \(Z_{t}\) state in a population to \(Z_{t+1}\) state. A species is also defined with recruitment and harvest functions, please see
speciesfor more information.- tlim
Number of simulation iterations (years) in the future. single int.
- equil_dist
Number of last n time for which the range difference should not exceed
equil_diffduring the equilibrium research. single int.- equil_diff
Difference threshold of the basal area under which equilibrium is assumed. single real.
- equil_time
Total maximum time simulation allowed in equilibrium research. Must be higher or equal to tlim and equil_dist. single int.
- harvest
Choice of harvest rules between default, Uneven and Even. This indicate what settings will be used. See Details.
- targetBA
BA value per ha that is targeted when using uneven harvesting. Single numeric in \(m^2\).
- targetRDI
RDI value that is targeted when using even harvesting. RDI is the ratio between the number of trees and the maximum number of trees given the self-thinning boundary for the corresponding mean diameter and species.
- targetKg
Kg value that is targeted when using even harvesting. Kg is the ratio between mean quadratic diameter of killed trees and mean quadratic diameter of trees before harvesting.
- final_harv
Final harvest time used when
harvestis set to "Even". This parameter drives the final cut time for even stands.- climate
Optional, climate matrix if climate variation along time is needed. Climate variation rely on species created with mu_gr class objects. This matrix require as many rows as time steps until equil_time. If the climate does not variate, a single row can given and will be reused.
- disturbance
![[Experimental]](figures/lifecycle-experimental.svg)
parameter.- correction
Choice of correction of the IPM between
"none"(default) and"cut". The second option set the last column to 0 in the IPM so that no individual can grow outside of the defines classes.- SurfEch
Value of plot size surface in ha
- verbose
Print message. FALSE by default
Value
Data.frame with long tidyverse format : a row for each observation and a column per variable. Columns are listed below, some may contains NA values, as for example species when there is a non-specific variable (BA).
- species
Name of the species.
- var
Variable of interest
- time
Time step of the simulation. If the equilibrium is the last time in
tliminput, this time will occur twice in the table.- mesh
Mesh class number, from 1 to n class.
- size
Size corresponding to the mesh class.
- equil
Logical if this time step is the equilibrium or last step of simulation
- value
Numeric values of the variables.
The variables are :
- n
Distribution of density by mesh along time per ha.
- N
Sum of density per ha. (colSums for n)
- BAsp
Basal area of the population per ha and species
- BAstand
Basal area of the population per ha and species when excluding size class below dth. See Harvesting vignette.
- h
Distribution of harvest density by mesh along time per ha.
- H
Sum of harvested density per ha. (colSums for h)