inst/examples/psi.R

In mattmar/dynamAedes: A Unified Mechanistic Model for the Population Dynamics of Invasive Aedes Mosquitoes

\donttest{
# Make a toy temperature time series
w <- matrix(seq(20,25,length.out=5),ncol=5)*1000

# Run the model
	mod <- dynamAedes.m(
	species="koreicus", 
	scale="ws",
	intro.eggs=10, 
	jhwv=2, 
	temps.matrix=w, 
	startd="2021-06-21", 
	endd="2021-06-25",
	lat=42,
	long=8,
	n.clusters=1, 
	iter=1,
	compressed.output=TRUE)

# Check proportion of viable populations on the last day of simulations
psi(input_sim = mod, eval_date = 3)
}