Neuro-fPCA-fMRI 0.0.6

Functional PCA on fMRI data to extract dominant temporal components.

Full Documentation

• Github Repository: https://github.com/yanivassf-lab/fMRI

• Documentation: https://fmri.readthedocs.io

• PyPI Package: https://pypi.org/project/Neuro-fPCA-fMRI

Overview

This module implements a functional principal component analysis (fPCA) pipeline designed for fMRI data. The overall procedure involves the following steps:

1. Data Loading and Preprocessing

   The 4D fMRI image and an associated 3D binary mask are loaded. The mask extracts voxel-specific time series, and the time axis is defined based on the number of timepoints.

2. B-Spline Basis Construction

   A set of B-spline basis functions is generated to represent smooth temporal dynamics. Each voxel’s time series is approximated as a linear combination of these basis functions.

3. Penalty Matrix Construction

   A penalty matrix is computed to ensure smoothness of the estimated coefficients. This penalty is defined via an integral involving products of derivatives (typically the second derivative) of the basis functions.

4. Regularized Spline Regression

   Each voxel time series is approximated by a spline expansion. The coefficients are estimated via a regularized regression problem, and a generalized cross-validation (GCV) criterion selects a voxel-specific regularization parameter (\(\lambda\)).

5. Functional Principal Component Analysis

   PCA is performed on the coefficient matrix. Multiplying eigenfunctions by coefficients yields voxel-specific importance scores mapped back to brain images, while the basis applied to each eigenfunction gives the temporal dynamics curve shared across voxels.

Output

For each eigenfunction, the following figures are generated:

• Voxel importance maps

  

• Intensity plot

  

• Signal intensity of best-scoring voxel

  

Mathematical background

This program implements the methodology from:

Roberto Viviani, Georg Grön and Manfred Spitzer. Functional Principal Component Analysis of fMRI Data.

Full mathematical derivation is available here: Detailed Math Docs

Credits

This package was created using Cookiecutter and the audreyr/cookiecutter-pypackage template.

License

• Free software: GNU GPLv3

Author

Refael Kohen <[e-mail]>, Yaniv Assaf Lab, Tel Aviv University.