sim=c("min", "mean", "max") option for plot.netsim for type = "network"

[smjenness]

Options would be sim=c("min", "mean", "max")

This would be able to automate this type of network plotting:

low <- which.min(sim1$epi$i.num[500, ])
high <- which.max(sim1$epi$i.num[500, ])

par(mfrow = c(1, 2))
plot(sim1, type = "network", col.status = TRUE, 
     at = 500, sim = low)
plot(sim1, type = "network", col.status = TRUE, 
     at = 500, sim = high)

[smjenness]

The day 4 NME tutorial should be updated accordingly with this

[smjenness]

This can actually be in the as.data.frame method for both ICM and network models, then the plot option can easily call this for the sim argument. UPDATE: nevermind, these calculations need to be time-specific, whereas as.data.frame is not.

[smjenness]

Also need to update Day 5 tutorial 1:

get_msim <- function(x, at) {
  which.min(abs(as.numeric(x$epi$i.num[at, ]) - 
              mean(as.numeric(x$epi$i.num[at, ]))))
}
med.si1 <- get_msim(si1, 500)
med.si2 <- get_msim(si2, 500)