R/examples/calc_hpd.R

In gilesjohnr/hmob: Tools for analyzing human mobility data

load("./output/simulations/13/influenza.Rdata")

x <- sim$B$wait.time   # basic gravity model
#x <- sim$R$wait.time   # gravity model with duration 

agg <- calc.wait.time(x)
hpd <- calc.hpd(agg, ci=c(0.5, 0.95))

par(mfrow=c(3,4))
for (i in 1:nrow(hpd)) {
     plot(agg[i,], type='l', main=paste("district ", hpd$district[i]), ylab='Waiting time density')
     abline(v=hpd$lo95[i], lty=2, col='red')
     abline(v=hpd$hi95[i], lty=2, col='red')
     abline(v=hpd$max[i], col='red')
}