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

Question

背景

我有两列数字（列＆＃34; dd ＆＃34;和列＆＃34; yy ＆＃34;，请参阅下文以访问此数据）。当我将＆＃34; dd＆＃34; （ X轴）与＆＃34; yy ＆＃34; （ Y轴），使用：

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

我明白了：

编码问题：

从＆＃34; -Inf＆＃34; 转移到＆＃34; 0＆＃34; u p，直到其他任何一点覆盖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