Risk by epidemiological unit

Compute the risk for each of the epidemiological units.

risk_unit(r, eu, fun = mean)

Arguments

r

RasterLayer object
eu

SpatialPolygons* object
fun

A summary function. By default the mean.

Value

numeric vector of the same length as the number of epidemiological units.

Details

Summarises the level of risk using a function such as the mean or median value across the epidemiological unit.

Examples

  cmr <- mapMCDA_datasets()
  cmr$cmr_admin3$rv <- risk_unit(cmr$animal.density, cmr$cmr_admin3)
#> Warning: CRS object has comment, which is lost in output
  sp::spplot(cmr$cmr_admin3[, "rv"], cuts = 3)