在R中集成非参数曲线

时间:2011-10-19 12:47:53

标签: r

只是一个警告,我一天前开始使用R ......如果有什么似乎是天真的简单,我道歉。

现在我试图让R接收带有影响的acelerometer数据的.txt文件并计算头部伤害标准测试。 HIC测试要求来自数据的曲线在特定间隔上进行积分。

等式位于下面的链接...我试图将其作为图像插入此处,但它不会让我。显然我需要一些声望点,才能让我这样做。

equation

a(t)是加速曲线。

到目前为止,我还没有在R中生成合适的曲线以匹配数据的问题。黄土功能工作得很好,正是我正在寻找的东西...我只是不知道如何整合它。据我所知,黄土是一个非参数回归,因此没有办法确定曲线的方程。有没有办法整合它呢?

如果没有,是否有另一种方法可以使用不同的功能来完成此任务?

非常感谢任何帮助或有见识的评论。

提前致谢,

韦斯

詹姆斯还有一个问题,在使用integrate()函数时,如何才能获得没有文本和错误的数字?

1 个答案:

答案 0 :(得分:7)

您可以使用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