This repository contains the Python code to simulate an fNIRS recording session for a block design, taking Gervain et al. (2012)'s design as an example.
Python >= 3.12 (earlier versions may work as well), and Python packages NumpPy, Matplotlib, and SciPy. To install them execute this on your console/terminal:
python -m pip install numpy matplotlib scipyimport numpy as np
import matplotlib.pyplot as plt
from src import HRF
sfreq = 10
time = np.arange(0, 20, 1/sfreq)
hbo, hbr = hrf(time, ratio=0.25)
plt.plot(hbo)
plt.plot(hbr)See this Jupyter notebook for a complete walkthrough.
Glover, G. H. (1999). Deconvolution of impulse response in event-related BOLD fMRI1. Neuroimage, 9(4), 416-429. https://doi.org/10.1006/nimg.1998.0419
Gervain, J., Berent, I., & Werker, J. F. (2012). Binding at birth: The newborn brain detects identity relations and sequential position in speech. Journal of Cognitive Neuroscience, 24(3), 564-574. https://doi.org/10.1162/jocn_a_00157