Summary
* The authors seek to develop an approach for modeling spatiotemporal brain dynamics, as measured with resting state functional MRI. For this, they have extended a transformer-based architecture to handle time-varying 3D scan data.
* They seek to predict subject traits and attributes like age, sex, and cognitive scores from the learned representations using three publicly available datasets.
* They have succeeded in outperforming state-of-the-art models in these tasks while reducing the computational complexity and increasing throughput in model evaluations.
* Additionally, they have demonstrated the usefulness of pre-training such models on downstream tasks.
* Lastly, the relationships captured by the model between the subject attributes and the brain regions are consistent with prior work.
Strengths
* The problem statement is well-motivated, and the limitations of prior works (i.e., ROI-based and two-step approaches) are detailed.
* The results were replicated across three datasets and various tasks, showcasing this architecture's robustness.
* The authors show the learned features are general enough to help downstream tasks through pre-training.
* The authors have performed computational complexity analysis. Because their model has almost double the throughput compared to the second-best-performing model, it paves the way to use these models in real-time settings.
* Regions of the brain used by the model (as in Section 4.4) to predict age align with previous work, which increases confidence that the model is learning meaningful features.
* The authors have submitted their full code, which assists replicability and helps the scientific community extend their work.
Weaknesses
* The novelty of the architecture vis-a-vis latest advances in transformer architectures ([23], [34], [35]) seems limited.
* The authors aimed to learn representations for brain dynamics but focussed only on predicting traits (age, sex etc), because they used only resting state scans. Dynamic attributes such as task-level performance measures -- e.g., reaction times and accuracies in HCP tasks), cognitive load in a working memory task etc -- are relevant to establish the generality of the findings. The model should be tested on these tasks also.
* While the authors have detailed the approach for the two contrastive losses if they could justify the choice of positive and negative samples through domain knowledge or established methods, it would help increase the confidence in the proposed approach.
* The authors have mentioned performing ablation studies to "substantiate their modeling choices." I expect that this includes studying the effect of hyperparameters like the number of layers, channel size, window size, etc., given that this is a novel architecture. Details of such studies should have been mentioned.
* Additionally, it would be helpful to see the loss curves of models for different datasets and targets. It will help understand the choice of a low number of epochs for training various models.
Questions
* Could the authors clarify if the parameters for z-scoring, i.e., mean and variance, were learned over the whole data or only on the training data (& applied to validation and test data)?
* I would expect that using the whole fMRI time series for any task prediction would provide better results than using sub-sequences. Could the authors provide more intuition about the seemingly counterintuitive results obtained in Section 4.6, especially in UKB (e.g., age prediction)?
* While we commend that the authors have shared their code, it would be helpful if they could also share their trained models (especially due to the high compute power required to generate these results). It would help the scientific community at large and also be in line with the spirit of Section 4.3.
Rating
6: Weak Accept: Technically solid, moderate-to-high impact paper, with no major concerns with respect to evaluation, resources, reproducibility, ethical considerations.
Confidence
5: You are absolutely certain about your assessment. You are very familiar with the related work and checked the math/other details carefully.
Limitations
* Explanation methods like IG are generally brittle. It would be helpful if one more explanation method could confirm the results in Section 4.4
* While a shifting window enables efficient feature extraction, it limits model expressivity. It would be useful to know how to extract features that span large spatial/temporal ranges. This could be particularly relevant for task-based fMRI data since signals of varying timescales could affect behavior in such paradigms.
* Authors have mentioned that they divided the fMRI data into sub-sequences due to memory constraints. Moreover, as the authors have used four Nvidia A100s (a total of 160GB of GPU memory with NvLink support), they are already on the higher end of compute power. It would be helpful to understand what those memory constraints were and suggestions, if any, on tackling them, as this will have implications in replicating and/or extending their work.