MTA (version 0.2.0)

mst: Multiscalar Typology

Description

Compute a multiscalar typology according to the three relative deviations (general: G, territorial: T and spatial: S). The elementary units are classified in eight classes according to their three relative positions.

Usage

mst(x, var1, var2, ref, key, order, dist, mat, xid, threshold,
  superior = FALSE)

Arguments

x

a sf object or a SpatialPolygonsDataFrame including var1 and var2.

var1

name of the numerator variable in x.

var2

name of the denominator variable in x.

ref

ratio of reference; if missing, the ratio of reference is the one of the whole study area (sum(var1) / sum(var2)).

key

aggregation key field for measuring the deviation (intermediate territorial level).

order

contiguity order.

dist

distance threshold defining the contiguity. The cartesian distance between units centroids is used by default ; use mat to apply different metrics.

mat

a distance matrix (road distance, travel time...) between x units. Row and column names must fit xid identifiers. (optional)

xid

identifier field in x (to be used for importing a personal distance matrix). Default to the first column.

threshold

defined to build the typology (100 is considered as the average).

superior

if TRUE, deviation values must be greater than threshold. If FALSE, deviation values must be lower than threshold.

Value

A dataframe including the ratio (var1/var2), the 3 relative deviations (G, T and S) and the resulting ordered typology (0 to 7).

• 0: none

• 1: G

• 2: T

• 3: G and T

• 4: S

• 5: G and S

• 6: T and S

• 7: G, T and S

Examples

# NOT RUN {
# load data
data("GrandParisMetropole")
# Multiscalar typology - Wealthiest territories  
row.names(com) <- com$LIBCOM
synthesis <- mst(x = com, var1 = "INC", var2 = "TH", key = "EPT",
                 order = 1, threshold = 125, superior = TRUE)
# Territories avbove 125% for the 3 deviations 
synthesis[synthesis$mst == 7, ]

# Multiscalar typology - Poorest territories
synthesis <- mst (x = com, var1 = "INC", var2 = "TH", key = "EPT",
                  order = 1, threshold = 75, superior = FALSE)
# Territories below 75 % for the three deviations
synthesis[synthesis$mst == 7, ]
# }