R › R help

time series line plot: Error in plot.window(...) : invalid 'xlim' value

Classic

List

Threaded

12 messages Options

ytra

time series line plot: Error in plot.window(...) : invalid 'xlim' value

Hello,

This should be pretty simple but I cannot get it right. Please point to the
right code. Thanks.

> last <- read.csv(file.path(dataDir,"plot1.csv"), as.is=T,stringsAsFactors
= FALSE)
> last
date r_wvht
1 8/6/2008 0.9766667
2 8/8/2008 0.7733333
3 8/11/2008 1.4833333
4 8/13/2008 1.5766667
5 8/14/2008 1.3900000
6 8/18/2008 0.7800000
7 8/20/2008 0.8383333
8 8/27/2008 1.7700000
9 8/28/2008 1.2950000
10 8/31/2008 2.4100000
11 9/2/2008 1.3166667
12 9/3/2008 1.3075000
13 9/4/2008 1.3900000
14 9/5/2008 1.6333333
15 9/8/2008 1.2416667
16 9/11/2008 1.3950000
17 10/12/2009 0.8633333
18 10/14/2009 2.7900000
19 10/19/2009 1.0325000
20 10/21/2009 1.9650000
21 10/26/2009 1.7800000
22 10/29/2009 1.5666667
23 10/30/2009 1.0500000
24 11/2/2009 1.4633333
25 11/3/2009 1.2400000
26 11/4/2009 1.0075000
27 11/11/2009 1.6050000
28 11/13/2009 1.8475000
> x<-as.vector(last$date)
> x
[1] "8/6/2008" "8/8/2008" "8/11/2008" "8/13/2008" "8/14/2008"
"8/18/2008" "8/20/2008" "8/27/2008"
[9] "8/28/2008" "8/31/2008" "9/2/2008" "9/3/2008" "9/4/2008"
"9/5/2008" "9/8/2008" "9/11/2008"
[17] "10/12/2009" "10/14/2009" "10/19/2009" "10/21/2009" "10/26/2009"
"10/29/2009" "10/30/2009" "11/2/2009"
[25] "11/3/2009" "11/4/2009" "11/11/2009" "11/13/2009"
> y<-as.vector(last$r_wvht)
> y
[1] 0.9766667 0.7733333 1.4833333 1.5766667 1.3900000 0.7800000 0.8383333
1.7700000 1.2950000 2.4100000
[11] 1.3166667 1.3075000 1.3900000 1.6333333 1.2416667 1.3950000 0.8633333
2.7900000 1.0325000 1.9650000
[21] 1.7800000 1.5666667 1.0500000 1.4633333 1.2400000 1.0075000 1.6050000
1.8475000
> plot(x,y)
Error in plot.window(...) : need finite 'xlim' values
In addition: Warning messages:
1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
2: In min(x) : no non-missing arguments to min; returning Inf
3: In max(x) : no non-missing arguments to max; returning -Inf
> plot(x,y,xlim=c("6/8/2008","11/13/2009"))
Error in plot.window(...) : invalid 'xlim' value
Y

[[alternative HTML version deleted]]

______________________________________________
[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.

arun kirshna

Re: time series line plot: Error in plot.window(...) : invalid 'xlim' value

This post has NOT been accepted by the mailing list yet.

Hi,

Try converting x to Date format

x1<-as.Date(x,format="%m/%d/%Y")
y<-as.vector(dat1$date)
plot(x1,y)

A.K.

Rui Barradas

Re: time series line plot: Error in plot.window(...) : invalid 'xlim' value

In reply to this post by ytra

Hello,

You could use dput(), it's not your first post...

last <- structure(list(date = c("8/6/2008", "8/8/2008", "8/11/2008",
"8/13/2008", "8/14/2008", "8/18/2008", "8/20/2008", "8/27/2008",
"8/28/2008", "8/31/2008", "9/2/2008", "9/3/2008", "9/4/2008",
"9/5/2008", "9/8/2008", "9/11/2008", "10/12/2009", "10/14/2009",
"10/19/2009", "10/21/2009", "10/26/2009", "10/29/2009", "10/30/2009",
"11/2/2009", "11/3/2009", "11/4/2009", "11/11/2009", "11/13/2009"
), r_wvht = c(0.9766667, 0.7733333, 1.4833333, 1.5766667, 1.39,
0.78, 0.8383333, 1.77, 1.295, 2.41, 1.3166667, 1.3075, 1.39,
1.6333333, 1.2416667, 1.395, 0.8633333, 2.79, 1.0325, 1.965,
1.78, 1.5666667, 1.05, 1.4633333, 1.24, 1.0075, 1.605, 1.8475
)), .Names = c("date", "r_wvht"), class = "data.frame", row.names = c("1",
"2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13",
"14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24",
"25", "26", "27", "28"))

last$date <- as.Date(last$date, format="%m/%d/%Y")

plot(r_wvht ~ date, data = last)

(You were forgetting that last$date was a character vector, not class Date.)

Hope this helps,

Rui Barradas

Em 31-07-2012 18:02, Yolande Tra escreveu:

> Hello,
>
> This should be pretty simple but I cannot get it right. Please point to the
> right code. Thanks.
>
>> last <- read.csv(file.path(dataDir,"plot1.csv"), as.is=T,stringsAsFactors
> = FALSE)
>> last
> date r_wvht
> 1 8/6/2008 0.9766667
> 2 8/8/2008 0.7733333
> 3 8/11/2008 1.4833333
> 4 8/13/2008 1.5766667
> 5 8/14/2008 1.3900000
> 6 8/18/2008 0.7800000
> 7 8/20/2008 0.8383333
> 8 8/27/2008 1.7700000
> 9 8/28/2008 1.2950000
> 10 8/31/2008 2.4100000
> 11 9/2/2008 1.3166667
> 12 9/3/2008 1.3075000
> 13 9/4/2008 1.3900000
> 14 9/5/2008 1.6333333
> 15 9/8/2008 1.2416667
> 16 9/11/2008 1.3950000
> 17 10/12/2009 0.8633333
> 18 10/14/2009 2.7900000
> 19 10/19/2009 1.0325000
> 20 10/21/2009 1.9650000
> 21 10/26/2009 1.7800000
> 22 10/29/2009 1.5666667
> 23 10/30/2009 1.0500000
> 24 11/2/2009 1.4633333
> 25 11/3/2009 1.2400000
> 26 11/4/2009 1.0075000
> 27 11/11/2009 1.6050000
> 28 11/13/2009 1.8475000
>> x<-as.vector(last$date)
>> x
> [1] "8/6/2008" "8/8/2008" "8/11/2008" "8/13/2008" "8/14/2008"
> "8/18/2008" "8/20/2008" "8/27/2008"
> [9] "8/28/2008" "8/31/2008" "9/2/2008" "9/3/2008" "9/4/2008"
> "9/5/2008" "9/8/2008" "9/11/2008"
> [17] "10/12/2009" "10/14/2009" "10/19/2009" "10/21/2009" "10/26/2009"
> "10/29/2009" "10/30/2009" "11/2/2009"
> [25] "11/3/2009" "11/4/2009" "11/11/2009" "11/13/2009"
>> y<-as.vector(last$r_wvht)
>> y
> [1] 0.9766667 0.7733333 1.4833333 1.5766667 1.3900000 0.7800000 0.8383333
> 1.7700000 1.2950000 2.4100000
> [11] 1.3166667 1.3075000 1.3900000 1.6333333 1.2416667 1.3950000 0.8633333
> 2.7900000 1.0325000 1.9650000
> [21] 1.7800000 1.5666667 1.0500000 1.4633333 1.2400000 1.0075000 1.6050000
> 1.8475000
>> plot(x,y)
> Error in plot.window(...) : need finite 'xlim' values
> In addition: Warning messages:
> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
> 2: In min(x) : no non-missing arguments to min; returning Inf
> 3: In max(x) : no non-missing arguments to max; returning -Inf
>> plot(x,y,xlim=c("6/8/2008","11/13/2009"))
> Error in plot.window(...) : invalid 'xlim' value
> Y
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> [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.

Michael Weylandt

Re: time series line plot: Error in plot.window(...) : invalid 'xlim' value

In reply to this post by ytra

On Tue, Jul 31, 2012 at 12:02 PM, Yolande Tra <[hidden email]> wrote:

Here's your problem -- "x" is all strings so you can't make an x-axis
out of them later. I'll show you how to diagnose that below. You
probably want as.Date() here -- I rarely (never?) have occasion to use
as.vector().

>> x
> [1] "8/6/2008" "8/8/2008" "8/11/2008" "8/13/2008" "8/14/2008"
> "8/18/2008" "8/20/2008" "8/27/2008"
> [9] "8/28/2008" "8/31/2008" "9/2/2008" "9/3/2008" "9/4/2008"
> "9/5/2008" "9/8/2008" "9/11/2008"
> [17] "10/12/2009" "10/14/2009" "10/19/2009" "10/21/2009" "10/26/2009"
> "10/29/2009" "10/30/2009" "11/2/2009"
> [25] "11/3/2009" "11/4/2009" "11/11/2009" "11/13/2009"
>> y<-as.vector(last$r_wvht)
>> y
> [1] 0.9766667 0.7733333 1.4833333 1.5766667 1.3900000 0.7800000 0.8383333
> 1.7700000 1.2950000 2.4100000
> [11] 1.3166667 1.3075000 1.3900000 1.6333333 1.2416667 1.3950000 0.8633333
> 2.7900000 1.0325000 1.9650000
> [21] 1.7800000 1.5666667 1.0500000 1.4633333 1.2400000 1.0075000 1.6050000
> 1.8475000
>> plot(x,y)
> Error in plot.window(...) : need finite 'xlim' values
> In addition: Warning messages:
> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
> 2: In min(x) : no non-missing arguments to min; returning Inf
> 3: In max(x) : no non-missing arguments to max; returning -Inf
>> plot(x,y,xlim=c("6/8/2008","11/13/2009"))
> Error in plot.window(...) : invalid 'xlim' value

Taking these in order --

1) we see something had to be coerced: I know it was the strings above
and all of them created NA's because they weren't numeric literals.
2) Since all our arguments were NA to min, we have no information
about the min() and return Inf.
3) Ibid, returning -Inf.

giving in turn the error:

x limits from -Inf to Inf don't make a well defined graph.

Big picture, you should probably use a zoo object here instead of a
data frame and awkward column conversion:

library(zoo)
x <- read.zoo( file_path) # Creates a "zoo" object x

plot(x) # Will use the index of x automatically as the x axis and the
values as the y because x is a zoo (read time series) object.

and it should all happen magically.

Michael

