calc.wait.time: Calculate aggregate waiting time distribution

In gilesjohnr/hmob: Tools for analyzing human mobility data

Description

This function aggregates the all simulated waiting time distributions into a single probability density for each district. The method uses simple linear pooling to get the weighted linear average of all waiting time probabilities (sometimes referred to as the 'linear opinion pool').

Usage

calc.wait.time(x)

Arguments

x	a three-dimensional array containing waiting time distributions for all districts and simulations

Value

a matrix (rows = districts, cols = epidemic generation)

Author(s)

John Giles

See Also

Other simulation: calc.hpd(), calc.prop.inf(), calc.prop.remain(), calc.timing.magnitude(), decay.func(), get.age.beta(), get.beta.params(), sim.TSIR.full(), sim.TSIR(), sim.combine.dual(), sim.combine(), sim.gravity.duration(), sim.gravity(), sim.lambda(), sim.pi(), sim.rho(), sim.tau()

Examples

load("./output/TSIR_sim_flu_high.Rdata")

x1 <- sim$B$wait.time
x2 <- sim$R$wait.time

agg1 <- calc.wait.time(x1)
agg2 <- calc.wait.time(x2)

plot(agg1[1,], type='l', col='blue', cex=2)
lines(agg2[1,], col='red', cex=2)

gilesjohnr/hmob documentation built on Aug. 8, 2020, 1:26 a.m.