Question

Just a warning, I started using R a day ago…my apologies if anything seems idiotically simple.

Right now im trying to have R take in a .txt file with acelerometer data of an impact and calculate a Head injury criterion test for it. The HIC test requires that curve from the data be integrated on a certain interval.

The equation is at the link below…i tried to insert it here as an image but it would not let me. Apparently i need some reputation points before it’ll let me do that.

a(t) is the aceleration curve.

So far i have not had an issue generating a suitable curve in R to match the data. The loess function worked quite well, and is exactly the kind of thing i was looking for…i just have no idea how to integrate it. As far as i can tell, loess is a non-parametric regression so there is no way to determine the equation of the curve iteslf. Is there a way to integrate it though?

If not, is there another way to acomplish this task using a different function?

Any help or insighful comments would be very much appreciated.

Thanks in advance,

Wes

One more question though James, how can i just get the number without the text and error when using the integrate() function?

Answer 1

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:

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