我在R中有两个相同长度的正值向量x,y.使用plot(x,y,"l",...),给出了我的有限向量x和y中的2维连续线图.有没有办法在R中的这个线图的某个范围内计算一个定积分？
edit1:我已经研究了R中的积分函数.但我不确定如何用两个向量传递函数来传递给它,因为我的向量都是有限的.
edit2:对于更多背景,x和y~10,000的长度.我写了一个函数来查找我正在观察的数据中的异常的周期[xi,xj].对于这些异常中的每一个,我都使用了情节来查看我的数据片段中发生了什么.现在我需要计算有关这些异常时期积分值的统计数据,所以我试图尽可能准确地得到与我的图形匹配的数字.X是一个时间变量,我花了很短的时间间隔.
你可以进行整合integrate().要创建一个功能你的矢量x和y,你需要的值之间进行插补.approxfun()就是这么做的.
integrate 需要一个函数和两个边界.
approxfun需要两个向量x,y就像你拥有的那样.
所以我的解决方案是:
integrate(approxfun(x,y), range(x)[1], range(x)[2])