我最近问了一个关于如何计算曲线下面积的问题:
How to calculate the area under each end of a sine curve
我收到了一些很好的答案,这些答案都涉及在线性回归模型中使用的预测函数。现在我想计算一个类似的曲线下面积,但是形成曲线的线是一个简单的平滑拟合(不是基于任何类型的回归分析)。
给定一个数据集:
y<-c(-13,16,35,40,28,36,43,33,40,33,22,-5,-27,-31,-29,-25,-26,-31,-26,-24,-25,-29,-23,4)
t<-1:24
我想为我的数据(上下两点的平均值)拟合平滑的线条,如下所示:
Smooth<-as.vector(filter(y,filter=rep(1,4))/4)
#I don't want to loose my first and last data points, so I cheat and add them back in.
Check<-data.frame(y,Smooth)
Smooth2<-ifelse(is.na(Check$Smooth),y, Check$Smooth)
plot(t,y)
abline(h=0,lty=2)
lines(t,Smooth2)
接下来,我想计算正曲线下面积以及负曲线下面积。我假设这将涉及找到3个根(线通过零的位置),描绘正曲线开始和停止的位置,负曲线开始和停止。