如何在R语言中获取两条交叉ECDF曲线的交点坐标？

ECDF是阶梯状的累积分布函数，两条ECDF的交点可能出现在阶梯的水平段或垂直过渡段，以下是具体实现步骤：

1. 提取ECDF核心数据
   从生成的ecdf对象中直接提取断点x值和对应的累积概率y值：

   # 提取两条ECDF的断点x和累积概率y
   x1 <- ecdf1$x
   y1 <- ecdf1$y
   x2 <- ecdf2$x
   y2 <- ecdf2$y
   

2. 生成密集评估点
   合并两条ECDF的所有断点并排序，再在相邻断点间插入密集中间点，避免漏掉水平段上的交点：

   # 合并并排序所有断点
   all_x <- sort(unique(c(x1, x2)))
   # 在相邻断点间插入100个中间点，生成密集x序列
   dense_x <- unlist(lapply(1:(length(all_x)-1), function(i) {
     seq(all_x[i], all_x[i+1], length.out = 100)
   }))
   dense_x <- c(dense_x, all_x[length(all_x)])
   

3. 计算密集点对应的ECDF值
   用ecdf函数计算所有密集x点对应的累积概率：

   y1_dense <- ecdf1(dense_x)
   y2_dense <- ecdf2(dense_x)
   

4. 定位并计算交点坐标
   通过判断两条ECDF值的差值符号变化找到交点区间，再根据阶梯函数特性计算坐标：

   # 计算差值，找出符号变化的位置
   diff_y <- y1_dense - y2_dense
   sign_changes <- which(diff(sign(diff_y)) != 0)
   
   # 逐个计算交点坐标
   intersections <- lapply(sign_changes, function(i) {
     x_left <- dense_x[i]
     x_right <- dense_x[i+1]
     y_l1 <- y1_dense[i]
     y_l2 <- y2_dense[i]
     y_r1 <- y1_dense[i+1]
     y_r2 <- y2_dense[i+1]
   
     # 情况1：第一条ECDF水平，第二条垂直上升
     if (y_l1 == y_r1 && y_l2 != y_r2) {
       data.frame(x = x_left, y = y_l1)
     # 情况2：第二条ECDF水平，第一条垂直上升
     } else if (y_l2 == y_r2 && y_l1 != y_r1) {
       data.frame(x = x_left, y = y_l2)
     # 情况3：两条ECDF同时阶梯变化，线性插值计算交点
     } else {
       slope1 <- (y_r1 - y_l1)/(x_right - x_left)
       slope2 <- (y_r2 - y_l2)/(x_right - x_left)
       x_inter <- (y_l2 - y_l1 + slope1*x_left - slope2*x_left)/(slope1 - slope2)
       y_inter <- y_l1 + slope1*(x_inter - x_left)
       data.frame(x = x_inter, y = y_inter)
     }
   })
   
   # 合并所有交点为数据框
   intersections_df <- do.call(rbind, intersections)
   

5. 验证结果（可选）
   绘制ECDF曲线并标记交点，确认结果正确性：

   plot(ecdf1, col = "red", main = "ECDF曲线与交点")
   lines(ecdf2, col = "blue")
   points(intersections_df, col = "darkgreen", pch = 16, cex = 1.2)
   

注意事项 #

• 如果数据量极大，可减少插入的中间点数量（比如将length.out设为50），在保证精度的同时提升效率。
• 该方法可捕捉所有交点（包括多次交叉的情况）。

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