[R] time series line plot: Error in plot.window(...) : invalid 'xlim' value
Rui Barradas
ruipbarradas at sapo.pt
Tue Jul 31 21:15:25 CEST 2012
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.
>>>>>>
>>>>>> ______________________________****________________
>>>>> 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.
>>>>>
>>>>>
More information about the R-help
mailing list