inst/examples/adci.R

In 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=100, 
	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)

# Derive interquarile of population egg abundance per cell
adci(mod, eval_date=1:3, breaks=c(0.25,0.50,0.75), st=1)
}