日期向量的分位数函数

我注意到基础 Rquantile函数不支持日期参数。

我理解为日期定义分位数需要在定义中小心（即，如果您有 6 个日期并要求第 25 个百分位数，则需要定义合适的舍入）。

是否有这样的分位数函数的有效实现，无论是作为基础的一部分还是另一个包。

下面的示例函数基本上实现了我感兴趣的功能（通过一些调整来处理第 0 个百分位数的情况），但我认为更有效的实现是可能的。

#Date quantile function.
dquantile <- function(x, probs){

  sx <- sort(x)

  pos <- round( probs * length(x) )

  return( sx[pos] )
}

# Example.
dates <- as.Date("01/01/1900", "%d/%m/%Y") + floor( 36500 * runif(100000) )

dquantile(dates, c(0.001, 0.025, 0.975, 0.999) )

The quantile函数确实支持日期，您只需指定type争论。您的问题可以通过以下方式解决：

dates <- as.Date("01/01/1900", "%d/%m/%Y") + floor( 36500 * runif(100000) )

quantile(dates, probs = c(0.001, 0.025, 0.975, 0.999), type = 1)

        0.1%         2.5%        97.5%        99.9% 
"1900-02-04" "1902-06-23" "1997-06-10" "1999-10-30"