> Y
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> [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.

Rui Barradas

Re: time series line plot: Error in plot.window(...) : invalid 'xlim' value

In reply to this post by Rui Barradas

Hello,

Sorry, I forgot the "time series" part of your question. You could use
instead one of

# 1. type = "l" gives a line plot
plot(r_wvht ~ date, data = last, type="l")

# 2. use time series object plot
library(zoo)
z <- zoo(last$r_wvht, order.by=last$date)
plot(z)

Rui Barradas

Em 31-07-2012 18:27, Rui Barradas escreveu:

> Hello,
>
> You could use dput(), it's not your first post...
>
>
> last <- structure(list(date = c("8/6/2008", "8/8/2008", "8/11/2008",
> "8/13/2008", "8/14/2008", "8/18/2008", "8/20/2008", "8/27/2008",
> "8/28/2008", "8/31/2008", "9/2/2008", "9/3/2008", "9/4/2008",
> "9/5/2008", "9/8/2008", "9/11/2008", "10/12/2009", "10/14/2009",
> "10/19/2009", "10/21/2009", "10/26/2009", "10/29/2009", "10/30/2009",
> "11/2/2009", "11/3/2009", "11/4/2009", "11/11/2009", "11/13/2009"
> ), r_wvht = c(0.9766667, 0.7733333, 1.4833333, 1.5766667, 1.39,
> 0.78, 0.8383333, 1.77, 1.295, 2.41, 1.3166667, 1.3075, 1.39,
> 1.6333333, 1.2416667, 1.395, 0.8633333, 2.79, 1.0325, 1.965,
> 1.78, 1.5666667, 1.05, 1.4633333, 1.24, 1.0075, 1.605, 1.8475
> )), .Names = c("date", "r_wvht"), class = "data.frame", row.names =
> c("1",
> "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13",
> "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24",
> "25", "26", "27", "28"))
>
> last$date <- as.Date(last$date, format="%m/%d/%Y")
>
> plot(r_wvht ~ date, data = last)
>
> (You were forgetting that last$date was a character vector, not class
> Date.)
>
> Hope this helps,
>
> Rui Barradas
>
> Em 31-07-2012 18:02, Yolande Tra escreveu:
>> Hello,
>>
>> This should be pretty simple but I cannot get it right. Please point
>> to the
>> right code. Thanks.
>>
>>> last <- read.csv(file.path(dataDir,"plot1.csv"),
>>> as.is=T,stringsAsFactors
>> = FALSE)
>>> last
>> date r_wvht
>> 1 8/6/2008 0.9766667
>> 2 8/8/2008 0.7733333
>> 3 8/11/2008 1.4833333
>> 4 8/13/2008 1.5766667
>> 5 8/14/2008 1.3900000
>> 6 8/18/2008 0.7800000
>> 7 8/20/2008 0.8383333
>> 8 8/27/2008 1.7700000
>> 9 8/28/2008 1.2950000
>> 10 8/31/2008 2.4100000
>> 11 9/2/2008 1.3166667
>> 12 9/3/2008 1.3075000
>> 13 9/4/2008 1.3900000
>> 14 9/5/2008 1.6333333
>> 15 9/8/2008 1.2416667
>> 16 9/11/2008 1.3950000
>> 17 10/12/2009 0.8633333
>> 18 10/14/2009 2.7900000
>> 19 10/19/2009 1.0325000
>> 20 10/21/2009 1.9650000
>> 21 10/26/2009 1.7800000
>> 22 10/29/2009 1.5666667
>> 23 10/30/2009 1.0500000
>> 24 11/2/2009 1.4633333
>> 25 11/3/2009 1.2400000
>> 26 11/4/2009 1.0075000
>> 27 11/11/2009 1.6050000
>> 28 11/13/2009 1.8475000
>>> x<-as.vector(last$date)
>>> x
>> [1] "8/6/2008" "8/8/2008" "8/11/2008" "8/13/2008" "8/14/2008"
>> "8/18/2008" "8/20/2008" "8/27/2008"
>> [9] "8/28/2008" "8/31/2008" "9/2/2008" "9/3/2008" "9/4/2008"
>> "9/5/2008" "9/8/2008" "9/11/2008"
>> [17] "10/12/2009" "10/14/2009" "10/19/2009" "10/21/2009" "10/26/2009"
>> "10/29/2009" "10/30/2009" "11/2/2009"
>> [25] "11/3/2009" "11/4/2009" "11/11/2009" "11/13/2009"
>>> y<-as.vector(last$r_wvht)
>>> y
>> [1] 0.9766667 0.7733333 1.4833333 1.5766667 1.3900000 0.7800000
>> 0.8383333
>> 1.7700000 1.2950000 2.4100000
>> [11] 1.3166667 1.3075000 1.3900000 1.6333333 1.2416667 1.3950000
>> 0.8633333
>> 2.7900000 1.0325000 1.9650000
>> [21] 1.7800000 1.5666667 1.0500000 1.4633333 1.2400000 1.0075000
>> 1.6050000
>> 1.8475000
>>> plot(x,y)
>> Error in plot.window(...) : need finite 'xlim' values
>> In addition: Warning messages:
>> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
>> 2: In min(x) : no non-missing arguments to min; returning Inf
>> 3: In max(x) : no non-missing arguments to max; returning -Inf
>>> plot(x,y,xlim=c("6/8/2008","11/13/2009"))
>> Error in plot.window(...) : invalid 'xlim' value
>> Y
>>
>> [[alternative HTML version deleted]]
>>
>> ______________________________________________
>> [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.

ytra

Re: time series line plot: Error in plot.window(...) : invalid 'xlim' value

Thank you everyone for the attempt to solve the problem
It is an irregular series and insert NAs when a date is missing
> library(zoo)
> z <- zoo(last$r_wvht, order.by=last$date)
> plot(z)
Error in plot.window(...) : need finite 'xlim' values
In addition: Warning messages:
1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
2: In min(x) : no non-missing arguments to min; returning Inf
3: In max(x) : no non-missing arguments to max; returning -Inf

I have also tried

> z <- read.zoo("plot1.csv", header = TRUE, sep = ",", format = "%m/%d/%Y")
> plot(z)
The plot does not look good. It does not display the individual dates.

Y
Y
On Tue, Jul 31, 2012 at 1:33 PM, Rui Barradas <[hidden email]> wrote:

> Hello,
>
> Sorry, I forgot the "time series" part of your question. You could use
> instead one of
>
> # 1. type = "l" gives a line plot
> plot(r_wvht ~ date, data = last, type="l")
>
> # 2. use time series object plot
> library(zoo)
> z <- zoo(last$r_wvht, order.by=last$date)
> plot(z)
>
> Rui Barradas
>
> Em 31-07-2012 18:27, Rui Barradas escreveu:
>
> Hello,
>>
>> You could use dput(), it's not your first post...
>>
>>
>> last <- structure(list(date = c("8/6/2008", "8/8/2008", "8/11/2008",
>> "8/13/2008", "8/14/2008", "8/18/2008", "8/20/2008", "8/27/2008",
>> "8/28/2008", "8/31/2008", "9/2/2008", "9/3/2008", "9/4/2008",
>> "9/5/2008", "9/8/2008", "9/11/2008", "10/12/2009", "10/14/2009",
>> "10/19/2009", "10/21/2009", "10/26/2009", "10/29/2009", "10/30/2009",
>> "11/2/2009", "11/3/2009", "11/4/2009", "11/11/2009", "11/13/2009"
>> ), r_wvht = c(0.9766667, 0.7733333, 1.4833333, 1.5766667, 1.39,
>> 0.78, 0.8383333, 1.77, 1.295, 2.41, 1.3166667, 1.3075, 1.39,
>> 1.6333333, 1.2416667, 1.395, 0.8633333, 2.79, 1.0325, 1.965,
>> 1.78, 1.5666667, 1.05, 1.4633333, 1.24, 1.0075, 1.605, 1.8475
>> )), .Names = c("date", "r_wvht"), class = "data.frame", row.names = c("1",
>> "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13",
>> "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24",
>> "25", "26", "27", "28"))
>>
>> last$date <- as.Date(last$date, format="%m/%d/%Y")
>>
>> plot(r_wvht ~ date, data = last)
>>
>> (You were forgetting that last$date was a character vector, not class
>> Date.)
>>
>> Hope this helps,
>>
>> Rui Barradas
>>
>> Em 31-07-2012 18:02, Yolande Tra escreveu:
>>
>>> Hello,
>>>
>>> This should be pretty simple but I cannot get it right. Please point to
>>> the
>>> right code. Thanks.
>>>
>>> last <- read.csv(file.path(dataDir,"**plot1.csv"), as.is
>>>> =T,stringsAsFactors
>>>>
>>> = FALSE)
>>>
>>>> last
>>>>
>>> date r_wvht
>>> 1 8/6/2008 0.9766667
>>> 2 8/8/2008 0.7733333
>>> 3 8/11/2008 1.4833333
>>> 4 8/13/2008 1.5766667
>>> 5 8/14/2008 1.3900000
>>> 6 8/18/2008 0.7800000
>>> 7 8/20/2008 0.8383333
>>> 8 8/27/2008 1.7700000
>>> 9 8/28/2008 1.2950000
>>> 10 8/31/2008 2.4100000
>>> 11 9/2/2008 1.3166667
>>> 12 9/3/2008 1.3075000
>>> 13 9/4/2008 1.3900000
>>> 14 9/5/2008 1.6333333
>>> 15 9/8/2008 1.2416667
>>> 16 9/11/2008 1.3950000
>>> 17 10/12/2009 0.8633333
>>> 18 10/14/2009 2.7900000
>>> 19 10/19/2009 1.0325000
>>> 20 10/21/2009 1.9650000
>>> 21 10/26/2009 1.7800000
>>> 22 10/29/2009 1.5666667
>>> 23 10/30/2009 1.0500000
>>> 24 11/2/2009 1.4633333
>>> 25 11/3/2009 1.2400000
>>> 26 11/4/2009 1.0075000
>>> 27 11/11/2009 1.6050000
>>> 28 11/13/2009 1.8475000
>>>
>>>> x<-as.vector(last$date)
>>>> x
>>>>
>>> [1] "8/6/2008" "8/8/2008" "8/11/2008" "8/13/2008" "8/14/2008"
>>> "8/18/2008" "8/20/2008" "8/27/2008"
>>> [9] "8/28/2008" "8/31/2008" "9/2/2008" "9/3/2008" "9/4/2008"
>>> "9/5/2008" "9/8/2008" "9/11/2008"
>>> [17] "10/12/2009" "10/14/2009" "10/19/2009" "10/21/2009" "10/26/2009"
>>> "10/29/2009" "10/30/2009" "11/2/2009"
>>> [25] "11/3/2009" "11/4/2009" "11/11/2009" "11/13/2009"
>>>
>>>> y<-as.vector(last$r_wvht)
>>>> y
>>>>
>>> [1] 0.9766667 0.7733333 1.4833333 1.5766667 1.3900000 0.7800000
>>> 0.8383333
>>> 1.7700000 1.2950000 2.4100000
>>> [11] 1.3166667 1.3075000 1.3900000 1.6333333 1.2416667 1.3950000
>>> 0.8633333
>>> 2.7900000 1.0325000 1.9650000
>>> [21] 1.7800000 1.5666667 1.0500000 1.4633333 1.2400000 1.0075000
>>> 1.6050000
>>> 1.8475000
>>>
>>>> plot(x,y)
>>>>
>>> Error in plot.window(...) : need finite 'xlim' values
>>> In addition: Warning messages:
>>> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
>>> 2: In min(x) : no non-missing arguments to min; returning Inf
>>> 3: In max(x) : no non-missing arguments to max; returning -Inf
>>>
>>>> plot(x,y,xlim=c("6/8/2008","**11/13/2009"))
>>>>
>>> Error in plot.window(...) : invalid 'xlim' value
>>> Y
>>>
>>> [[alternative HTML version deleted]]
>>>
>>> ______________________________**________________
>>> [hidden email] mailing list
>>> https://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>
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>
>> ______________________________**________________
>> [hidden email] mailing list
>> https://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>
>> and provide commented, minimal, self-contained, reproducible code.
>>
>
>

[[alternative HTML version deleted]]

______________________________________________
[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.

Rui Barradas

Re: time series line plot: Error in plot.window(...) : invalid 'xlim' value

Hello,

It only gives that error if you don't

last$date <- as.Date(last$date, format="%m/%d/%Y")

You must have dates, not character values.
Try it, then make a zoo object, then plot it.

Rui Barradas

Em 31-07-2012 18:54, Yolande Tra escreveu:

> Thank you everyone for the attempt to solve the problem
> It is an irregular series and insert NAs when a date is missing
>> library(zoo)
>> z <- zoo(last$r_wvht, order.by=last$date)
>> plot(z)
> Error in plot.window(...) : need finite 'xlim' values
> In addition: Warning messages:
> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
> 2: In min(x) : no non-missing arguments to min; returning Inf
> 3: In max(x) : no non-missing arguments to max; returning -Inf
>
> I have also tried
>
>> z <- read.zoo("plot1.csv", header = TRUE, sep = ",", format = "%m/%d/%Y")
>> plot(z)
> The plot does not look good. It does not display the individual dates.
>
> Y
> Y
> On Tue, Jul 31, 2012 at 1:33 PM, Rui Barradas <[hidden email]> wrote:
>
>> Hello,
>>
>> Sorry, I forgot the "time series" part of your question. You could use
>> instead one of
>>
>> # 1. type = "l" gives a line plot
>> plot(r_wvht ~ date, data = last, type="l")
>>
>> # 2. use time series object plot
>> library(zoo)
>> z <- zoo(last$r_wvht, order.by=last$date)
>> plot(z)
>>
>> Rui Barradas
>>
>> Em 31-07-2012 18:27, Rui Barradas escreveu:
>>
>> Hello,
>>> You could use dput(), it's not your first post...
>>>
>>>
>>> last <- structure(list(date = c("8/6/2008", "8/8/2008", "8/11/2008",
>>> "8/13/2008", "8/14/2008", "8/18/2008", "8/20/2008", "8/27/2008",
>>> "8/28/2008", "8/31/2008", "9/2/2008", "9/3/2008", "9/4/2008",
>>> "9/5/2008", "9/8/2008", "9/11/2008", "10/12/2009", "10/14/2009",
>>> "10/19/2009", "10/21/2009", "10/26/2009", "10/29/2009", "10/30/2009",
>>> "11/2/2009", "11/3/2009", "11/4/2009", "11/11/2009", "11/13/2009"
>>> ), r_wvht = c(0.9766667, 0.7733333, 1.4833333, 1.5766667, 1.39,
>>> 0.78, 0.8383333, 1.77, 1.295, 2.41, 1.3166667, 1.3075, 1.39,
>>> 1.6333333, 1.2416667, 1.395, 0.8633333, 2.79, 1.0325, 1.965,
>>> 1.78, 1.5666667, 1.05, 1.4633333, 1.24, 1.0075, 1.605, 1.8475
>>> )), .Names = c("date", "r_wvht"), class = "data.frame", row.names = c("1",
>>> "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13",
>>> "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24",
>>> "25", "26", "27", "28"))
>>>
>>> last$date <- as.Date(last$date, format="%m/%d/%Y")
>>>
>>> plot(r_wvht ~ date, data = last)
>>>
>>> (You were forgetting that last$date was a character vector, not class
>>> Date.)
>>>
>>> Hope this helps,
>>>
>>> Rui Barradas
>>>
>>> Em 31-07-2012 18:02, Yolande Tra escreveu:
>>>
>>>> Hello,
>>>>
>>>> This should be pretty simple but I cannot get it right. Please point to
>>>> the
>>>> right code. Thanks.
>>>>
>>>> last <- read.csv(file.path(dataDir,"**plot1.csv"), as.is
>>>>> =T,stringsAsFactors
>>>>>
>>>> = FALSE)
>>>>
>>>>> last
>>>>>
>>>> date r_wvht
>>>> 1 8/6/2008 0.9766667
>>>> 2 8/8/2008 0.7733333
>>>> 3 8/11/2008 1.4833333
>>>> 4 8/13/2008 1.5766667
>>>> 5 8/14/2008 1.3900000
>>>> 6 8/18/2008 0.7800000
>>>> 7 8/20/2008 0.8383333
>>>> 8 8/27/2008 1.7700000
>>>> 9 8/28/2008 1.2950000
>>>> 10 8/31/2008 2.4100000
>>>> 11 9/2/2008 1.3166667
>>>> 12 9/3/2008 1.3075000
>>>> 13 9/4/2008 1.3900000
>>>> 14 9/5/2008 1.6333333
>>>> 15 9/8/2008 1.2416667
>>>> 16 9/11/2008 1.3950000
>>>> 17 10/12/2009 0.8633333
>>>> 18 10/14/2009 2.7900000
>>>> 19 10/19/2009 1.0325000
>>>> 20 10/21/2009 1.9650000
>>>> 21 10/26/2009 1.7800000
>>>> 22 10/29/2009 1.5666667
>>>> 23 10/30/2009 1.0500000
>>>> 24 11/2/2009 1.4633333
>>>> 25 11/3/2009 1.2400000
>>>> 26 11/4/2009 1.0075000
>>>> 27 11/11/2009 1.6050000
>>>> 28 11/13/2009 1.8475000
>>>>
>>>>> x<-as.vector(last$date)
>>>>> x
>>>>>
>>>> [1] "8/6/2008" "8/8/2008" "8/11/2008" "8/13/2008" "8/14/2008"
>>>> "8/18/2008" "8/20/2008" "8/27/2008"
>>>> [9] "8/28/2008" "8/31/2008" "9/2/2008" "9/3/2008" "9/4/2008"
>>>> "9/5/2008" "9/8/2008" "9/11/2008"
>>>> [17] "10/12/2009" "10/14/2009" "10/19/2009" "10/21/2009" "10/26/2009"
>>>> "10/29/2009" "10/30/2009" "11/2/2009"
>>>> [25] "11/3/2009" "11/4/2009" "11/11/2009" "11/13/2009"
>>>>
>>>>> y<-as.vector(last$r_wvht)
>>>>> y
>>>>>
>>>> [1] 0.9766667 0.7733333 1.4833333 1.5766667 1.3900000 0.7800000
>>>> 0.8383333
>>>> 1.7700000 1.2950000 2.4100000
>>>> [11] 1.3166667 1.3075000 1.3900000 1.6333333 1.2416667 1.3950000
>>>> 0.8633333
>>>> 2.7900000 1.0325000 1.9650000
>>>> [21] 1.7800000 1.5666667 1.0500000 1.4633333 1.2400000 1.0075000
>>>> 1.6050000
>>>> 1.8475000
>>>>
>>>>> plot(x,y)
>>>>>
>>>> Error in plot.window(...) : need finite 'xlim' values
>>>> In addition: Warning messages:
>>>> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
>>>> 2: In min(x) : no non-missing arguments to min; returning Inf
>>>> 3: In max(x) : no non-missing arguments to max; returning -Inf
>>>>
>>>>> plot(x,y,xlim=c("6/8/2008","**11/13/2009"))
>>>>>
>>>> Error in plot.window(...) : invalid 'xlim' value
>>>> Y
>>>>
>>>> [[alternative HTML version deleted]]
>>>>
>>>> ______________________________**________________
>>>> [hidden email] mailing list
>>>> https://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>
>>>> and provide commented, minimal, self-contained, reproducible code.
>>>>
>>> ______________________________**________________
>>> [hidden email] mailing list
>>> https://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>
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>

ytra

Re: time series line plot: Error in plot.window(...) : invalid 'xlim' value

Thanks.
This is an irregular time series. The line plot does not look good because
of the gap between 9/11/2008 and 10/12/2009. I think two plots would be
better.
How would you include the date on the x-axis. Right now it only gives one
tick mark, 2009.

Yolande

On Tue, Jul 31, 2012 at 2:01 PM, Rui Barradas <[hidden email]> wrote:

> Hello,
>
> It only gives that error if you don't
>
>
> last$date <- as.Date(last$date, format="%m/%d/%Y")
>
> You must have dates, not character values.
> Try it, then make a zoo object, then plot it.
>
> Rui Barradas
>
> Em 31-07-2012 18:54, Yolande Tra escreveu:
>
> Thank you everyone for the attempt to solve the problem
>> It is an irregular series and insert NAs when a date is missing
>>
>>> library(zoo)
>>> z <- zoo(last$r_wvht, order.by=last$date)
>>> plot(z)
>>>
>> Error in plot.window(...) : need finite 'xlim' values
>> In addition: Warning messages:
>> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
>> 2: In min(x) : no non-missing arguments to min; returning Inf
>> 3: In max(x) : no non-missing arguments to max; returning -Inf
>>
>> I have also tried
>>
>> z <- read.zoo("plot1.csv", header = TRUE, sep = ",", format = "%m/%d/%Y")
>>> plot(z)
>>>
>> The plot does not look good. It does not display the individual dates.
>>
>> Y
>> Y
>> On Tue, Jul 31, 2012 at 1:33 PM, Rui Barradas <[hidden email]>
>> wrote:
>>
>> Hello,
>>>
>>> Sorry, I forgot the "time series" part of your question. You could use
>>> instead one of
>>>
>>> # 1. type = "l" gives a line plot
>>> plot(r_wvht ~ date, data = last, type="l")
>>>
>>> # 2. use time series object plot
>>> library(zoo)
>>> z <- zoo(last$r_wvht, order.by=last$date)
>>> plot(z)
>>>
>>> Rui Barradas
>>>
>>> Em 31-07-2012 18:27, Rui Barradas escreveu:
>>>
>>> Hello,
>>>
>>>> You could use dput(), it's not your first post...
>>>>
>>>>
>>>> last <- structure(list(date = c("8/6/2008", "8/8/2008", "8/11/2008",
>>>> "8/13/2008", "8/14/2008", "8/18/2008", "8/20/2008", "8/27/2008",
>>>> "8/28/2008", "8/31/2008", "9/2/2008", "9/3/2008", "9/4/2008",
>>>> "9/5/2008", "9/8/2008", "9/11/2008", "10/12/2009", "10/14/2009",
>>>> "10/19/2009", "10/21/2009", "10/26/2009", "10/29/2009", "10/30/2009",
>>>> "11/2/2009", "11/3/2009", "11/4/2009", "11/11/2009", "11/13/2009"
>>>> ), r_wvht = c(0.9766667, 0.7733333, 1.4833333, 1.5766667, 1.39,
>>>> 0.78, 0.8383333, 1.77, 1.295, 2.41, 1.3166667, 1.3075, 1.39,
>>>> 1.6333333, 1.2416667, 1.395, 0.8633333, 2.79, 1.0325, 1.965,
>>>> 1.78, 1.5666667, 1.05, 1.4633333, 1.24, 1.0075, 1.605, 1.8475
>>>> )), .Names = c("date", "r_wvht"), class = "data.frame", row.names =
>>>> c("1",
>>>> "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13",
>>>> "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24",
>>>> "25", "26", "27", "28"))
>>>>
>>>> last$date <- as.Date(last$date, format="%m/%d/%Y")
>>>>
>>>> plot(r_wvht ~ date, data = last)
>>>>
>>>> (You were forgetting that last$date was a character vector, not class
>>>> Date.)
>>>>
>>>> Hope this helps,
>>>>
>>>> Rui Barradas
>>>>
>>>> Em 31-07-2012 18:02, Yolande Tra escreveu:
>>>>
>>>> Hello,
>>>>>
>>>>> This should be pretty simple but I cannot get it right. Please point to
>>>>> the
>>>>> right code. Thanks.
>>>>>
>>>>> last <- read.csv(file.path(dataDir,"****plot1.csv"), as.is
>>>>>
>>>>> =T,stringsAsFactors
>>>>>>
>>>>>> = FALSE)
>>>>>
>>>>> last
>>>>>>
>>>>>> date r_wvht
>>>>> 1 8/6/2008 0.9766667
>>>>> 2 8/8/2008 0.7733333
>>>>> 3 8/11/2008 1.4833333
>>>>> 4 8/13/2008 1.5766667
>>>>> 5 8/14/2008 1.3900000
>>>>> 6 8/18/2008 0.7800000
>>>>> 7 8/20/2008 0.8383333
>>>>> 8 8/27/2008 1.7700000
>>>>> 9 8/28/2008 1.2950000
>>>>> 10 8/31/2008 2.4100000
>>>>> 11 9/2/2008 1.3166667
>>>>> 12 9/3/2008 1.3075000
>>>>> 13 9/4/2008 1.3900000
>>>>> 14 9/5/2008 1.6333333
>>>>> 15 9/8/2008 1.2416667
>>>>> 16 9/11/2008 1.3950000
>>>>> 17 10/12/2009 0.8633333
>>>>> 18 10/14/2009 2.7900000
>>>>> 19 10/19/2009 1.0325000
>>>>> 20 10/21/2009 1.9650000
>>>>> 21 10/26/2009 1.7800000
>>>>> 22 10/29/2009 1.5666667
>>>>> 23 10/30/2009 1.0500000
>>>>> 24 11/2/2009 1.4633333
>>>>> 25 11/3/2009 1.2400000
>>>>> 26 11/4/2009 1.0075000
>>>>> 27 11/11/2009 1.6050000
>>>>> 28 11/13/2009 1.8475000
>>>>>
>>>>> x<-as.vector(last$date)
>>>>>> x
>>>>>>
>>>>>> [1] "8/6/2008" "8/8/2008" "8/11/2008" "8/13/2008" "8/14/2008"
>>>>> "8/18/2008" "8/20/2008" "8/27/2008"
>>>>> [9] "8/28/2008" "8/31/2008" "9/2/2008" "9/3/2008" "9/4/2008"
>>>>> "9/5/2008" "9/8/2008" "9/11/2008"
>>>>> [17] "10/12/2009" "10/14/2009" "10/19/2009" "10/21/2009" "10/26/2009"
>>>>> "10/29/2009" "10/30/2009" "11/2/2009"
>>>>> [25] "11/3/2009" "11/4/2009" "11/11/2009" "11/13/2009"
>>>>>
>>>>> y<-as.vector(last$r_wvht)
>>>>>> y
>>>>>>
>>>>>> [1] 0.9766667 0.7733333 1.4833333 1.5766667 1.3900000 0.7800000
>>>>> 0.8383333
>>>>> 1.7700000 1.2950000 2.4100000
>>>>> [11] 1.3166667 1.3075000 1.3900000 1.6333333 1.2416667 1.3950000
>>>>> 0.8633333
>>>>> 2.7900000 1.0325000 1.9650000
>>>>> [21] 1.7800000 1.5666667 1.0500000 1.4633333 1.2400000 1.0075000
>>>>> 1.6050000
>>>>> 1.8475000
>>>>>
>>>>> plot(x,y)
>>>>>>
>>>>>> Error in plot.window(...) : need finite 'xlim' values
>>>>> In addition: Warning messages:
>>>>> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
>>>>> 2: In min(x) : no non-missing arguments to min; returning Inf
>>>>> 3: In max(x) : no non-missing arguments to max; returning -Inf
>>>>>
>>>>> plot(x,y,xlim=c("6/8/2008","****11/13/2009"))
>>>>>>
>>>>>> Error in plot.window(...) : invalid 'xlim' value
>>>>> Y
>>>>>
>>>>> [[alternative HTML version deleted]]
>>>>>
>>>>> ______________________________****________________
>>>>> [hidden email] 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>
>>>>> >
>>>>> PLEASE do read the posting guide http://www.R-project.org/**<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.
>>>>>
>>>>> ______________________________****________________
>>>> [hidden email] 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>
>>>> >
>>>> PLEASE do read the posting guide http://www.R-project.org/**<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.
>>>>
>>>>
>>>
>

Rui Barradas

Re: time series line plot: Error in plot.window(...) : invalid 'xlim' value

Hello,

Adapted from the help page for plot.zoo

#library(zoo)
#z <- zoo(last$r_wvht, order.by=last$date)
plot(z, xaxt = "n")
tt <- time(z)
ix <- seq(1, length(tt), length.out=8)
axis(side = 1, at = tt[ix], labels = FALSE)
labs <- format(tt, "%Y-%b-%d")
axis(side = 1, at = tt[ix], labels = labs[ix], tcl = -0.7, cex.axis =
0.7, las=2)

It looks better, but I wouldn't risk "good".

Rui Barradas

Em 31-07-2012 19:10, Yolande Tra escreveu:

> Thanks.
> This is an irregular time series. The line plot does not look good because
> of the gap between 9/11/2008 and 10/12/2009. I think two plots would be
> better.
> How would you include the date on the x-axis. Right now it only gives one
> tick mark, 2009.
>
> Yolande
>
> On Tue, Jul 31, 2012 at 2:01 PM, Rui Barradas <[hidden email]> wrote:
>
>> Hello,
>>
>> It only gives that error if you don't
>>
>>
>> last$date <- as.Date(last$date, format="%m/%d/%Y")
>>
>> You must have dates, not character values.
>> Try it, then make a zoo object, then plot it.
>>
>> Rui Barradas
>>
>> Em 31-07-2012 18:54, Yolande Tra escreveu:
>>
>> Thank you everyone for the attempt to solve the problem
>>> It is an irregular series and insert NAs when a date is missing
>>>
>>>> library(zoo)
>>>> z <- zoo(last$r_wvht, order.by=last$date)
>>>> plot(z)
>>>>
>>> Error in plot.window(...) : need finite 'xlim' values
>>> In addition: Warning messages:
>>> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
>>> 2: In min(x) : no non-missing arguments to min; returning Inf
>>> 3: In max(x) : no non-missing arguments to max; returning -Inf
>>>
>>> I have also tried
>>>
>>> z <- read.zoo("plot1.csv", header = TRUE, sep = ",", format = "%m/%d/%Y")
>>>> plot(z)
>>>>
>>> The plot does not look good. It does not display the individual dates.
>>>
>>> Y
>>> Y
>>> On Tue, Jul 31, 2012 at 1:33 PM, Rui Barradas <[hidden email]>
>>> wrote:
>>>
>>> Hello,
>>>> Sorry, I forgot the "time series" part of your question. You could use
>>>> instead one of
>>>>
>>>> # 1. type = "l" gives a line plot
>>>> plot(r_wvht ~ date, data = last, type="l")
>>>>
>>>> # 2. use time series object plot
>>>> library(zoo)
>>>> z <- zoo(last$r_wvht, order.by=last$date)
>>>> plot(z)
>>>>
>>>> Rui Barradas
>>>>
>>>> Em 31-07-2012 18:27, Rui Barradas escreveu:
>>>>
>>>> Hello,
>>>>
>>>>> You could use dput(), it's not your first post...
>>>>>
>>>>>
>>>>> last <- structure(list(date = c("8/6/2008", "8/8/2008", "8/11/2008",
>>>>> "8/13/2008", "8/14/2008", "8/18/2008", "8/20/2008", "8/27/2008",
>>>>> "8/28/2008", "8/31/2008", "9/2/2008", "9/3/2008", "9/4/2008",
>>>>> "9/5/2008", "9/8/2008", "9/11/2008", "10/12/2009", "10/14/2009",
>>>>> "10/19/2009", "10/21/2009", "10/26/2009", "10/29/2009", "10/30/2009",
>>>>> "11/2/2009", "11/3/2009", "11/4/2009", "11/11/2009", "11/13/2009"
>>>>> ), r_wvht = c(0.9766667, 0.7733333, 1.4833333, 1.5766667, 1.39,
>>>>> 0.78, 0.8383333, 1.77, 1.295, 2.41, 1.3166667, 1.3075, 1.39,
>>>>> 1.6333333, 1.2416667, 1.395, 0.8633333, 2.79, 1.0325, 1.965,
>>>>> 1.78, 1.5666667, 1.05, 1.4633333, 1.24, 1.0075, 1.605, 1.8475
>>>>> )), .Names = c("date", "r_wvht"), class = "data.frame", row.names =
>>>>> c("1",
>>>>> "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13",
>>>>> "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24",
>>>>> "25", "26", "27", "28"))
>>>>>
>>>>> last$date <- as.Date(last$date, format="%m/%d/%Y")
>>>>>
>>>>> plot(r_wvht ~ date, data = last)
>>>>>
>>>>> (You were forgetting that last$date was a character vector, not class
>>>>> Date.)
>>>>>
>>>>> Hope this helps,
>>>>>
>>>>> Rui Barradas
>>>>>
>>>>> Em 31-07-2012 18:02, Yolande Tra escreveu:
>>>>>
>>>>> Hello,
>>>>>> This should be pretty simple but I cannot get it right. Please point to
>>>>>> the
>>>>>> right code. Thanks.
>>>>>>
>>>>>> last <- read.csv(file.path(dataDir,"****plot1.csv"), as.is
>>>>>>
>>>>>> =T,stringsAsFactors
>>>>>>> = FALSE)
>>>>>> last
>>>>>>> date r_wvht
>>>>>> 1 8/6/2008 0.9766667
>>>>>> 2 8/8/2008 0.7733333
>>>>>> 3 8/11/2008 1.4833333
>>>>>> 4 8/13/2008 1.5766667
>>>>>> 5 8/14/2008 1.3900000
>>>>>> 6 8/18/2008 0.7800000
>>>>>> 7 8/20/2008 0.8383333
>>>>>> 8 8/27/2008 1.7700000
>>>>>> 9 8/28/2008 1.2950000
>>>>>> 10 8/31/2008 2.4100000
>>>>>> 11 9/2/2008 1.3166667
>>>>>> 12 9/3/2008 1.3075000
>>>>>> 13 9/4/2008 1.3900000
>>>>>> 14 9/5/2008 1.6333333
>>>>>> 15 9/8/2008 1.2416667
>>>>>> 16 9/11/2008 1.3950000
>>>>>> 17 10/12/2009 0.8633333
>>>>>> 18 10/14/2009 2.7900000
>>>>>> 19 10/19/2009 1.0325000
>>>>>> 20 10/21/2009 1.9650000
>>>>>> 21 10/26/2009 1.7800000
>>>>>> 22 10/29/2009 1.5666667
>>>>>> 23 10/30/2009 1.0500000
>>>>>> 24 11/2/2009 1.4633333
>>>>>> 25 11/3/2009 1.2400000
>>>>>> 26 11/4/2009 1.0075000
>>>>>> 27 11/11/2009 1.6050000
>>>>>> 28 11/13/2009 1.8475000
>>>>>>
>>>>>> x<-as.vector(last$date)
>>>>>>> x
>>>>>>>
>>>>>>> [1] "8/6/2008" "8/8/2008" "8/11/2008" "8/13/2008" "8/14/2008"
>>>>>> "8/18/2008" "8/20/2008" "8/27/2008"
>>>>>> [9] "8/28/2008" "8/31/2008" "9/2/2008" "9/3/2008" "9/4/2008"
>>>>>> "9/5/2008" "9/8/2008" "9/11/2008"
>>>>>> [17] "10/12/2009" "10/14/2009" "10/19/2009" "10/21/2009" "10/26/2009"
>>>>>> "10/29/2009" "10/30/2009" "11/2/2009"
>>>>>> [25] "11/3/2009" "11/4/2009" "11/11/2009" "11/13/2009"
>>>>>>
>>>>>> y<-as.vector(last$r_wvht)
>>>>>>> y
>>>>>>>
>>>>>>> [1] 0.9766667 0.7733333 1.4833333 1.5766667 1.3900000 0.7800000
>>>>>> 0.8383333
>>>>>> 1.7700000 1.2950000 2.4100000
>>>>>> [11] 1.3166667 1.3075000 1.3900000 1.6333333 1.2416667 1.3950000
>>>>>> 0.8633333
>>>>>> 2.7900000 1.0325000 1.9650000
>>>>>> [21] 1.7800000 1.5666667 1.0500000 1.4633333 1.2400000 1.0075000
>>>>>> 1.6050000
>>>>>> 1.8475000
>>>>>>
>>>>>> plot(x,y)
>>>>>>> Error in plot.window(...) : need finite 'xlim' values
>>>>>> In addition: Warning messages:
>>>>>> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
>>>>>> 2: In min(x) : no non-missing arguments to min; returning Inf
>>>>>> 3: In max(x) : no non-missing arguments to max; returning -Inf
>>>>>>
>>>>>> plot(x,y,xlim=c("6/8/2008","****11/13/2009"))
>>>>>>> Error in plot.window(...) : invalid 'xlim' value
>>>>>> Y
>>>>>>
>>>>>> [[alternative HTML version deleted]]
>>>>>>
>>>>>> ______________________________****________________
>>>>>> [hidden email] 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>
>>>>>> PLEASE do read the posting guide http://www.R-project.org/**<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.
>>>>>>
>>>>>> ______________________________****________________
>>>>> [hidden email] 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>
>>>>> PLEASE do read the posting guide http://www.R-project.org/**<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.
>>>>>
>>>>>

ytra

Re: time series line plot: Error in plot.window(...) : invalid 'xlim' value

It does look better indeed. Thanks,

Yolande

On Tue, Jul 31, 2012 at 3:15 PM, Rui Barradas <[hidden email]> wrote:

> Hello,
>
> Adapted from the help page for plot.zoo
>
> #library(zoo)
> #z <- zoo(last$r_wvht, order.by=last$date)
> plot(z, xaxt = "n")
> tt <- time(z)
> ix <- seq(1, length(tt), length.out=8)
> axis(side = 1, at = tt[ix], labels = FALSE)
> labs <- format(tt, "%Y-%b-%d")
> axis(side = 1, at = tt[ix], labels = labs[ix], tcl = -0.7, cex.axis = 0.7,
> las=2)
>
>
> It looks better, but I wouldn't risk "good".
>
> Rui Barradas
>
> Em 31-07-2012 19:10, Yolande Tra escreveu:
>
>> Thanks.
>> This is an irregular time series. The line plot does not look good because
>> of the gap between 9/11/2008 and 10/12/2009. I think two plots would be
>> better.
>> How would you include the date on the x-axis. Right now it only gives one
>> tick mark, 2009.
>>
>> Yolande
>>
>> On Tue, Jul 31, 2012 at 2:01 PM, Rui Barradas <[hidden email]>
>> wrote:
>>
>> Hello,
>>>
>>> It only gives that error if you don't
>>>
>>>
>>> last$date <- as.Date(last$date, format="%m/%d/%Y")
>>>
>>> You must have dates, not character values.
>>> Try it, then make a zoo object, then plot it.
>>>
>>> Rui Barradas
>>>
>>> Em 31-07-2012 18:54, Yolande Tra escreveu:
>>>
>>> Thank you everyone for the attempt to solve the problem
>>>
>>>> It is an irregular series and insert NAs when a date is missing
>>>>
>>>> library(zoo)
>>>>> z <- zoo(last$r_wvht, order.by=last$date)
>>>>> plot(z)
>>>>>
>>>>> Error in plot.window(...) : need finite 'xlim' values
>>>> In addition: Warning messages:
>>>> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
>>>> 2: In min(x) : no non-missing arguments to min; returning Inf
>>>> 3: In max(x) : no non-missing arguments to max; returning -Inf
>>>>
>>>> I have also tried
>>>>
>>>> z <- read.zoo("plot1.csv", header = TRUE, sep = ",", format =
>>>> "%m/%d/%Y")
>>>>
>>>>> plot(z)
>>>>>
>>>>> The plot does not look good. It does not display the individual dates.
>>>>
>>>> Y
>>>> Y
>>>> On Tue, Jul 31, 2012 at 1:33 PM, Rui Barradas <[hidden email]>
>>>> wrote:
>>>>
>>>> Hello,
>>>>
>>>>> Sorry, I forgot the "time series" part of your question. You could
>>>>> use
>>>>> instead one of
>>>>>
>>>>> # 1. type = "l" gives a line plot
>>>>> plot(r_wvht ~ date, data = last, type="l")
>>>>>
>>>>> # 2. use time series object plot
>>>>> library(zoo)
>>>>> z <- zoo(last$r_wvht, order.by=last$date)
>>>>> plot(z)
>>>>>
>>>>> Rui Barradas
>>>>>
>>>>> Em 31-07-2012 18:27, Rui Barradas escreveu:
>>>>>
>>>>> Hello,
>>>>>
>>>>> You could use dput(), it's not your first post...
>>>>>>
>>>>>>
>>>>>> last <- structure(list(date = c("8/6/2008", "8/8/2008", "8/11/2008",
>>>>>> "8/13/2008", "8/14/2008", "8/18/2008", "8/20/2008", "8/27/2008",
>>>>>> "8/28/2008", "8/31/2008", "9/2/2008", "9/3/2008", "9/4/2008",
>>>>>> "9/5/2008", "9/8/2008", "9/11/2008", "10/12/2009", "10/14/2009",
>>>>>> "10/19/2009", "10/21/2009", "10/26/2009", "10/29/2009", "10/30/2009",
>>>>>> "11/2/2009", "11/3/2009", "11/4/2009", "11/11/2009", "11/13/2009"
>>>>>> ), r_wvht = c(0.9766667, 0.7733333, 1.4833333, 1.5766667, 1.39,
>>>>>> 0.78, 0.8383333, 1.77, 1.295, 2.41, 1.3166667, 1.3075, 1.39,
>>>>>> 1.6333333, 1.2416667, 1.395, 0.8633333, 2.79, 1.0325, 1.965,
>>>>>> 1.78, 1.5666667, 1.05, 1.4633333, 1.24, 1.0075, 1.605, 1.8475
>>>>>> )), .Names = c("date", "r_wvht"), class = "data.frame", row.names =
>>>>>> c("1",
>>>>>> "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13",
>>>>>> "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24",
>>>>>> "25", "26", "27", "28"))
>>>>>>
>>>>>> last$date <- as.Date(last$date, format="%m/%d/%Y")
>>>>>>
>>>>>> plot(r_wvht ~ date, data = last)
>>>>>>
>>>>>> (You were forgetting that last$date was a character vector, not class
>>>>>> Date.)
>>>>>>
>>>>>> Hope this helps,
>>>>>>
>>>>>> Rui Barradas
>>>>>>
>>>>>> Em 31-07-2012 18:02, Yolande Tra escreveu:
>>>>>>
>>>>>> Hello,
>>>>>>
>>>>>>> This should be pretty simple but I cannot get it right. Please
>>>>>>> point to
>>>>>>> the
>>>>>>> right code. Thanks.
>>>>>>>
>>>>>>> last <- read.csv(file.path(dataDir,"******plot1.csv"), as.is
>>>>>>>
>>>>>>>
>>>>>>> =T,stringsAsFactors
>>>>>>>
>>>>>>>> = FALSE)
>>>>>>>>
>>>>>>> last
>>>>>>>
>>>>>>>> date r_wvht
>>>>>>>>
>>>>>>> 1 8/6/2008 0.9766667
>>>>>>> 2 8/8/2008 0.7733333
>>>>>>> 3 8/11/2008 1.4833333
>>>>>>> 4 8/13/2008 1.5766667
>>>>>>> 5 8/14/2008 1.3900000
>>>>>>> 6 8/18/2008 0.7800000
>>>>>>> 7 8/20/2008 0.8383333
>>>>>>> 8 8/27/2008 1.7700000
>>>>>>> 9 8/28/2008 1.2950000
>>>>>>> 10 8/31/2008 2.4100000
>>>>>>> 11 9/2/2008 1.3166667
>>>>>>> 12 9/3/2008 1.3075000
>>>>>>> 13 9/4/2008 1.3900000
>>>>>>> 14 9/5/2008 1.6333333
>>>>>>> 15 9/8/2008 1.2416667
>>>>>>> 16 9/11/2008 1.3950000
>>>>>>> 17 10/12/2009 0.8633333
>>>>>>> 18 10/14/2009 2.7900000
>>>>>>> 19 10/19/2009 1.0325000
>>>>>>> 20 10/21/2009 1.9650000
>>>>>>> 21 10/26/2009 1.7800000
>>>>>>> 22 10/29/2009 1.5666667
>>>>>>> 23 10/30/2009 1.0500000
>>>>>>> 24 11/2/2009 1.4633333
>>>>>>> 25 11/3/2009 1.2400000
>>>>>>> 26 11/4/2009 1.0075000
>>>>>>> 27 11/11/2009 1.6050000
>>>>>>> 28 11/13/2009 1.8475000
>>>>>>>
>>>>>>> x<-as.vector(last$date)
>>>>>>>
>>>>>>>> x
>>>>>>>>
>>>>>>>> [1] "8/6/2008" "8/8/2008" "8/11/2008" "8/13/2008"
>>>>>>>> "8/14/2008"
>>>>>>>>
>>>>>>> "8/18/2008" "8/20/2008" "8/27/2008"
>>>>>>> [9] "8/28/2008" "8/31/2008" "9/2/2008" "9/3/2008" "9/4/2008"
>>>>>>> "9/5/2008" "9/8/2008" "9/11/2008"
>>>>>>> [17] "10/12/2009" "10/14/2009" "10/19/2009" "10/21/2009" "10/26/2009"
>>>>>>> "10/29/2009" "10/30/2009" "11/2/2009"
>>>>>>> [25] "11/3/2009" "11/4/2009" "11/11/2009" "11/13/2009"
>>>>>>>
>>>>>>> y<-as.vector(last$r_wvht)
>>>>>>>
>>>>>>>> y
>>>>>>>>
>>>>>>>> [1] 0.9766667 0.7733333 1.4833333 1.5766667 1.3900000 0.7800000
>>>>>>>>
>>>>>>> 0.8383333
>>>>>>> 1.7700000 1.2950000 2.4100000
>>>>>>> [11] 1.3166667 1.3075000 1.3900000 1.6333333 1.2416667 1.3950000
>>>>>>> 0.8633333
>>>>>>> 2.7900000 1.0325000 1.9650000
>>>>>>> [21] 1.7800000 1.5666667 1.0500000 1.4633333 1.2400000 1.0075000
>>>>>>> 1.6050000
>>>>>>> 1.8475000
>>>>>>>
>>>>>>> plot(x,y)
>>>>>>>
>>>>>>>> Error in plot.window(...) : need finite 'xlim' values
>>>>>>>>
>>>>>>> In addition: Warning messages:
>>>>>>> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by
>>>>>>> coercion
>>>>>>> 2: In min(x) : no non-missing arguments to min; returning Inf
>>>>>>> 3: In max(x) : no non-missing arguments to max; returning -Inf
>>>>>>>
>>>>>>> plot(x,y,xlim=c("6/8/2008","******11/13/2009"))
>>>>>>>
>>>>>>> Error in plot.window(...) : invalid 'xlim' value
>>>>>>>>
>>>>>>> Y
>>>>>>>
>>>>>>> [[alternative HTML version deleted]]
>>>>>>>
>>>>>>> ______________________________******________________
>>>>>>> [hidden email] 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/**<http://www.r-project.org/**>
>>>>>>> <ht**tp://www.r-project.org/** <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.
>>>>>>>
>>>>>>> ______________________________******________________
>>>>>>>
>>>>>> [hidden email] 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/**<http://www.r-project.org/**>
>>>>>> <ht**tp://www.r-project.org/** <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.
>>>>>>
>>>>>>
>>>>>>
>

Rui Barradas

Re: time series line plot: Error in plot.window(...) : invalid 'xlim' value

In reply to this post by Rui Barradas

Hello,

Inline.
Em 31-07-2012 21:45, Yolande Tra escreveu:

> I am sorry but I still got some errors. May be it is because of the package
> zoo. Please be patient with me. Thanks.
>
> Yolande
>
>> last1 <- read.csv(file.path(dataDir,"plot1_2008.csv"), as.is=T,stringsAsFactors
> = FALSE)
>> library(zoo)
> Attaching package: ‘zoo’
> The following object(s) are masked from ‘package:base’:
> as.Date, as.Date.numeric

Before creating the zoo object, you MUST change the class of last1$date
to Date:

last1$date <- as.Date(last1$Date, format="%m/%d/%Y")

If not, it's a character vector, NOT a date.
And plot.zoo doesn't recognize valid x axis values.

Rui Barradas

>> u<- zoo(last1$r_wvht, order.by=last1$date)
>> plot(u, xaxt="n")
> Error in plot.window(...) : need finite 'xlim' values
> In addition: Warning messages:
> 1: In min(x) : no non-missing arguments to min; returning Inf
> 2: In max(x) : no non-missing arguments to max; returning -Inf
> 3: In min(x) : no non-missing arguments to min; returning Inf
> 4: In max(x) : no non-missing arguments to max; returning -Inf
> 5: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
> 6: In min(x) : no non-missing arguments to min; returning Inf
> 7: In max(x) : no non-missing arguments to max; returning -Inf
>
>
> On Tue, Jul 31, 2012 at 3:15 PM, Rui Barradas <[hidden email]> wrote:
>
>> Hello,
>>
>> Adapted from the help page for plot.zoo
>>
>> #library(zoo)
>> #z <- zoo(last$r_wvht, order.by=last$date)
>> plot(z, xaxt = "n")
>> tt <- time(z)
>> ix <- seq(1, length(tt), length.out=8)
>> axis(side = 1, at = tt[ix], labels = FALSE)
>> labs <- format(tt, "%Y-%b-%d")
>> axis(side = 1, at = tt[ix], labels = labs[ix], tcl = -0.7, cex.axis = 0.7,
>> las=2)
>>
>>
>> It looks better, but I wouldn't risk "good".
>>
>> Rui Barradas
>>
>> Em 31-07-2012 19:10, Yolande Tra escreveu:
>>
>>> Thanks.
>>> This is an irregular time series. The line plot does not look good because
>>> of the gap between 9/11/2008 and 10/12/2009. I think two plots would be
>>> better.
>>> How would you include the date on the x-axis. Right now it only gives one
>>> tick mark, 2009.
>>>
>>> Yolande
>>>
>>> On Tue, Jul 31, 2012 at 2:01 PM, Rui Barradas <[hidden email]>
>>> wrote:
>>>
>>> Hello,
>>>> It only gives that error if you don't
>>>>
>>>>
>>>> last$date <- as.Date(last$date, format="%m/%d/%Y")
>>>>
>>>> You must have dates, not character values.
>>>> Try it, then make a zoo object, then plot it.
>>>>
>>>> Rui Barradas
>>>>
>>>> Em 31-07-2012 18:54, Yolande Tra escreveu:
>>>>
>>>> Thank you everyone for the attempt to solve the problem
>>>>
>>>>> It is an irregular series and insert NAs when a date is missing
>>>>>
>>>>> library(zoo)
>>>>>> z <- zoo(last$r_wvht, order.by=last$date)
>>>>>> plot(z)
>>>>>>
>>>>>> Error in plot.window(...) : need finite 'xlim' values
>>>>> In addition: Warning messages:
>>>>> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
>>>>> 2: In min(x) : no non-missing arguments to min; returning Inf
>>>>> 3: In max(x) : no non-missing arguments to max; returning -Inf
>>>>>
>>>>> I have also tried
>>>>>
>>>>> z <- read.zoo("plot1.csv", header = TRUE, sep = ",", format =
>>>>> "%m/%d/%Y")
>>>>>
>>>>>> plot(z)
>>>>>>
>>>>>> The plot does not look good. It does not display the individual dates.
>>>>> Y
>>>>> Y
>>>>> On Tue, Jul 31, 2012 at 1:33 PM, Rui Barradas <[hidden email]>
>>>>> wrote:
>>>>>
>>>>> Hello,
>>>>>
>>>>>> Sorry, I forgot the "time series" part of your question. You could
>>>>>> use
>>>>>> instead one of
>>>>>>
>>>>>> # 1. type = "l" gives a line plot
>>>>>> plot(r_wvht ~ date, data = last, type="l")
>>>>>>
>>>>>> # 2. use time series object plot
>>>>>> library(zoo)
>>>>>> z <- zoo(last$r_wvht, order.by=last$date)
>>>>>> plot(z)
>>>>>>
>>>>>> Rui Barradas
>>>>>>
>>>>>> Em 31-07-2012 18:27, Rui Barradas escreveu:
>>>>>>
>>>>>> Hello,
>>>>>>
>>>>>> You could use dput(), it's not your first post...
>>>>>>>
>>>>>>> last <- structure(list(date = c("8/6/2008", "8/8/2008", "8/11/2008",
>>>>>>> "8/13/2008", "8/14/2008", "8/18/2008", "8/20/2008", "8/27/2008",
>>>>>>> "8/28/2008", "8/31/2008", "9/2/2008", "9/3/2008", "9/4/2008",
>>>>>>> "9/5/2008", "9/8/2008", "9/11/2008", "10/12/2009", "10/14/2009",
>>>>>>> "10/19/2009", "10/21/2009", "10/26/2009", "10/29/2009", "10/30/2009",
>>>>>>> "11/2/2009", "11/3/2009", "11/4/2009", "11/11/2009", "11/13/2009"
>>>>>>> ), r_wvht = c(0.9766667, 0.7733333, 1.4833333, 1.5766667, 1.39,
>>>>>>> 0.78, 0.8383333, 1.77, 1.295, 2.41, 1.3166667, 1.3075, 1.39,
>>>>>>> 1.6333333, 1.2416667, 1.395, 0.8633333, 2.79, 1.0325, 1.965,
>>>>>>> 1.78, 1.5666667, 1.05, 1.4633333, 1.24, 1.0075, 1.605, 1.8475
>>>>>>> )), .Names = c("date", "r_wvht"), class = "data.frame", row.names =
>>>>>>> c("1",
>>>>>>> "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13",
>>>>>>> "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24",
>>>>>>> "25", "26", "27", "28"))
>>>>>>>
>>>>>>> last$date <- as.Date(last$date, format="%m/%d/%Y")
>>>>>>>
>>>>>>> plot(r_wvht ~ date, data = last)
>>>>>>>
>>>>>>> (You were forgetting that last$date was a character vector, not class
>>>>>>> Date.)
>>>>>>>
>>>>>>> Hope this helps,
>>>>>>>
>>>>>>> Rui Barradas
>>>>>>>
>>>>>>> Em 31-07-2012 18:02, Yolande Tra escreveu:
>>>>>>>
>>>>>>> Hello,
>>>>>>>
>>>>>>>> This should be pretty simple but I cannot get it right. Please
>>>>>>>> point to
>>>>>>>> the
>>>>>>>> right code. Thanks.
>>>>>>>>
>>>>>>>> last <- read.csv(file.path(dataDir,"******plot1.csv"), as.is
>>>>>>>>
>>>>>>>>
>>>>>>>> =T,stringsAsFactors
>>>>>>>>
>>>>>>>>> = FALSE)
>>>>>>>>>
>>>>>>>> last
>>>>>>>>
>>>>>>>>> date r_wvht
>>>>>>>>>
>>>>>>>> 1 8/6/2008 0.9766667
>>>>>>>> 2 8/8/2008 0.7733333
>>>>>>>> 3 8/11/2008 1.4833333
>>>>>>>> 4 8/13/2008 1.5766667
>>>>>>>> 5 8/14/2008 1.3900000
>>>>>>>> 6 8/18/2008 0.7800000
>>>>>>>> 7 8/20/2008 0.8383333
>>>>>>>> 8 8/27/2008 1.7700000
>>>>>>>> 9 8/28/2008 1.2950000
>>>>>>>> 10 8/31/2008 2.4100000
>>>>>>>> 11 9/2/2008 1.3166667
>>>>>>>> 12 9/3/2008 1.3075000
>>>>>>>> 13 9/4/2008 1.3900000
>>>>>>>> 14 9/5/2008 1.6333333
>>>>>>>> 15 9/8/2008 1.2416667
>>>>>>>> 16 9/11/2008 1.3950000
>>>>>>>> 17 10/12/2009 0.8633333
>>>>>>>> 18 10/14/2009 2.7900000
>>>>>>>> 19 10/19/2009 1.0325000
>>>>>>>> 20 10/21/2009 1.9650000
>>>>>>>> 21 10/26/2009 1.7800000
>>>>>>>> 22 10/29/2009 1.5666667
>>>>>>>> 23 10/30/2009 1.0500000
>>>>>>>> 24 11/2/2009 1.4633333
>>>>>>>> 25 11/3/2009 1.2400000
>>>>>>>> 26 11/4/2009 1.0075000
>>>>>>>> 27 11/11/2009 1.6050000
>>>>>>>> 28 11/13/2009 1.8475000
>>>>>>>>
>>>>>>>> x<-as.vector(last$date)
>>>>>>>>
>>>>>>>>> x
>>>>>>>>>
>>>>>>>>> [1] "8/6/2008" "8/8/2008" "8/11/2008" "8/13/2008"
>>>>>>>>> "8/14/2008"
>>>>>>>>>
>>>>>>>> "8/18/2008" "8/20/2008" "8/27/2008"
>>>>>>>> [9] "8/28/2008" "8/31/2008" "9/2/2008" "9/3/2008" "9/4/2008"
>>>>>>>> "9/5/2008" "9/8/2008" "9/11/2008"
>>>>>>>> [17] "10/12/2009" "10/14/2009" "10/19/2009" "10/21/2009" "10/26/2009"
>>>>>>>> "10/29/2009" "10/30/2009" "11/2/2009"
>>>>>>>> [25] "11/3/2009" "11/4/2009" "11/11/2009" "11/13/2009"
>>>>>>>>
>>>>>>>> y<-as.vector(last$r_wvht)
>>>>>>>>
>>>>>>>>> y
>>>>>>>>>
>>>>>>>>> [1] 0.9766667 0.7733333 1.4833333 1.5766667 1.3900000 0.7800000
>>>>>>>>>
>>>>>>>> 0.8383333
>>>>>>>> 1.7700000 1.2950000 2.4100000
>>>>>>>> [11] 1.3166667 1.3075000 1.3900000 1.6333333 1.2416667 1.3950000
>>>>>>>> 0.8633333
>>>>>>>> 2.7900000 1.0325000 1.9650000
>>>>>>>> [21] 1.7800000 1.5666667 1.0500000 1.4633333 1.2400000 1.0075000
>>>>>>>> 1.6050000
>>>>>>>> 1.8475000
>>>>>>>>
>>>>>>>> plot(x,y)
>>>>>>>>
>>>>>>>>> Error in plot.window(...) : need finite 'xlim' values
>>>>>>>>>
>>>>>>>> In addition: Warning messages:
>>>>>>>> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by
>>>>>>>> coercion
>>>>>>>> 2: In min(x) : no non-missing arguments to min; returning Inf
>>>>>>>> 3: In max(x) : no non-missing arguments to max; returning -Inf
>>>>>>>>
>>>>>>>> plot(x,y,xlim=c("6/8/2008","******11/13/2009"))
>>>>>>>>
>>>>>>>> Error in plot.window(...) : invalid 'xlim' value
>>>>>>>> Y
>>>>>>>>
>>>>>>>> [[alternative HTML version deleted]]
>>>>>>>>
>>>>>>>> ______________________________******________________
>>>>>>>> [hidden email] 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/**<http://www.r-project.org/**>
>>>>>>>> <ht**tp://www.r-project.org/** <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.
>>>>>>>>
>>>>>>>> ______________________________******________________
>>>>>>>>
>>>>>>> [hidden email] 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/**<http://www.r-project.org/**>
>>>>>>> <ht**tp://www.r-project.org/** <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.
>>>>>>>
>>>>>>>
>>>>>>>

ytra

Re: time series line plot: Error in plot.window(...) : invalid 'xlim' value

It is not working yet

> last1$date <- as.Date(last1$date, format="%m/%d/%Y")
> u<- zoo(last1$r_wvht, order.by=last1$date)
> plot(u, xaxt="n")
Error in plot.window(...) : need finite 'ylim' values
In addition: Warning messages:
1: In min(x) : no non-missing arguments to min; returning Inf
2: In max(x) : no non-missing arguments to max; returning -Inf
3: In min(x) : no non-missing arguments to min; returning Inf
4: In max(x) : no non-missing arguments to max; returning -Inf
Instead I read the help page on zoo (as you pointed out) and create a zoo
object as follows

>z <-
zoo(c(0.98,0.77,1.48,1.58,1.39,0.78,0.84,1.77,1.30,2.41,1.32,1.31,1.39,1.63,1.24,1.40),

as.Date(c("2008-08-06","2008-08-08","2008-08-11","2008-08-13","2008-08-14","2008-08-18","2008-08-20",
"2008-08-27", "2008-08-28", "2008-08-31","2008-09-02",
"2008-09-03","2008-09-04","2008-09-05","2008-09-08","2008-09-11")))
plot(z, xaxt="n")
> plot(z, xaxt="n",xlab="",ylab="")

And now there was no warning.
Thanks,
Yolande

On Tue, Jul 31, 2012 at 5:06 PM, Rui Barradas <[hidden email]> wrote:

> Hello,
>
> Inline.
> Em 31-07-2012 21:45, Yolande Tra escreveu:
>
> I am sorry but I still got some errors. May be it is because of the package
>> zoo. Please be patient with me. Thanks.
>>
>> Yolande
>>
>> last1 <- read.csv(file.path(dataDir,"**plot1_2008.csv"), as.is
>>> =T,stringsAsFactors
>>>
>> = FALSE)
>>
>>> library(zoo)
>>>
>> Attaching package: zoo
>> The following object(s) are masked from package:base:
>> as.Date, as.Date.numeric
>>
>
> Before creating the zoo object, you MUST change the class of last1$date to
> Date:
>
> last1$date <- as.Date(last1$Date, format="%m/%d/%Y")
>
> If not, it's a character vector, NOT a date.
> And plot.zoo doesn't recognize valid x axis values.
>
> Rui Barradas
>
>> u<- zoo(last1$r_wvht, order.by=last1$date)
>>> plot(u, xaxt="n")
>>>
>> Error in plot.window(...) : need finite 'xlim' values
>> In addition: Warning messages:
>> 1: In min(x) : no non-missing arguments to min; returning Inf
>> 2: In max(x) : no non-missing arguments to max; returning -Inf
>> 3: In min(x) : no non-missing arguments to min; returning Inf
>> 4: In max(x) : no non-missing arguments to max; returning -Inf
>> 5: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
>> 6: In min(x) : no non-missing arguments to min; returning Inf
>> 7: In max(x) : no non-missing arguments to max; returning -Inf
>>
>>
>> On Tue, Jul 31, 2012 at 3:15 PM, Rui Barradas <[hidden email]>
>> wrote:
>>
>> Hello,
>>>
>>> Adapted from the help page for plot.zoo
>>>
>>> #library(zoo)
>>> #z <- zoo(last$r_wvht, order.by=last$date)
>>> plot(z, xaxt = "n")
>>> tt <- time(z)
>>> ix <- seq(1, length(tt), length.out=8)
>>> axis(side = 1, at = tt[ix], labels = FALSE)
>>> labs <- format(tt, "%Y-%b-%d")
>>> axis(side = 1, at = tt[ix], labels = labs[ix], tcl = -0.7, cex.axis =
>>> 0.7,
>>> las=2)
>>>
>>>
>>> It looks better, but I wouldn't risk "good".
>>>
>>> Rui Barradas
>>>
>>> Em 31-07-2012 19:10, Yolande Tra escreveu:
>>>
>>> Thanks.
>>>> This is an irregular time series. The line plot does not look good
>>>> because
>>>> of the gap between 9/11/2008 and 10/12/2009. I think two plots would be
>>>> better.
>>>> How would you include the date on the x-axis. Right now it only gives
>>>> one
>>>> tick mark, 2009.
>>>>
>>>> Yolande
>>>>
>>>> On Tue, Jul 31, 2012 at 2:01 PM, Rui Barradas <[hidden email]>
>>>> wrote:
>>>>
>>>> Hello,
>>>>
>>>>> It only gives that error if you don't
>>>>>
>>>>>
>>>>> last$date <- as.Date(last$date, format="%m/%d/%Y")
>>>>>
>>>>> You must have dates, not character values.
>>>>> Try it, then make a zoo object, then plot it.
>>>>>
>>>>> Rui Barradas
>>>>>
>>>>> Em 31-07-2012 18:54, Yolande Tra escreveu:
>>>>>
>>>>> Thank you everyone for the attempt to solve the problem
>>>>>
>>>>> It is an irregular series and insert NAs when a date is missing
>>>>>>
>>>>>> library(zoo)
>>>>>>
>>>>>>> z <- zoo(last$r_wvht, order.by=last$date)
>>>>>>> plot(z)
>>>>>>>
>>>>>>> Error in plot.window(...) : need finite 'xlim' values
>>>>>>>
>>>>>> In addition: Warning messages:
>>>>>> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by
>>>>>> coercion
>>>>>> 2: In min(x) : no non-missing arguments to min; returning Inf
>>>>>> 3: In max(x) : no non-missing arguments to max; returning -Inf
>>>>>>
>>>>>> I have also tried
>>>>>>
>>>>>> z <- read.zoo("plot1.csv", header = TRUE, sep = ",", format =
>>>>>> "%m/%d/%Y")
>>>>>>
>>>>>> plot(z)
>>>>>>>
>>>>>>> The plot does not look good. It does not display the individual
>>>>>>> dates.
>>>>>>>
>>>>>> Y
>>>>>> Y
>>>>>> On Tue, Jul 31, 2012 at 1:33 PM, Rui Barradas <[hidden email]>
>>>>>> wrote:
>>>>>>
>>>>>> Hello,
>>>>>>
>>>>>> Sorry, I forgot the "time series" part of your question. You could
>>>>>>> use
>>>>>>> instead one of
>>>>>>>
>>>>>>> # 1. type = "l" gives a line plot
>>>>>>> plot(r_wvht ~ date, data = last, type="l")
>>>>>>>
>>>>>>> # 2. use time series object plot
>>>>>>> library(zoo)
>>>>>>> z <- zoo(last$r_wvht, order.by=last$date)
>>>>>>> plot(z)
>>>>>>>
>>>>>>> Rui Barradas
>>>>>>>
>>>>>>> Em 31-07-2012 18:27, Rui Barradas escreveu:
>>>>>>>
>>>>>>> Hello,
>>>>>>>
>>>>>>> You could use dput(), it's not your first post...
>>>>>>>
>>>>>>>>
>>>>>>>> last <- structure(list(date = c("8/6/2008", "8/8/2008", "8/11/2008",
>>>>>>>> "8/13/2008", "8/14/2008", "8/18/2008", "8/20/2008", "8/27/2008",
>>>>>>>> "8/28/2008", "8/31/2008", "9/2/2008", "9/3/2008", "9/4/2008",
>>>>>>>> "9/5/2008", "9/8/2008", "9/11/2008", "10/12/2009", "10/14/2009",
>>>>>>>> "10/19/2009", "10/21/2009", "10/26/2009", "10/29/2009",
>>>>>>>> "10/30/2009",
>>>>>>>> "11/2/2009", "11/3/2009", "11/4/2009", "11/11/2009", "11/13/2009"
>>>>>>>> ), r_wvht = c(0.9766667, 0.7733333, 1.4833333, 1.5766667, 1.39,
>>>>>>>> 0.78, 0.8383333, 1.77, 1.295, 2.41, 1.3166667, 1.3075, 1.39,
>>>>>>>> 1.6333333, 1.2416667, 1.395, 0.8633333, 2.79, 1.0325, 1.965,
>>>>>>>> 1.78, 1.5666667, 1.05, 1.4633333, 1.24, 1.0075, 1.605, 1.8475
>>>>>>>> )), .Names = c("date", "r_wvht"), class = "data.frame", row.names =
>>>>>>>> c("1",
>>>>>>>> "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13",
>>>>>>>> "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24",
>>>>>>>> "25", "26", "27", "28"))
>>>>>>>>
>>>>>>>> last$date <- as.Date(last$date, format="%m/%d/%Y")
>>>>>>>>
>>>>>>>> plot(r_wvht ~ date, data = last)
>>>>>>>>
>>>>>>>> (You were forgetting that last$date was a character vector, not
>>>>>>>> class
>>>>>>>> Date.)
>>>>>>>>
>>>>>>>> Hope this helps,
>>>>>>>>
>>>>>>>> Rui Barradas
>>>>>>>>
>>>>>>>> Em 31-07-2012 18:02, Yolande Tra escreveu:
>>>>>>>>
>>>>>>>> Hello,
>>>>>>>>
>>>>>>>> This should be pretty simple but I cannot get it right. Please
>>>>>>>>> point to
>>>>>>>>> the
>>>>>>>>> right code. Thanks.
>>>>>>>>>
>>>>>>>>> last <- read.csv(file.path(dataDir,"********plot1.csv"), as.is
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> =T,stringsAsFactors
>>>>>>>>>
>>>>>>>>> = FALSE)
>>>>>>>>>>
>>>>>>>>>> last
>>>>>>>>>
>>>>>>>>> date r_wvht
>>>>>>>>>>
>>>>>>>>>> 1 8/6/2008 0.9766667
>>>>>>>>> 2 8/8/2008 0.7733333
>>>>>>>>> 3 8/11/2008 1.4833333
>>>>>>>>> 4 8/13/2008 1.5766667
>>>>>>>>> 5 8/14/2008 1.3900000
>>>>>>>>> 6 8/18/2008 0.7800000
>>>>>>>>> 7 8/20/2008 0.8383333
>>>>>>>>> 8 8/27/2008 1.7700000
>>>>>>>>> 9 8/28/2008 1.2950000
>>>>>>>>> 10 8/31/2008 2.4100000
>>>>>>>>> 11 9/2/2008 1.3166667
>>>>>>>>> 12 9/3/2008 1.3075000
>>>>>>>>> 13 9/4/2008 1.3900000
>>>>>>>>> 14 9/5/2008 1.6333333
>>>>>>>>> 15 9/8/2008 1.2416667
>>>>>>>>> 16 9/11/2008 1.3950000
>>>>>>>>> 17 10/12/2009 0.8633333
>>>>>>>>> 18 10/14/2009 2.7900000
>>>>>>>>> 19 10/19/2009 1.0325000
>>>>>>>>> 20 10/21/2009 1.9650000
>>>>>>>>> 21 10/26/2009 1.7800000
>>>>>>>>> 22 10/29/2009 1.5666667
>>>>>>>>> 23 10/30/2009 1.0500000
>>>>>>>>> 24 11/2/2009 1.4633333
>>>>>>>>> 25 11/3/2009 1.2400000
>>>>>>>>> 26 11/4/2009 1.0075000
>>>>>>>>> 27 11/11/2009 1.6050000
>>>>>>>>> 28 11/13/2009 1.8475000
>>>>>>>>>
>>>>>>>>> x<-as.vector(last$date)
>>>>>>>>>
>>>>>>>>> x
>>>>>>>>>>
>>>>>>>>>> [1] "8/6/2008" "8/8/2008" "8/11/2008" "8/13/2008"
>>>>>>>>>> "8/14/2008"
>>>>>>>>>>
>>>>>>>>>> "8/18/2008" "8/20/2008" "8/27/2008"
>>>>>>>>> [9] "8/28/2008" "8/31/2008" "9/2/2008" "9/3/2008"
>>>>>>>>> "9/4/2008"
>>>>>>>>> "9/5/2008" "9/8/2008" "9/11/2008"
>>>>>>>>> [17] "10/12/2009" "10/14/2009" "10/19/2009" "10/21/2009"
>>>>>>>>> "10/26/2009"
>>>>>>>>> "10/29/2009" "10/30/2009" "11/2/2009"
>>>>>>>>> [25] "11/3/2009" "11/4/2009" "11/11/2009" "11/13/2009"
>>>>>>>>>
>>>>>>>>> y<-as.vector(last$r_wvht)
>>>>>>>>>
>>>>>>>>> y
>>>>>>>>>>
>>>>>>>>>> [1] 0.9766667 0.7733333 1.4833333 1.5766667 1.3900000
>>>>>>>>>> 0.7800000
>>>>>>>>>>
>>>>>>>>>> 0.8383333
>>>>>>>>> 1.7700000 1.2950000 2.4100000
>>>>>>>>> [11] 1.3166667 1.3075000 1.3900000 1.6333333 1.2416667 1.3950000
>>>>>>>>> 0.8633333
>>>>>>>>> 2.7900000 1.0325000 1.9650000
>>>>>>>>> [21] 1.7800000 1.5666667 1.0500000 1.4633333 1.2400000 1.0075000
>>>>>>>>> 1.6050000
>>>>>>>>> 1.8475000
>>>>>>>>>
>>>>>>>>> plot(x,y)
>>>>>>>>>
>>>>>>>>> Error in plot.window(...) : need finite 'xlim' values
>>>>>>>>>>
>>>>>>>>>> In addition: Warning messages:
>>>>>>>>> 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by
>>>>>>>>> coercion
>>>>>>>>> 2: In min(x) : no non-missing arguments to min; returning Inf
>>>>>>>>> 3: In max(x) : no non-missing arguments to max; returning -Inf
>>>>>>>>>
>>>>>>>>> plot(x,y,xlim=c("6/8/2008","********11/13/2009"))
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Error in plot.window(...) : invalid 'xlim' value
>>>>>>>>> Y
>>>>>>>>>
>>>>>>>>> [[alternative HTML version deleted]]
>>>>>>>>>
>>>>>>>>> ______________________________********________________
>>>>>>>>> [hidden email] 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://stat.ethz.ch/mailman/****listinfo/r-help>
>>>>>>>>> <ht**tps://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**://ethz.ch/mailman/listinfo/r-****help<http://ethz.ch/mailman/listinfo/r-**help>
>>>>>>>>> >
>>>>>>>>> <http**s://stat.ethz.ch/**mailman/**listinfo/r-help<http://stat.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/**<http://www.r-project.org/**>
>>>>>>>>> <ht**tp://www.r-project.org/** <http://www.r-project.org/**>>
>>>>>>>>> <ht**tp://www.r-project.org/** <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>
>>>>>>>>> >
>>>>>>>>> <http://www.**r-project.org/**posting-guide.**html<http://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.
>>>>>>>>>
>>>>>>>>> ______________________________********________________
>>>>>>>>>
>>>>>>>>> [hidden email] 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://stat.ethz.ch/mailman/****listinfo/r-help>
>>>>>>>> <ht**tps://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**://ethz.ch/mailman/listinfo/r-****help<http://ethz.ch/mailman/listinfo/r-**help>
>>>>>>>> >
>>>>>>>> <http**s://stat.ethz.ch/**mailman/**listinfo/r-help<http://stat.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/**<http://www.r-project.org/**>
>>>>>>>> <ht**tp://www.r-project.org/** <http://www.r-project.org/**>>
>>>>>>>> <ht**tp://www.r-project.org/** <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>
>>>>>>>> >
>>>>>>>> <http://www.**r-project.org/**posting-guide.**html<http://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.
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>>
>