Add covariates to model

[sempwn]

For more flexible modeling, would be useful to include a covariate matrix for locations and predictors that can be used to explain variance of R0 for each location.

For example, we may have difference facility types we would want to use to adjust R0 estimates, or type of ventilation, age of building, density etc.