# Load data
library(sf)
com <- st_read(system.file("metroparis.gpkg", package = "MTA"), layer = "com", quiet = TRUE)
# general absolute deviation
com$gdevabs <- gdev(x = com, var1 = "INC",var2 = "TH", type = "abs")
# Territorial absolute deviation calculation
com$mdevabs <- tdev(x = com, var1 = "INC", var2 = "TH", type = "abs",
key = "LIBEPT")
# Local absolute deviation calculation redistribution
com$ldevabs <- sdev(x = com, xid = "DEPCOM", var1 = "INC", var2 = "TH",
order = 1, type = "abs")
# Compute the synthesis DataFrame (absolute deviations)
mas(x = com,
gdevabs = "gdevabs",
tdevabs = "mdevabs",
sdevabs = "ldevabs",
num = "INC")
Run the code above in your browser using DataLab