median.default was changed between 2.7.1 and 2.8.1 to
call sum(...)/2 instead of mean(...) and that causes
the problem for POSIXct objects (sum fails but mean
works for them).
Bill Dunlap
TIBCO Software Inc - Spotfire Division
wdunlap tibco.com
> -----Original Message-----
> From: William Dunlap
> Sent: Tuesday, March 10, 2009 11:37 AM
> To: 'r-help@r-project.org'
> Subject: Re: [R] puzzled by math on date-time objects
>
> Your problem arises in R 2.8.1 (and 2.9.0-devel, but not 2.7.0) when
> length(POSIXct object) is even, because median(POSIXct object)
> passes a POSIXct object to median.default, which calls
> sum() in the even length case.
>
> > median( as.POSIXct(Sys.time()))
> [1] "2009-03-10 10:28:46 PDT"
> > median( as.POSIXct(rep(Sys.time(),2)))
> Error in Summary.POSIXct(c(1236706132.54740,
> 1236706132.54740), na.rm = FALSE) :
> 'sum' not defined for "POSIXt" objects
> > traceback()
> 4: stop(gettextf("'%s' not defined for \"POSIXt\" objects", .Generic),
> domain = NA)
> 3: Summary.POSIXct(c(1236706132.54740, 1236706132.54740),
> na.rm = FALSE)
> 2: median.default(as.POSIXct(rep(Sys.time(), 2)))
> 1: median(as.POSIXct(rep(Sys.time(), 2)))
> > version
> _
> platform i686-pc-linux-gnu
> arch i686
> os linux-gnu
> system i686, linux-gnu
> status
> major 2
> minor 8.1
> year 2008
> month 12
> day 22
> svn rev 47281
> language R
> version.string R version 2.8.1 (2008-12-22)
>
> Bill Dunlap
> TIBCO Software Inc - Spotfire Division
> wdunlap tibco.com
>
> ------------------------------------------------------------
> [R] puzzled by math on date-time objects
>
> Denis Chabot chabotd at globetrotter.net
> Tue Mar 10 16:44:07 CET 2009
> Previous message: [R] nonmetric clustering
> Next message: [R] perform subgroup meta-analysis and create
> forest plot displaying subgroups
> Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]
> Hi,
>
> I don't understand the following. When I create a small
> artificial set
> of date information in class POSIXct, I can calculate the
> mean and the
> median:
>
> a = as.POSIXct(Sys.time())
> a = a + 60*0:10; a
>
> [1] "2009-03-10 11:30:16 EDT" "2009-03-10 11:31:16 EDT"
> "2009-03-10
> 11:32:16 EDT"
> [4] "2009-03-10 11:33:16 EDT" "2009-03-10 11:34:16 EDT"
> "2009-03-10
> 11:35:16 EDT"
> [7] "2009-03-10 11:36:16 EDT" "2009-03-10 11:37:16 EDT"
> "2009-03-10
> 11:38:16 EDT"
> [10] "2009-03-10 11:39:16 EDT" "2009-03-10 11:40:16 EDT"
>
> median(a)
> [1] "2009-03-10 11:35:16 EDT"
> mean(a)
> [1] "2009-03-10 11:35:16 EDT"
>
>
> But for real data (for this post, a short subset is in object
> c) that
> I have converted into a POSIXct object, I cannot calculate
> the median
> with median(), though I do get it with summary():
>
> c
> [1] "2009-02-24 14:51:18 EST" "2009-02-24 14:51:19 EST"
> "2009-02-24
> 14:51:19 EST"
> [4] "2009-02-24 14:51:20 EST" "2009-02-24 14:51:20 EST"
> "2009-02-24
> 14:51:21 EST"
> [7] "2009-02-24 14:51:21 EST" "2009-02-24 14:51:22 EST"
> "2009-02-24
> 14:51:22 EST"
> [10] "2009-02-24 14:51:22 EST"
>
> class(c)
> [1] "POSIXt" "POSIXct"
>
> median(c)
> Erreur dans Summary.POSIXct(c(1235505080.6, 1235505081.1), na.rm =
> FALSE) :
> 'sum' not defined for "POSIXt" objects
>
> One difference is that in my own date-time series, some events are
> repeated (the original data contained fractions of seconds).
> But then,
> why can I get a median through summary()?
>
> summary(c)
> Min. 1st
> Qu. Median
> "2009-02-24 14:51:18 EST" "2009-02-24 14:51:19 EST" "2009-02-24
> 14:51:20 EST"
> Mean 3rd
> Qu. Max.
> "2009-02-24 14:51:20 EST" "2009-02-24 14:51:21 EST" "2009-02-24
> 14:51:22 EST"
>
> Thanks in advance,
>
>
> Denis Chabot
>
> sessionInfo()
> R version 2.8.1 Patched (2009-01-19 r47650)
> i386-apple-darwin9.6.0
>
> locale:
> fr_CA.UTF-8/fr_CA.UTF-8/C/C/fr_CA.UTF-8/fr_CA.UTF-8
>
> attached base packages:
> [1] stats graphics grDevices utils datasets methods base
>
> other attached packages:
> [1] doBy_3.7 chron_2.3-30
>
> loaded via a namespace (and not attached):
> [1] Hmisc_3.5-2 cluster_1.11.12 grid_2.8.1 lattice_0.17-20
> tools_2.8.1
>
______________________________________________
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.
