Please access the Jupyter notebook at: https://drive.google.com/file/d/1cYMEdBIuAPpO-5OQoFhP34RvgsZ-cMre/view?usp=drive_link.
This will open a tab within Google Drive; to run the Jupyter notebook in Google Colab, please save a copy of the .ipynb
file to your Google Drive and open it.
This should automatica
5809
lly open the file in a Google Colab instance, which will allow you to compile the code and view the outputs directly from your browser without installing anything on your computer.
Today we will be exploring time-series analysis with two comprehensive libraries in python:
We'll be working with neuroimaging time-series from a case study described in our original pyspi publication. Briefly, this entails resting-state functional magnetic resonance imaging (fMRI) data, which approximates neural activity over time based on hemodynamic coupling. Typically, researchers will extract time series from parcellated cortical and subcortical brain regions, and will then focus either on local dynamics or pairwise coupling:
From these time series, we can then summarise properties of local dynamics with hctsa and properties of pairwise coupling between brain regions using pyspi:
In the original pyspi publication, we asked how well each of 237 statistics for pairwise interactions (SPIs) in pyspi could distinguish fMRI scans obtained while participants were at rest (i.e., chilling in the scanner not actively performing a task) vs. while viewing a film.
We found several SPIs that made this distinction with statistically significant accuracy, many of which surpassed the more typically used Pearson correlation coefficient: