Ensemble average over fixed-length segments

Averages a (n_segments, n_samples) matrix down its rows. Standard synchronous-averaging recipe for repeated stimulus responses (e.g. evoked potentials).

Usage

morie_dsp_ensemble_average(segments)

Arguments

segments

Numeric matrix (rows = trials, cols = samples).

Value

Row mean as a length-ncol vector.

References

Rangayyan & Krishnan (2015), Ch. 3, sec. 3.3.

Examples

set.seed(1)
t <- seq.int(0, 99) / 100
segs <- t(replicate(8, sin(2 * pi * 5 * t) + rnorm(100, sd = 0.3)))
y <- morie_dsp_ensemble_average(segs)
sd(y) < sd(segs[1, ])  # averaging reduces noise
#> [1] TRUE