From 1b4a959701ba054fb1bb7d30d1d8204ea13dffd7 Mon Sep 17 00:00:00 2001 From: Paige Miller Date: Tue, 26 Nov 2024 19:57:33 +0000 Subject: [PATCH] updates --- scratch/ngm-structured-sir.R | 6 ++---- 1 file changed, 2 insertions(+), 4 deletions(-) diff --git a/scratch/ngm-structured-sir.R b/scratch/ngm-structured-sir.R index bf197d3..cc33bf6 100644 --- a/scratch/ngm-structured-sir.R +++ b/scratch/ngm-structured-sir.R @@ -8,8 +8,8 @@ ngm_sir <- function(N, V, K, VE, p_severe) { K <- K * S / N eigenvalues <- eigen(K) - r_effective <- max(eigenvalues$values) - infections <- eigenvalues$vectors[, which.max(eigenvalues$values)] * N + r_effective <- max(abs(eigenvalues$values)) + infections <- eigenvalues$vectors[, which.max(abs(eigenvalues$values))] severe_infections <- infections * p_severe return(list(r_e = r_effective, infections = infections, severe_infections = severe_infections)) } @@ -22,8 +22,6 @@ K <- matrix(c( 1, 1, 1, 1 # general ), nrow = 4, ncol = 4) - - N <- c(100, 100, 10, 790) # pop size: core, kids, travelers, general V <- c(10, 0, 0, 0) # doses VE <- 0.7 # vaccine efficacy