R › R help

Trouble with Spatial Data Example Script

Classic

List

Threaded

6 messages Options

Lorenzo Isella

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")

###########################################

______________________________________________
[hidden email] 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.

David Carlson

Re: Trouble with Spatial Data Example Script

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: [hidden email] [mailto:r-help-bounces@r-
> project.org] On Behalf Of Lorenzo Isella
> Sent: Friday, August 10, 2012 3:21 PM
> To: [hidden email]
> 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")
>
> ###########################################
>
> ______________________________________________
> [hidden email] 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.

______________________________________________
[hidden email] 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.

John Kane

Re: Trouble with Spatial Data Example Script

In reply to this post by Lorenzo Isella

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: [hidden email]
> Sent: Fri, 10 Aug 2012 22:20:56 +0200
> To: [hidden email]
> 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")
>
> ###########################################
>
> ______________________________________________
> [hidden email] 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

______________________________________________
[hidden email] 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.

Lorenzo Isella

Re: Trouble with Spatial Data Example Script

In reply to this post by David Carlson

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 <[hidden email]>
wrote:

> Can you get other plotting commands to work? E.g.
>
> plot(1:5, 1:5)
>
>> -----Original Message-----
>> From: Lorenzo Isella [mailto:[hidden email]]
>> Sent: Saturday, August 11, 2012 3:20 AM
>> To: [hidden email]
>> 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 <[hidden email]>
>> 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: [hidden email] [mailto:r-help-bounces@r-
>> >> project.org] On Behalf Of Lorenzo Isella
>> >> Sent: Friday, August 10, 2012 3:21 PM
>> >> To: [hidden email]
>> >> 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")
>> >>
>> >> ###########################################
>> >>
>> >> ______________________________________________
>> >> [hidden email] 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.

Michael Weylandt

Re: Trouble with Spatial Data Example Script

I have no experience with spplot() or maps in general, but might this
be more or less tied to R FAQ 7.22?

Michael

On Sat, Aug 11, 2012 at 12:00 PM, Lorenzo Isella
<[hidden email]> wrote:

> 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 <[hidden email]>
> wrote:
>
>> Can you get other plotting commands to work? E.g.
>>
>> plot(1:5, 1:5)
>>
>>> -----Original Message-----
>>> From: Lorenzo Isella [mailto:[hidden email]]
>>> Sent: Saturday, August 11, 2012 3:20 AM
>>> To: [hidden email]
>>> 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 <[hidden email]>
>>> 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: [hidden email] [mailto:r-help-bounces@r-
>>> >> project.org] On Behalf Of Lorenzo Isella
>>> >> Sent: Friday, August 10, 2012 3:21 PM
>>> >> To: [hidden email]
>>> >> 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")
>>> >>
>>> >> ###########################################
>>> >>
>>> >> ______________________________________________
>>> >> [hidden email] 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.
>
> ______________________________________________
> [hidden email] 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.

David Carlson

Re: Trouble with Spatial Data Example Script

That is it. If you are going to use source() change

spplot(gadm, "language", col.regions=col, main="Swiss Language
Regions")

To

print(spplot(gadm, "language", col.regions=col, main="Swiss Language
Regions"))

-------
David

> -----Original Message-----
> From: R. Michael Weylandt [mailto:[hidden email]]
> Sent: Saturday, August 11, 2012 12:27 PM
> To: Lorenzo Isella
> Cc: [hidden email]; John Kane; [hidden email]
> Subject: Re: [R] Trouble with Spatial Data Example Script
>
> I have no experience with spplot() or maps in general, but might this
> be more or less tied to R FAQ 7.22?
>
> Michael
>
> On Sat, Aug 11, 2012 at 12:00 PM, Lorenzo Isella
> <[hidden email]> wrote:
> > 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
> <[hidden email]>
> > wrote:
> >
> >> Can you get other plotting commands to work? E.g.
> >>
> >> plot(1:5, 1:5)
> >>
> >>> -----Original Message-----
> >>> From: Lorenzo Isella [mailto:[hidden email]]
> >>> Sent: Saturday, August 11, 2012 3:20 AM
> >>> To: [hidden email]
> >>> 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
> <[hidden email]>
> >>> 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: [hidden email] [mailto:r-help-bounces@r-
> >>> >> project.org] On Behalf Of Lorenzo Isella
> >>> >> Sent: Friday, August 10, 2012 3:21 PM
> >>> >> To: [hidden email]
> >>> >> 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")
> >>> >>
> >>> >> ###########################################
> >>> >>
> >>> >> ______________________________________________
> >>> >> [hidden email] 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.
> >
> > ______________________________________________
> > [hidden email] 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.