只是一个警告,我一天前开始使用R ......如果有什么似乎是天真的简单,我道歉。
现在我试图让R接收带有影响的acelerometer数据的.txt文件并计算头部伤害标准测试。 HIC测试要求来自数据的曲线在特定间隔上进行积分。
等式位于下面的链接...我试图将其作为图像插入此处,但它不会让我。显然我需要一些声望点,才能让我这样做。
a(t)是加速曲线。
到目前为止,我还没有在R中生成合适的曲线以匹配数据的问题。黄土功能工作得很好,正是我正在寻找的东西...我只是不知道如何整合它。据我所知,黄土是一个非参数回归,因此没有办法确定曲线的方程。有没有办法整合它呢?
如果没有,是否有另一种方法可以使用不同的功能来完成此任务?
非常感谢任何帮助或有见识的评论。
提前致谢,
韦斯
詹姆斯还有一个问题,在使用integrate()函数时,如何才能获得没有文本和错误的数字?
答案 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