r - 计算R中线图的积分

Question

我在 R 中有两个长度相同的正值向量 x,y。使用 plot(x, y, "l",...)，在我的有限向量中给出一个二维连续线图x 和 y。有没有办法在 R 中的这条线图的某个范围内计算定积分？

edit1:我研究了 R 中的积分函数。但是我不确定如何从两个向量中创建一个函数以传递给它，因为我的向量都是有限的。

edit2:对于更多背景，x 和 y 的长度 ~ 10,000。我编写了一个函数来查找我正在观察的数据中异常的周期 [xi, xj]。对于这些异常中的每一个，我都使用 plot 来查看这些数据片段中发生了什么。现在我需要计算有关这些异常时期积分值的统计数据，因此我试图获得尽可能准确的数字以与我的图表匹配。 X 是一个时间变量，我已经采取了非常好的时间间隔。

Answer 1

您可以使用 integrate() 进行集成。要从向量 x 和 y 创建函数，您需要在值之间进行插值。 approxfun() 正是这样做的。

integrate 接受一个函数和两个边界。

approxfun 采用两个向量 x 和 y 就像你拥有的那样。

所以我的解决方案是:

integrate(approxfun(x,y), range(x)[1], range(x)[2])