flows (version 1.1.1)

plotMapDomFlows: Dominant Flows Map

Description

This function plots a dominant flows map.

Usage

plotMapDomFlows(mat, spdf, spdfid, w, wid, wvar, wcex = 0.05, legend.flows.pos = "topright", legend.flows.title = "flow intensity", legend.nodes.pos = "topleft", legend.node.txt = c("Dominant", "Intermediary", "Dominated", "Size proportional\nto sum of inflows"), add = FALSE)

Arguments

mat

A square matrix of dominant flows (see domflows).

spdf

A SpatialPolygonsDataFrame or a SpatialPointsDataFrame of units.

spdfid

Name of the unique identifier variable in the spdf data.frame.

w

A data.frame which contains the weight variable used to plot units sizes on the map.

wid

Name of the unique identifier variable in w.

wvar

Name of the weight variable in w.

wcex

Share of the surface of the map occupied by circles (0.02 is 2%).

legend.flows.pos

Position of the flows legend, one of "topleft", "top", "topright", "left", "right", "bottomleft", "bottom", "bottomright".

legend.flows.title

Title of the flows legend.

legend.nodes.pos

Position of the nodes legend, one of "topleft", "top", "topright", "left", "right", "bottomleft", "bottom", "bottomright".

legend.node.txt

Text of the nodes legend.

add

A boolean, if TRUE, add the layer to an existing plot.

See Also

domflows, plotDomFlows

Examples

# Import data
data(nav)
myflows <- prepflows(mat = nav, i = "i", j = "j", fij = "fij")

# Remove the matrix diagonal
diag(myflows) <- 0

# Select the dominant flows (incoming flows criterion)
flowSel1 <- domflows(mat = myflows, w = colSums(myflows), k = 1)
# Select the first flows
flowSel2 <- firstflows(mat = myflows, method = "nfirst", ties.method = "first",
                       k = 1)

# Combine selections
flowSel <- myflows * flowSel1 * flowSel2

# Node weights
inflows <- data.frame(id = colnames(myflows), w = colSums(myflows))

# Plot dominant flows map
opar <- par(mar = c(0,0,2,0))
sp::plot(GE, col = "#cceae7", border = NA)
plotMapDomFlows(mat = flowSel, spdf = UA, spdfid = "ID", w = inflows, wid = "id",
                wvar = "w", wcex = 0.05, add = TRUE,
                legend.flows.pos = "bottomleft",
                legend.flows.title = "Nb. of commuters")
title("Dominant Flows of Commuters")
mtext(text = "INSEE, 2011", side = 4, line = -1, adj = 0.01, cex = 0.8)
par(opar)