[R] Hist with Dates

Jeff Newmiller jdnewmil at dcn.davis.CA.us
Tue Jan 15 18:30:59 CET 2013


You should read the help for hist carefully:

?hist.Date

hist has very specific rules for applying breaks that are often different than tabulation (tapply or aggregate, with or without using the cut function first).
---------------------------------------------------------------------------
Jeff Newmiller                        The     .....       .....  Go Live...
DCN:<jdnewmil at dcn.davis.ca.us>        Basics: ##.#.       ##.#.  Live Go...
                                      Live:   OO#.. Dead: OO#..  Playing
Research Engineer (Solar/Batteries            O.O#.       #.O#.  with
/Software/Embedded Controllers)               .OO#.       .OO#.  rocks...1k
--------------------------------------------------------------------------- 
Sent from my phone. Please excuse my brevity.

Dominic Roye <dominic.roye at gmail.com> wrote:

>Hello R-Helpers,
>
>
>I would like to generate the frequencies of dates (data.RData).
>
>str(data)
>Date[1:43748], format: "2001-01-01" "2001-01-01" "2001-01-01"
>"2001-01-01"
>...
>
>
>test <- hist(data,"days",freq=T)
>> test$counts
>[1] 27 17 12 15  6  4  5 13  5 13 15  9  7 10 11 21 10 10  9 11 15  7
>10
>14
>
>
>but the count is wrong, because i have only 9 times 01/01/2001 etc.
>
>> data[1:30]
> [1] "2001-01-01" "2001-01-01" "2001-01-01" "2001-01-01" "2001-01-01"
> [6] "2001-01-01" "2001-01-01" "2001-01-01" "2001-01-01"
>
>"2001-01-02"
>[11] "2001-01-02" "2001-01-02" "2001-01-02" "2001-01-02" "2001-01-02"
>[16] "2001-01-02" "2001-01-02" "2001-01-02" "2001-01-02" "2001-01-02"
>[21] "2001-01-03" "2001-01-03" "2001-01-03" "2001-01-03" "2001-01-03"
>[26] "2001-01-03" "2001-01-03" "2001-01-03" "2001-01-03" "2001-01-03"
>
>
>
>Now, when i give R the breaks.
>
>> br <-
>seq(as.Date("31.12.2000",format="%d.%m.%Y"),as.Date("01.01.2012",format="%d.%m.%Y"),by=1)
>> test2<- hist(data,br,freq=T)
>> test2$counts[1:10]
> [1] 12 15 17 12 15  6  4  5 13  5
>
>
>In this case, neither get the right counts but limiting i get the
>correct
>counts:
>
>> test3<- hist(data[1:50],br[1:6],freq=T)
>> test3$counts
>[1]  9 11 13  9  8
>
>
>I hope someone can help me with this problem.
>
>
>Thank you very much!!
>
>Best regards,
>
>Dominic
>
>
>------------------------------------------------------------------------
>
>______________________________________________
>R-help at 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.



More information about the R-help mailing list