Compute jerk from position data

Computes jerk (third derivative of position) from position assays.

Usage

computeJerk(pe, method = "central", sampling_rate = NULL)

Arguments

pe

A PhysioExperiment object containing position assays.

method

Finite difference method: "central" (default), "forward", or "backward".

sampling_rate

Sampling rate in Hz. If NULL (default), uses samplingRate(pe).

Value

PhysioExperiment with new jerk_x, jerk_y, jerk_z assays added.

References

Winter DA (2009). "Biomechanics and Motor Control of Human Movement." 4th ed. John Wiley & Sons.

See also

computeAcceleration() for second-order derivatives, computeVelocity() for first-order derivatives, differentiate() for general-purpose numerical differentiation.

Examples

demo <- demoMoCapData(seed = 1)
pe_jerk <- computeJerk(demo$mocap)
SummarizedExperiment::assayNames(pe_jerk)
#> [1] "position_x" "position_y" "position_z" "jerk_x"     "jerk_y"    
#> [6] "jerk_z"