[R] mapping data to a geographic map of Europe

Bert Gunter gunter.berton at gene.com
Fri Nov 1 20:45:48 CET 2013


... but you may be interested in this:

http://andywoodruff.com/blog/why-are-choropleth-mercator-maps-bad-because-we-said-so/

Cheers,
Bert

On Fri, Nov 1, 2013 at 12:18 PM, Adams, Jean <jvadams at usgs.gov> wrote:
> Claudia,
>
> I have not worked through the example myself.  Since you seem to be getting
> errors, perhaps a different example would help.  Here are some more
> choropleth maps (although these use US states rather than European
> countries).
>
> http://blog.revolutionanalytics.com/2009/11/choropleth-challenge-result.html
>
> Jean
>
>
>
> On Fri, Nov 1, 2013 at 12:04 PM, <paladini at trustindata.de> wrote:
>
>> Hi Jean,
>> thanks again for your response. As I told you  I did the downloads and
>> double checked if I selected the right directory.
>> But I noticed right now what happend:
>> The command in the example is :
>>
>> eurMap <- readShapePoly(fn="NUTS_2010_**60M_SH/Shape/data/NUTS_RG_60M_**
>> 2010")
>>
>> But it should be :
>> eurMap <- readShapePoly(fn="NUTS_2010_**60M_SH/data/ggg/NUTS_RG_60M_**
>> 2010")
>>
>> Because if you do the downloads and unzip these data there is no such
>> think as a "Shape" directory.
>>
>> Now eurMap <- readShapePoly(fn="NUTS_2010_**60M_SH/data/NUTS_RG_60M_2010")
>> works.
>>
>>
>> What happens now is that after typing:
>> eurEduMapDf <- merge(eurMapDf, eurEdu, by.x="id", by.y="GEO")
>> I get another error message because "eurMapDf" is unknown.
>>
>> So I supposed it should be:
>> eurEduMapDf <- merge(eurMap, eurEdu, by.x="id", by.y="GEO")
>>
>> But this doesn't work either.The error message this time is: undefined
>> column selected.
>>
>> Did I do something wrong?
>>
>>
>>
>> Best regards
>>
>> Claudia
>>
>>
>>
>>
>>
>>
>>
>>
>> Zitat von "Adams, Jean" <jvadams at usgs.gov>:
>>
>>  Claudia,
>>>
>>> You should cc r-help on all correspondence so that others can follow the
>>> thread.
>>>
>>> In the second paragraph of the link I sent you
>>>      http://www.r-bloggers.com/**maps-in-r-choropleth-maps/<http://www.r-bloggers.com/maps-in-r-choropleth-maps/>
>>> a link is provided for the NUTS data,
>>>      "The polygons for drawing the administrative boundaries were obtained
>>> from this link. In particular, the NUTS 2010 shapefile in the 1:60 million
>>> scale was downloaded and used. The other available scales would allow the
>>> drawing of better defined maps, but at a computational cost. The zipped
>>> file has to be extracted in a folder of choice for using it later."
>>>
>>> http://epp.eurostat.ec.europa.**eu/portal/page/portal/gisco_**
>>> Geographical_information_maps/**popups/references/**administrative_units_
>>> **statistical_units_1<http://epp.eurostat.ec.europa.eu/portal/page/portal/gisco_Geographical_information_maps/popups/references/administrative_units_statistical_units_1>
>>>
>>> If you want to follow the example, you will need to download this data to
>>> your computer and then make sure that you refer to the appropriate
>>> directory when using the readShapePoly() function.
>>>
>>> Jean
>>>
>>>
>>>
>>> On Wed, Oct 30, 2013 at 11:14 AM, <paladini at trustindata.de> wrote:
>>>
>>>  Hi Jean,
>>>> thank you for your advice.
>>>> The page looks quite interesting and I tried the example in  GNU R. I did
>>>> all the downloads.
>>>> But
>>>> Just in the beginnig after typing
>>>>
>>>> eurMap <- readShapePoly(fn="NUTS_2010_****60M_SH/Shape/data/NUTS_RG_60M_
>>>> ****
>>>>
>>>> 2010")
>>>> I get the following error message:
>>>>
>>>> Error in getinfo.shape(filen) : Error opening SHP file
>>>>
>>>> To you have an idea what I did wrong?
>>>>
>>>> Thanks a lot and best regards
>>>>
>>>> Claudia
>>>>
>>>>
>>>>
>>>>
>>>> Zitat von "Adams, Jean" <jvadams at usgs.gov>:
>>>>
>>>>  Check out this link for some examples
>>>>
>>>>>     http://www.r-bloggers.com/****maps-in-r-choropleth-maps/<http://www.r-bloggers.com/**maps-in-r-choropleth-maps/>
>>>>> <htt**p://www.r-bloggers.com/maps-**in-r-choropleth-maps/<http://www.r-bloggers.com/maps-in-r-choropleth-maps/>
>>>>> >
>>>>>
>>>>>
>>>>> Jean
>>>>>
>>>>>
>>>>> On Tue, Oct 29, 2013 at 12:02 PM, <paladini at trustindata.de> wrote:
>>>>>
>>>>>  Hello,
>>>>>
>>>>>> I would like to draw a map of Europe. Each country should be colored
>>>>>> depending on how it scores in an index called GPIndex.
>>>>>> Say a dark red for real bad countries a light red for those which are
>>>>>> not
>>>>>> so bad, light blue for the fairly good ones and so on up to the really
>>>>>> good
>>>>>> ones in a dark blue.
>>>>>> I never worked with geographic maps before so I tried library maps but
>>>>>> I
>>>>>> didn't get far,- especially because all examples I found only seem to
>>>>>> work
>>>>>> for the United states. So I'm a bit lost.
>>>>>> I would be nice if somebody could help me.
>>>>>>
>>>>>> Thanking you in anticipation!
>>>>>>
>>>>>> Best regards
>>>>>>
>>>>>> Claudia
>>>>>>
>>>>>> ______________________________******________________
>>>>>> R-help at r-project.org mailing list
>>>>>> https://stat.ethz.ch/mailman/******listinfo/r-help<https://stat.ethz.ch/mailman/****listinfo/r-help>
>>>>>> <https://**stat.ethz.ch/mailman/****listinfo/r-help<https://stat.ethz.ch/mailman/**listinfo/r-help>
>>>>>> >
>>>>>> <https://stat.**ethz.ch/**mailman/listinfo/r-**help<http://ethz.ch/mailman/listinfo/r-**help>
>>>>>> <http**s://stat.ethz.ch/mailman/**listinfo/r-help<https://stat.ethz.ch/mailman/listinfo/r-help>
>>>>>> >
>>>>>>
>>>>>> >
>>>>>> PLEASE do read the posting guide http://www.R-project.org/**
>>>>>> posting-guide.html <http://www.R-project.org/****posting-guide.html<http://www.R-project.org/**posting-guide.html>
>>>>>> <http://www.**R-project.org/posting-guide.**html<http://www.R-project.org/posting-guide.html>
>>>>>> >
>>>>>>
>>>>>> >
>>>>>>
>>>>>> and provide commented, minimal, self-contained, reproducible code.
>>>>>>
>>>>>>
>>>>>>
>>>>
>>>>
>>
>>
>
>         [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.



-- 

Bert Gunter
Genentech Nonclinical Biostatistics

(650) 467-7374



More information about the R-help mailing list