搜索以找到两个覆盖曲线下面积95％的数字

Question

背景：

我有两列数字（列“ dd ”和列“ yy ”， 请参见下文以访问此数据 ）。 当我将“ dd” （ X轴 ）相对于“ yy ”（ Y轴 ）作图时，请使用：

plot(dd, yy, ty="l", xlab = "dd", ylab ="yy", lwd=2, xaxs="i")

我得到这个：

编码问题：

从“ -Inf”变为“ 0”， 直到还有其他点覆盖此曲线下面积的95％？ 换句话说，如果箭头要覆盖此曲线下95％的面积，我们想查看与两个箭头提示相对应的两个数字吗？

这是R代码，可以轻松地将这些小数据导入Rstudio：

id <- "0B5V8AyEFBTmXM1VIYUYxSG5tSjQ"

Data <- read.csv(paste0("https://docs.google.com/uc?id=",id,"&export=download")) 

Answer 1

编辑：我以前的近似不是很准确。 原来， pracma包具有这样的功能， cumtrapz ，近似于积分与梯形，并在从左边开始的每个点输出的总和。 使用该功能可以提供更准确的答案。

library(pracma)
points <- read.csv(paste0("https://docs.google.com/uc?id=",id,"&export=download")) 
area   <- cumtrapz(points$dd, points$yy)
ind    <- rev(which(area[length(area)] - area > .95))[1]

point  <- points$dd[ind]
point

# [1] -0.06406406