Dear All,
I need to do something relatively simple: generate a map of Europe and
paint the various states with different colors (only 3-4 are needed)
according to a rule.
I would like to keep it as simple as possible and the script resorting to
the sp package that I found at
http://bit.ly/Oc71ub
is exactly what I am looking for (I need to repeat the exercise with
Europe instead of Switzerland and to have a legend).
However, when I try to run the script (which I also paste at the end of
the email), well, nothing happens and I certainly do not see a map of
Switzerland.
I am not familiar with the sp package, but am I missing the obvious?
Any help is appreciated.
Best Regards
Lorenzo
########################################?
library(sp)
con <- url("http://gadm.org/data/rda/CHE_adm1.RData")
print(load(con))
close(con)
language <- c("german", "german",
"german","german",
"german","german","french", "french",
"german","german","french", "french",
"german", "french","german","german",
"german","german","german", "german",
"german","italian","german","french",
"french","german","german")
gadm$language <- as.factor(language)
col = rainbow(length(levels(gadm$language)))
spplot(gadm, "language", col.regions=col, main="Swiss Language
Regions")
###########################################
The script is working for me. Did you get any error or warning messages? What version of R are you using? What operating system does your computer use? When the command, print(load(con)), runs, what message do you see? ---------------------------------------------- David L Carlson Associate Professor of Anthropology Texas A&M University College Station, TX 77843-4352> -----Original Message----- > From: r-help-bounces at r-project.org [mailto:r-help-bounces at r- > project.org] On Behalf Of Lorenzo Isella > Sent: Friday, August 10, 2012 3:21 PM > To: r-help at stat.math.ethz.ch > Subject: [R] Trouble with Spatial Data Example Script > > Dear All, > I need to do something relatively simple: generate a map of Europe and > paint the various states with different colors (only 3-4 are needed) > according to a rule. > I would like to keep it as simple as possible and the script resorting > to > the sp package that I found at > > http://bit.ly/Oc71ub > > is exactly what I am looking for (I need to repeat the exercise with > Europe instead of Switzerland and to have a legend). > However, when I try to run the script (which I also paste at the end of > the email), well, nothing happens and I certainly do not see a map of > Switzerland. > I am not familiar with the sp package, but am I missing the obvious? > Any help is appreciated. > Best Regards > > Lorenzo > > ########################################? > > library(sp) > con <- url("http://gadm.org/data/rda/CHE_adm1.RData") > print(load(con)) > close(con) > > > > language <- c("german", "german", "german","german", > "german","german","french", "french", > "german","german","french", "french", > "german", "french","german","german", > "german","german","german", "german", > "german","italian","german","french", > "french","german","german") > gadm$language <- as.factor(language) > col = rainbow(length(levels(gadm$language))) > spplot(gadm, "language", col.regions=col, main="Swiss Language > Regions") > > ########################################### > > ______________________________________________ > 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.
As David Carlson mentioned earlier your script is working well for me. sessionInfo() R version 2.15.1 (2012-06-22) Platform: i686-pc-linux-gnu (32-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_CA.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=en_CA.UTF-8 LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=C LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=en_CA.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] sp_0.9-99 loaded via a namespace (and not attached): [1] fortunes_1.5-0 grid_2.15.1 lattice_0.20-6 John Kane Kingston ON Canada> -----Original Message----- > From: lorenzo.isella at gmail.com > Sent: Fri, 10 Aug 2012 22:20:56 +0200 > To: r-help at stat.math.ethz.ch > Subject: [R] Trouble with Spatial Data Example Script > > Dear All, > I need to do something relatively simple: generate a map of Europe and > paint the various states with different colors (only 3-4 are needed) > according to a rule. > I would like to keep it as simple as possible and the script resorting to > the sp package that I found at > > http://bit.ly/Oc71ub > > is exactly what I am looking for (I need to repeat the exercise with > Europe instead of Switzerland and to have a legend). > However, when I try to run the script (which I also paste at the end of > the email), well, nothing happens and I certainly do not see a map of > Switzerland. > I am not familiar with the sp package, but am I missing the obvious? > Any help is appreciated. > Best Regards > > Lorenzo > > ########################################? > > library(sp) > con <- url("http://gadm.org/data/rda/CHE_adm1.RData") > print(load(con)) > close(con) > > > > language <- c("german", "german", "german","german", > "german","german","french", "french", > "german","german","french", "french", > "german", "french","german","german", > "german","german","german", "german", > "german","italian","german","french", > "french","german","german") > gadm$language <- as.factor(language) > col = rainbow(length(levels(gadm$language))) > spplot(gadm, "language", col.regions=col, main="Swiss Language Regions") > > ########################################### > > ______________________________________________ > 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.____________________________________________________________ GET FREE SMILEYS FOR YOUR IM & EMAIL - Learn more at http://www.inbox.com/smileys Works with AIM?, MSN? Messenger, Yahoo!? Messenger, ICQ?, Google Talk? and most webmails
Hello,
The situation is a bit odd.
If I open an R session and I paste there the script, then it does work.
However, if I save it as a text file and then I run it with the
source("myfile.R") command, then I fall back to the case that started
this
thread i.e. I see nothing.
Does anybody have a clue at what is going on?
I paste below the output of sessionInfo(), which is the same no matter
whether I paste the commands in a shell or use a .R file.
sessionInfo()
R version 2.15.1 (2012-06-22)
Platform: x86_64-pc-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
[5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
[7] LC_PAPER=C LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] sp_0.9-99
loaded via a namespace (and not attached):
[1] grid_2.15.1 lattice_0.20-6
Many thanks
Lorenzo
On Sat, 11 Aug 2012 17:18:43 +0200, David L Carlson <dcarlson at tamu.edu>
wrote:
> Can you get other plotting commands to work? E.g.
>
> plot(1:5, 1:5)
>
>> -----Original Message-----
>> From: Lorenzo Isella [mailto:lorenzo.isella at gmail.com]
>> Sent: Saturday, August 11, 2012 3:20 AM
>> To: dcarlson at tamu.edu
>> Subject: Re: [R] Trouble with Spatial Data Example Script
>>
>> Thanks for your reply.
>> I am running R 2.15.1 on debian testing, amd64 architecture.
>>
>>
>>
>>
>> ~$ R
>>
>> R version 2.15.1 (2012-06-22) -- "Roasted Marshmallows"
>> Copyright (C) 2012 The R Foundation for Statistical Computing
>> ISBN 3-900051-07-0
>> Platform: x86_64-pc-linux-gnu (64-bit)
>>
>>
>> and when I run the script, this is what I see
>>
>>
>>
>> source('map-switzerland.R')
>> [1] "gadm"
>>
>> gadm contains the data, but I do not see any map.
>> Any hint about what is going on?
>> Cheers
>>
>> Lorenzo
>>
>>
>>
>>
>>
>> On Fri, 10 Aug 2012 23:43:34 +0200, David L Carlson <dcarlson at
tamu.edu>
>> wrote:
>>
>> > The script is working for me. Did you get any error or warning
>> messages?
>> > What version of R are you using? What operating system does your
>> > computer use? When the command, print(load(con)), runs, what
message
>> do
>> > you see?
>> >
>> > ----------------------------------------------
>> > David L Carlson
>> > Associate Professor of Anthropology
>> > Texas A&M University
>> > College Station, TX 77843-4352
>> >
>> >
>> >> -----Original Message-----
>> >> From: r-help-bounces at r-project.org [mailto:r-help-bounces
at r-
>> >> project.org] On Behalf Of Lorenzo Isella
>> >> Sent: Friday, August 10, 2012 3:21 PM
>> >> To: r-help at stat.math.ethz.ch
>> >> Subject: [R] Trouble with Spatial Data Example Script
>> >>
>> >> Dear All,
>> >> I need to do something relatively simple: generate a map of
Europe
>> and
>> >> paint the various states with different colors (only 3-4 are
needed)
>> >> according to a rule.
>> >> I would like to keep it as simple as possible and the script
>> resorting
>> >> to
>> >> the sp package that I found at
>> >>
>> >> http://bit.ly/Oc71ub
>> >>
>> >> is exactly what I am looking for (I need to repeat the
exercise with
>> >> Europe instead of Switzerland and to have a legend).
>> >> However, when I try to run the script (which I also paste at
the end
>> of
>> >> the email), well, nothing happens and I certainly do not see a
map
>> of
>> >> Switzerland.
>> >> I am not familiar with the sp package, but am I missing the
obvious?
>> >> Any help is appreciated.
>> >> Best Regards
>> >>
>> >> Lorenzo
>> >>
>> >> ########################################?
>> >>
>> >> library(sp)
>> >> con <-
url("http://gadm.org/data/rda/CHE_adm1.RData")
>> >> print(load(con))
>> >> close(con)
>> >>
>> >>
>> >>
>> >> language <- c("german", "german",
"german","german",
>> >> "german","german","french",
"french",
>> >> "german","german","french",
"french",
>> >> "german",
"french","german","german",
>> >> "german","german","german",
"german",
>> >>
"german","italian","german","french",
>> >> "french","german","german")
>> >> gadm$language <- as.factor(language)
>> >> col = rainbow(length(levels(gadm$language)))
>> >> spplot(gadm, "language", col.regions=col,
main="Swiss Language
>> >> Regions")
>> >>
>> >> ###########################################
>> >>
>> >> ______________________________________________
>> >> 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.