Scaling Relationships

Seasonality

library(ramp.xds)
library(ramp.library)
library(ramp.work)

load("sis.rda")

sin_par1 <- makepar_F_sin(floor=0.2)
F_s1 = make_function(sin_par1)
sis_seas1 <- xds_setup_cohort(Xname = "SIS", F_season = F_s1) 

xde_scaling_eir(sis_seas1, 25) -> sis_seas1

require(viridis)

Loading required package: viridis

Loading required package: viridisLite

clrs = turbo(25)

plot_eirpr(sis, grey(0.5))
lines_eirpr(sis_seas1, "darkblue")

with(sis_seas1$output$eirpr,{
  points(aeir, pr, col = clrs)
  lines(scaling[[5]]$aeir, scaling[[5]]$pr, col = clrs[5])
  lines(scaling[[10]]$aeir, scaling[[10]]$pr, col = clrs[10])
  lines(scaling[[15]]$aeir, scaling[[15]]$pr, col = clrs[15])
  lines(scaling[[20]]$aeir, scaling[[20]]$pr, col = clrs[20])
})