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 <ruipbarra...@sapo.pt> 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 <ruipbarra...@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 <ruipbarra...@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 <ruipbarra...@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@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://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.
>>>>>>>>>
>>>>>>>>> ______________________________********________________
>>>>>>>>>
>>>>>>>>> R-help@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://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.
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>>
>
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org 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.
