Preprocessing Functions for PhysioExperiment Objects • PhysioPreprocess

Preprocessing Functions for PhysioExperiment Objects

PhysioPreprocess provides a comprehensive preprocessing toolkit for multi-modal physiological signal data. Built on top of PhysioCore, it delivers 30 exported functions covering digital filtering (Butterworth, FIR, notch), resampling and interpolation, ICA-based artifact removal, bad channel detection, baseline correction, re-referencing, and reproducible preprocessing pipelines – all operating directly on PhysioExperiment objects.

Installation

You can install PhysioPreprocess from r-universe:

install.packages("PhysioPreprocess",
  repos = c("https://x-biosignal.r-universe.dev", "https://cloud.r-project.org"))

Or install the development version from GitHub:

# install.packages("remotes")
remotes::install_github("x-biosignal/PhysioPreprocess")

Quick Start

library(PhysioPreprocess)

# Create sample EEG data (4 seconds, 4 channels, 250 Hz)
signal_matrix <- matrix(rnorm(1000 * 4), nrow = 1000, ncol = 4)
colnames(signal_matrix) <- c("Fz", "Cz", "Pz", "Oz")

pe <- PhysioExperiment(
  assays = list(raw = signal_matrix),
  samplingRate = 250
)

# Apply a 1-40 Hz bandpass Butterworth filter
pe <- filterSignals(pe, lowcut = 1, highcut = 40, order = 4)

# Remove 50 Hz line noise with a notch filter
pe <- notchFilter(pe, freq = 50)

# ICA-based artifact removal
ica_result <- runICA(pe, n_components = 4)
pe <- removeICAComponents(pe, ica_result, components = c(1))

# Re-reference to average
pe <- rereference(pe, ref = "average")

# Check reference status
isAverageReferenced(pe)  # TRUE

Features

Digital Filters

A full suite of frequency-domain filters for physiological signals:

Butterworth: filterSignals(), butterworthFilter() – bandpass, lowpass, and highpass IIR filters with configurable order
FIR: firFilter() – finite impulse response filters for linear-phase filtering
Notch: notchFilter() – remove power line interference (50/60 Hz) and harmonics
Detrending: detrendSignal(), detrendSignals() – remove linear or polynomial trends from signals

Resampling

Flexible sample rate conversion and multi-rate signal support:

Rate conversion: resample(), decimate(), interpolate() – change sampling rates with anti-aliasing
Multi-rate support: assaySamplingRates(), setAssaySamplingRate() – manage per-assay sampling rates for mixed-rate recordings

Artifact Handling

Automated and semi-automated artifact detection and removal:

ICA decomposition: icaDecompose(), runICA() – decompose signals into independent components for artifact identification
ICA removal: icaRemove(), removeICAComponents() – remove artifact components (e.g., eye blinks, muscle activity) and reconstruct clean signals
Bad channels: detectBadChannels(), interpolateBadChannels() – identify noisy or flat channels and interpolate from neighbors
Epoch rejection: rejectBadEpochs() – reject epochs exceeding amplitude or variance thresholds

Baseline Correction

Remove baseline drift and DC offsets:

baselineCorrect() – subtract mean of a specified baseline window from each epoch
removeBaseline() – flexible baseline removal with configurable time windows and methods

Re-referencing

Electrode re-referencing for EEG and related modalities:

rereference() – re-reference to average, specific electrode(s), or linked mastoids
getCurrentReference() – query the current reference scheme
isAverageReferenced() – check whether average reference has been applied

Pipeline API

Build reproducible, shareable preprocessing chains:

createPipeline() – define an ordered sequence of preprocessing steps with parameters
applyPipeline() – apply a saved pipeline to new data for consistent processing across datasets

# Build a reusable preprocessing pipeline
pipeline <- createPipeline(
  list(filterSignals, lowcut = 1, highcut = 40),
  list(notchFilter, freq = 50),
  list(rereference, ref = "average")
)

# Apply to any PhysioExperiment
pe_clean <- applyPipeline(pe, pipeline)

Dependencies

R (>= 4.2)
PhysioCore – core data structures
signal – DSP primitives
methods, SummarizedExperiment, S4Vectors, stats
Suggests: fastICA, testthat, knitr, rmarkdown

PhysioExperiment Ecosystem

PhysioPreprocess is part of the PhysioExperiment ecosystem, a suite of R packages for multi-modal physiological signal analysis:

Package	Description
PhysioCore	Core data structures and accessors
PhysioIO	File I/O (EDF, HDF5, BIDS, CSV, MAT)
PhysioPreprocess	Signal preprocessing and artifact removal
PhysioAnalysis	Spectral analysis, epoching, statistics, visualization
PhysioMoCap	Motion capture data processing
PhysioOpenSim	OpenSim biomechanical modeling integration

Visit the r-universe page to browse all available packages.

License

MIT License. See LICENSE for details.

Author

Yusuke Matsui

Governance & support

Part of the Physio ecosystem. Community and policy documents (Code of Conduct, Contributing, Governance, Support, Security, and Deprecation & lifecycle) live in the umbrella repository.