2025-11-11T16:04:09.679573

A Mixed Model Approach for Estimating Regional Functional Connectivity from Voxel-level BOLD Signals

Liu, Zhang, Tran et al.

Resting-state brain functional connectivity quantifies the synchrony between activity patterns of different brain regions. In functional magnetic resonance imaging (fMRI), each region comprises a set of spatially contiguous voxels at which blood-oxygen-level-dependent signals are acquired. The ubiquitous Correlation of Averages (CA) estimator, and other similar metrics, are computed from spatially aggregated signals within each region, and remain the quantifications of inter-regional connectivity most used by neuroscientists despite their bias that stems from intra-regional correlation and measurement error. We leverage the framework of linear mixed-effects models to isolate different sources of variability in the voxel-level signals, including both inter-regional and intra-regional correlation and measurement error. A novel computational pipeline, focused on subject-level inter-regional correlation parameters of interest, is developed to address the challenges of applying maximum (or restricted maximum) likelihood estimation to such structured, high-dimensional spatiotemporal data. Simulation results demonstrate the reliability of correlation estimates and their large sample standard error approximations, and their superiority relative to CA. The proposed method is applied to two public fMRI data sets. First, we analyze scans of a dead rat to assess false positive performance when connectivity is absent. Second, individual human brain networks are constructed for subjects from a Human Connectome Project test-retest database. Concordance between inter-regional correlation estimates for test-retest scans of the same subject are shown to be higher for the proposed method relative to CA.

academic

A Mixed Model Approach for Estimating Regional Functional Connectivity from Voxel-level BOLD Signals

基本信息

论文ID: 2211.02192
标题: A Mixed Model Approach for Estimating Regional Functional Connectivity from Voxel-level BOLD Signals
作者: Ruobin Liu, Chao Zhang, Chau Tran, Sophie Achard, Wendy Meiring, Alexander Petersen
分类: stat.ME (Statistics - Methodology), stat.AP (Statistics - Applications)
发表时间: 2022年11月（arXiv预印本，2025年11月更新）
论文链接: https://arxiv.org/abs/2211.02192

摘要

本文针对功能磁共振成像(fMRI)中静息态脑功能连接估计问题，提出了一种基于线性混合效应模型的新方法。传统的"平均相关性"(Correlation of Averages, CA)估计器虽然被广泛使用，但由于区域内相关性和测量误差导致的偏差问题一直存在。本文通过线性混合效应模型框架，分离体素级信号中的不同变异性来源，包括区域间和区域内相关性以及测量误差。研究开发了一个新颖的计算管道，专注于个体水平的区域间相关参数估计，并采用最大似然估计方法处理高维时空数据的挑战。仿真结果证明了相关性估计的可靠性及其相对于CA的优越性。

研究背景与动机

问题定义

核心问题: 如何从体素级BOLD信号准确估计脑区间功能连接性，避免传统方法的偏差问题
技术挑战:
- fMRI数据具有复杂的时空依赖性
- 体素数量远大于时间维度，导致计算挑战
- 区域内空间相关性和测量误差影响连接性估计的准确性

研究重要性

功能连接性是研究神经退行性疾病和意识障碍等病理学的基础
准确的连接性估计对个体和群体水平的神经科学研究至关重要
现有方法的偏差可能影响疾病诊断和个体特征化研究

现有方法局限性

传统CA估计器存在的问题：

偏差问题: 由区域内相关性αⱼ和噪声信号比βⱼ约束，导致估计值偏向0
依赖性忽略: 未考虑体素级信号的时空依赖性
参数约束: 连接性参数受采样方案和机器噪声影响，缺乏内在性

核心贡献

提出新的统计模型: 基于线性混合效应模型的体素级BOLD信号建模框架，明确区分区域间和区域内变异性
开发高效估计方法:
- 两阶段估计策略，结合限制最大似然估计(ReML)
- 首次在功能连接建模中应用Vecchia似然近似方法
理论保证: 提供估计器的大样本性质和渐近推断理论
实证验证: 在仿真和真实数据(死鼠扫描、HCP测试-重测数据)上验证方法的优越性

X = Zμ + Uη + γ + ε

其中：

μⱼ: 区域固定效应(区域均值)
ηⱼₘ: 区域随机效应(诱导区域间依赖性)
γⱼₗₘ: 体素级随机效应(诱导区域内依赖性)
εⱼₗₘ: 测量误差

协方差结构参数化

区域间相关性: Var(η) = (SRS) ⊗ A，其中R = {ρⱼⱼ'}为目标相关矩阵
区域内结构: Λⱼ = Cⱼ ⊗ Bⱼ (可分离的空间-时间协方差)
核函数定义:
- 空间核: Matérn核 K(d; ν, φ)
- 时间核: 高斯核 H(|m-m'|; τ)

技术创新点

两阶段估计策略

阶段1: 区域特异性参数估计

使用ReML估计每个区域的参数 θⱼ = kᵧⱼ, σ²ᵧⱼ, φᵧⱼ, τᵧⱼ
通过限制似然消除区域效应的影响

阶段2: 全局和区域间参数估计

估计区域间相关性参数 θ = τη, kη, ρ₁₂, σ²η
固定阶段1的估计值，专注于连接性参数

Vecchia近似

为解决计算复杂度问题(O(N³)时间，O(N²)内存)，采用Vecchia似然近似：

p(X) ≈ p(X_π(1)) ∏ᵢ₌₂ᴺ p(X_π(i) | X_π(j), j ∈ Jᵢ)

通过小的条件集|Jᵢ|=100实现计算效率提升

实验设置

数据集

仿真数据:
- J=3个脑区，M=60个小波系数
- 使用活鼠实验的空间坐标(L₁=41, L₂=25, L₃=77个体素)
- 变化信号强度δⱼ ∈ {0.1, 0.5, 0.7}和空间协方差ψⱼ ∈ {0.2, 0.5, 0.8}
真实数据:
- 死鼠扫描数据(验证假阳性率)
- HCP测试-重测数据库(42个被试，J=92个默认模式网络区域)