Perform event detection as in Betzel et al (2021): Individualized event structure drives individual differences in whole-brain functional connectivity

• Example parcellated resting fMRI BOLD time series from here.
• System labels for each brain region.
• Analysis script.

1. Reads in the parcel time series data, generates edge time series, and calculates the root sum squared (RSS) amplitude at each frame.
2. Uses a circular-shift null model to generate null time series and repeats the procedure from Step 1, resulting in a null distribution of RSS values.
3. Identifies statistically significant frames (observed RSS > than that of null) and partitions time series into segments of temporally contiguous supra-threshold frames.
4. Identifies peak RSS within each segment and extracts pattern of activity (parcel time series) and edge time series (co-fluctuation) at that instant.
5. Makes a couple figures.

If you use this code, please cite: R Betzel, S Cutts, S Greenwell, O Sporns (2021). Individualized event structure drives individual differences in whole-brain functional connectivity. bioRxiv link to paper