shapefile to produce a linked micromap in R

问题

I am trying to replicate the example provided in the micromap R package of a linked map with the poverty rate for the 50 US states. I am trying to achieve the same results with a Shapefile from Germany. Here is my code:

library(micromap)
library(ggplot2)

setwd ('C:/Users/Jesus/Dropbox/linked maps/Chapter4')
ger<-readShapePoly("germany3.shp")
edPov<-read.csv('gerpoverty.csv')

statePolys <- create_map_table(ger, 'VARNAME_1')
head(statePolys)
lmplot(stat.data=edPov, 
       map.data=statePolys, 
       panel.types=c('labels', 'dot','map'),
       panel.data=list('Id1','poverty',NA),
       ord.by='poverty',
       grouping=5, median.row=F,
       map.link=c('Id1','VARNAME_1'))

When I tried to parse it into a flat table for use with ggplot2 the procedure does not work (instead of having 16 observations it has 8000 plus observations). Also when I try to plot the linked map, the following error appears on screen:

Error in `[.data.frame`(DF, , ord.by) : undefined columns selected

Here is the link to the files in a zip file:

https://www.dropbox.com/s/c43k755aadvu2z6/germany.rar

Any ideas or suggestions?

Thanks,

回答1:

The problem is that create_map_table(...) takes whatever you give it as the id column (VARNAME_1 in your case), and puts that in a column called ID in the output (statePolys in your case). So now there is a statePolys$ID column with whatever you had in VARNAME_1. When you tell lmplot(...) to join on edPov$Id1 and statePolys$VARNAME_1, it can't find the latter column. Another way to say this is that the second element of map.link must always be ID. So this code works (I used the ID_1 column in both edPov and ger here, since that seems to contain a German State ID).

library(micromap)
library(ggplot2)

ger<-readShapePoly("germany3.shp")
edPov<-read.csv('gerpoverty.csv')

statePolys <- create_map_table(ger, 'ID_1')  # ID_1 stored in statePolys$ID
head(statePolys)
lmplot(stat.data=edPov, 
       map.data=statePolys, 
       panel.types=c('labels', 'dot','map'),
       panel.data=list('Id1','poverty',NA),
       ord.by='poverty',
       grouping=5, median.row=F,
       map.link=c('ID_1','ID'))

Here, ID_1 refers to edPoly$ID_1, and ID refers to statePolys$ID.

来源：https://stackoverflow.com/questions/21651985/shapefile-to-produce-a-linked-micromap-in-r

标签

maps

gis