[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