calc.decay.stats: Calculate model summary and validation statistics for trip...
In gilesjohnr/hmob: Tools for analyzing human mobility data

Description Usage Arguments Value Author(s) See Also Examples

This function takes the output from a the jags.model function and calculates summaries of the posterior parameter estimates, convergence diagnostics, and sample sizes for each route.

1	calc.decay.stats(x, y.route, districts, pop.dens, dist.mat)

`x`	JAGS model object
`y.route`	matrix of route-level trip duration counts
`districts`	vector of district names
`pop.dens`	vector of district population densites
`dist.mat`	distance matrix

a list containing a dataframe of parameter estimates and diagnostics, and a matrix of actual and predicted trip duration counts

John Giles

Other model processing: get.param.vals(), mob.data.array.pop.level(), mob.data.array.route.level()

load('./data/duration_data_arrays_5day_full.Rdata')                    # y.month, y.route, y.pop
load('./MARCC/output/decay_model_route_hierarchical_5day_full.Rdata')  # out # Decay model parameters (Lambda)
load('./data/distance_matrix_named.rdata')                             # distance.matrix
pop <- read.csv("./data/population_data.csv")                          # population densiity

districts <- dimnames(y.route)$origin

dens <- vector()
for (i in seq_along(districts)) {
     dens <- c(dens, pop$density[as.numeric(which(pop$ID_2 == districts[i]))])
}

mod.decay <- calc.decay.stats(x=out, 
                              y.route=y.route, 
                              districts=districts, 
                              pop.dens=dens, 
                              dist.mat=distance.matrix[districts, districts])