library(lookr)
library(robfilter)
gaze <- Gazedata("inst/tests/data/Coartic_WFFArea_1a/001P00XS1/Coartic_Block1_001P00XS1.gazedata")
gaze4 <- gaze[gaze$TrialNo == 4, ]
# Filtering window of approximately 150 ms
smoothing_width <- floor(150 / lwl_constants$ms_per_frame)
xmean <- med.filter(gaze4$XMean, width = smoothing_width)
#> 27 out of 467 time series values in y are missing.
# Plot raw x-gazes and median filtered values
plot(xmean)
How nice would this be...