Estou tentando criar um mapa temático mostrando todos os 50 estados dos EUA, mas estou tendo problemas para realocar o Alasca e o Havaí de maneira confiável. Eu tenho algumas idéias, mas nenhuma delas funciona bem. Eu vou demonstrá-los agora.

Primeiro precisamos importar os dados; usando os dados nomaps pacote não é suficiente porque não inclui o Havaí e o Alasca.

setwd(tempdir())
download.file("https://dl.dropbox.com/s/wl0z5rpygtowqbf/states_21basic.zip?dl=1", 
              "usmapdata.zip", 
              method = "curl")
# This is a mirror of http://www.arcgis.com/home/item.html?
# id=f7f805eb65eb4ab787a0a3e1116ca7e5
unzip("usmapdata.zip")

require(rgdal)
all_states <- readOGR("states_21basic/", "states")

require(ggplot2); require(maptools); require(rgeos); require(mapproj);
all_states <- fortify(all_states, region = "STATE_NAME")

Agora nós definimos algumas estéticas da trama:

p <- ggplot() + geom_polygon( 
  aes(x=long, y=lat, group = group, fill = as.numeric(as.factor(id))), 
  colour="white", size = 0.25
) + coord_map(projection="azequalarea") + 
scale_fill_gradient(limits = c(1,50))

Agora, removemos todo o segundo plano, etc., para que não entrem em conflito quando nos sobrepormos aos estados não contíguos:

p <-   p + theme(axis.line=element_blank(),
            axis.text.x=element_blank(),
            axis.text.y=element_blank(),
            axis.ticks=element_blank(),
            axis.title.x=element_blank(),
            axis.title.y=element_blank(),
            panel.background=element_blank(),
            panel.border=element_blank(),
            panel.grid.major=element_blank(),
            panel.grid.minor=element_blank(),
            plot.background=element_blank())

Usando viewports

A primeira ideia que tive foi usar viewports:

AK <- p %+% subset(all_states, id == "Alaska") + theme(legend.position = "none")
HI <- p %+% subset(all_states, id == "Hawaii") + theme(legend.position = "none")
contiguous <- p %+% subset(all_states, id != "Alaska" & id != "Hawaii")

grid.newpage()
vp <- viewport(width = 1, height = 1)
print(contiguous, vp = vp)
subvp1 <- viewport(width = 0.25, height = 0.25, x = 0.18, y = 0.33)
print(AK, vp = subvp1)
subvp2 <- viewport(width = 0.12, height = 0.12, x = 0.32, y = 0.27)
print(HI, vp = subvp2)

Isso parece bom, mas não é satisfatório, pois é muito sensível a pequenas alterações na figura, por exemplo, redimensionamento ou alterações no tamanho e na forma da legenda.

Movendo manualmente o Alasca e o Havaí

all_states_AKHImoved <- within(all_states, {
  lat[id == "Alaska"] <- lat[id == "Alaska"] - 45
  long[id == "Alaska"] <- long[id == "Alaska"] + 40
  lat[id == "Hawaii"] <- lat[id == "Hawaii"] + 0
  long[id == "Hawaii"] <- long[id == "Hawaii"] + 70
})
p %+% all_states_AKHImoved

Isso não é satisfatório porque o Alasca geralmente não é escalável na maioria dos mapas dos EUA, então parece muito grande. Também a mudança do Alasca e do Havaí altera a distorção introduzida pela projeção do mapa.

Questão

Alguém tem alguma abordagem melhor?