r - 在 R 中积分非参数曲线

Question

只是一个警告，我一天前开始使用 R ......如果有什么事情看起来很简单，我很抱歉。

现在我正在尝试让 R 接收一个带有撞击加速度计数据的 .txt 文件，并为其计算头部损伤标准测试。 HIC 测试要求数据的曲线在某个时间间隔内进行积分。

方程式在下面的链接中......我试图将它作为图像插入这里，但它不会让我。显然我需要一些声望点才能让我这样做。



a(t) 是加速度曲线。

到目前为止，我还没有遇到在 R 中生成合适的曲线以匹配数据的问题。 loess 函数工作得很好，正是我正在寻找的那种东西……我只是不知道如何整合它。据我所知，loess 是一种非参数回归，因此无法确定曲线本身的方程。有没有办法整合它？

如果没有，是否有另一种方法可以使用不同的函数完成此任务？

任何帮助或有见地的评论将不胜感激。

提前致谢，

韦斯

詹姆斯还有一个问题，我怎样才能在使用集成()函数时获得没有文本和错误的数字？

Answer 1

您可以使用 predict您的 loess 上的功能模型来创建一个与 integrate 一起使用的函数.

# using the inbuilt dataset "pressure"
plot(pressure,type="l")

# create loess object and prediction function
l <- loess(pressure~temperature,pressure)
f <- function(x) predict(l,newdata=x)

# perform integration
integrate(f,0,360)
40176.5 with absolute error < 4.6

并单独提取值:

integrate(f,0,360)$value
[1] 40176.5