Scaling Relationships

Exploring the PfEIR and the PfPR

Load the required packages:

library(ramp.xds)
library(ramp.library)
library(ramp.work)
library(deSolve)
library(rootSolve)
library(graphics)

xde_scaling

The function xde_scaling() defines the relationship between the EIR and the PR, and it outputs stable orbits for each value of aEIR in a mesh from \(10^{-1}\) up to \(10^{3}\)

xds_setup_cohort(Xname = "SIS") -> sis

xde_scaling_eir(sis, 25) -> sis

plot_eirpr(sis)