R: creazione di una mappa della rete mondiale

Question

Vorrei produrre una mappa della rete mondiale in R simile a the one below from this page.

Sono stato alla ricerca di un pacchetto R che mi permetterà di fare questo, ma non sono stato in grado di trovare uno. C'è D3 JavaScript Network Graphs from R ma non sono riuscito a trovare un esempio di mappa della rete mondiale.
Come posso creare qualcosa di simile in R?

Answer 1

Bene, FWIW: ecco un breve & modo sporco per tracciare i vertici ("città") su una mappa e collegare i bordi tra i vertici con le frecce:

library(maps) 
library(diagram) 
library(plotrix) 
palette(rainbow(20)) 
data("world.cities") 

pdf(tf <- tempfile(fileext = ".pdf"), width = 40, height = 20) 

map('world', fill = TRUE, col = "lightgray", mar = rep(0, 4)) 

nodes <- transform(with(world.cities, world.cities[pop > 5e6,]), country.etc = as.factor(country.etc)) 
with(nodes, points(long, lat, col=country.etc, pch=19, cex=rescale(pop, c(1, 8)))) 

set.seed(1) 
edges <- subset(data.frame(from = sample(nodes$name, 20, replace = TRUE), to = sample(nodes$name, 20, replace = TRUE), stringsAsFactors = F), from != to) 
edges <- merge(merge(edges, nodes[, c("name", "long", "lat")], by.x = "from", by.y = "name"), nodes[, c("name", "long", "lat")], by.x = "to", by.y = "name") 
edges$col <- as.integer(nodes$country.etc[match(edges$from, nodes$name)]) 

apply(edges[, -(1:2)], 1, function(x) curvedarrow(to=x[3:4], from=x[1:2], lcol=x[5], curve=.1, arr.pos = 1, lwd=.5)) 

dev.off()   
shell.exec(tf) 

Answer 2

Ecco una soluzione utilizzando il pacchetto geosphere e maps. Utilizzando la funzione gcIntermediate è possibile "punti intermedi su un grande cerchio tra i due punti utilizzati per definire il grande cerchio".

Ecco un esempio che mostra conncetions viaggio dall'aeroporto JKF (dati di esempio dal pacchetto nycflights13 filtrata con dplyr):

library(maps) 
library(geosphere) 
library(dplyr) 
library(nycflights13) 


usairports <- filter(airports, lat < 48.5) 
usairports <- filter(usairports, lon > -130) 
usairports <- filter(usairports, faa!="JFK") 
jfk <- filter(airports, faa=="JFK") 

map("world", regions=c("usa"), fill=T, col="grey8", bg="grey15", ylim=c(21.0,50.0), xlim=c(-130.0,-65.0)) 


for (i in (1:dim(usairports)[1])) { 

inter <- gcIntermediate(c(jfk$lon[1], jfk$lat[1]), c(usairports$lon[i], usairports$lat[i]), n=200) 

lines(inter, lwd=0.1, col="turquoise2")  
} 

points(usairports$lon,usairports$lat, pch=3, cex=0.1, col="chocolate1") 

Questo si basa su un tutorial pubblicato sul mio blog.