当前位置：文江博客话题详情

在 R 中积分非参数曲线

发布于 2024-12-10 21:33:33 字数 546 浏览 0 评论 0原文

只是警告，我一天前开始使用 R...如果有任何事情看起来非常简单，我深表歉意。

现在，我正在尝试让 R 接收包含撞击加速度计数据的 .txt 文件，并为其计算头部受伤标准测试。 HIC 测试要求将数据曲线按一定间隔进行积分。

该方程位于下面的链接...我尝试将其作为图像插入此处，但它不允许我这样做。显然我需要一些声誉点才能让我这样做。

equation

a(t) 是加速度曲线。

到目前为止，我在 R 中生成合适的曲线来匹配数据还没有遇到问题。 loess 函数运行得很好，这正是我一直在寻找的东西......我只是不知道如何集成它。据我所知，黄土是非参数回归，因此无法确定曲线的方程。有没有办法整合它？

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

任何帮助或有见地的评论将非常感激。

提前致谢，

韦斯

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

原文

分享到QQ

分享到微博

如果你对这篇内容有疑问，欢迎到本站社区发帖提问参与讨论，获取更多帮助，或者扫码二维码加入 Web 技术交流群。

发布评论

需要登录才能够评论，你可以免费注册一个本站的账号。

半暖夏伤 2024-12-17 21:33:33

您可以在 loess 模型上使用 predict 函数来创建与 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

You can use the predict function on your loess model to create a function to use with 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

And to extract the value alone: