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