Estimate the distribution of final outbreak sizes by group using stochastic simulations of multi-group model

Usage

getFinalSizeDist(n, transmrates, recoveryrate, popsize, initR, initI, initV)

Arguments

n: the number of simulations to run
transmrates: matrix of group-to-group (column-to-row) transmission rates
recoveryrate: inverse of mean infectious period
popsize: the population size of each group
initR: initial number of each group already infected and removed for each simulation (included in final size result); if a matrix, row i is used for simulation i
initI: initial number of each group infectious for each simulation; if a matrix, row i is used for simulation i
initV: initial number of each group immune due to vaccination or prior-outbreak infection for each simulation (not included in final size result); if a matrix, row i is used for simulation i

Value

a matrix with the final number infected from each group (column) in each simulation (row)

Examples

getFinalSizeDist(n = 10, transmrates = matrix(0.2, 2 ,2), recoveryrate = 0.3,
popsize = c(100, 150), initR = c(0, 0), initI = c(0, 1), initV = c(10, 10))
#>       [,1] [,2]
#>  [1,]    0    1
#>  [2,]    0    1
#>  [3,]   27   40
#>  [4,]   30   41
#>  [5,]   26   47
#>  [6,]    0    1
#>  [7,]    0    1
#>  [8,]   41   62
#>  [9,]    0    1
#> [10,]   40   68