在R中针对lubridate格式时长计算ECDF的问题求助

核心思路 #

ecdf()函数仅支持数值型输入，因此需要先将lubridate的分秒时长转换为总秒数，再通过ECDF函数或直接统计比例得到目标结论。

步骤1：将时长转换为数值型（总秒数） #

使用lubridate::period_to_seconds()把分秒格式的时长转为总秒数，这样就能用常规数值方法计算：

library(lubridate)

# 假设你的数据框是dat
dat$Value_sec <- period_to_seconds(dat$Value)

步骤2：计算特定个体的ECDF对应比例 #

以计算A的排名为例，有两种等价方法：

方法1：使用ECDF函数

# 创建ECDF函数
pace_ecdf <- ecdf(dat$Value_sec)

# 获取A对应的总秒数
a_sec <- dat$Value_sec[dat$Name == "A"]

# 计算：比A慢（时长更长）的参与者比例，即A超过的人数占比
p_faster <- 1 - pace_ecdf(a_sec)

# 格式化输出结论
cat(sprintf("A比%.0f%%的测试参与者更快", p_faster * 100))

方法2：直接统计比例

无需创建ECDF函数，直接用均值计算时长大于A的比例：

a_sec <- dat$Value_sec[dat$Name == "A"]
p_faster <- mean(dat$Value_sec > a_sec)

cat(sprintf("A比%.0f%%的测试参与者更快", p_faster * 100))

关键说明 #

• 时长越短速度越快，因此“比A更快”的是时长≤A的参与者，而我们要的“A比XX%的人更快”对应的是时长>A的参与者比例。
• 如果需要包含时长等于A的情况，把>换成>=即可，结论改为“A至少比XX%的测试参与者更快”。

内容的提问来源于stack exchange，提问作者csarvf_01