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

Tue, Jul 31, 2012 12:15 PM

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 <ruipbarradas at sapo.pt> 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 <ruipbarradas at sapo.pt>
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]]

______________________________****________________
R-help at r-project.org mailing list
https://stat.ethz.ch/mailman/****listinfo/r-help<https://stat.ethz.ch/mailman/**listinfo/r-help>
<https://stat.**ethz.ch/mailman/listinfo/r-**help<https://stat.ethz.ch/mailman/listinfo/r-help>
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.

______________________________****________________

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.

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

Thread (9 messages)