SPM/连接

一个会话（也称为“运行”）是在扫描仪中收集数据的一段时间。在会话之间，数据采集会暂停，通常是为了让参与者休息并能够与实验者交流。SPM 默认情况下会为单变量推断指定多会话 GLM。但是，您可能希望将所有会话连接成一个，以便进行后续（特别是连接性）分析。

为什么要连接？

对于基于 GLM 的单变量推断，连接不是必需的。此外，它可能会导致与多会话“完整模型”相比产生不同的估计值，因为 1) 如果一个条件太靠近会话结束，在与 HRF 卷积后它可能会错误地对来自后续会话的数据进行建模（在每个会话结束之前获取额外的扫描可以有效地减轻这种影响）和 2) 它通过在整个连接的会话中为每个条件提供一个 beta 来禁用条件 × 会话交互。

然而，SPM 时间序列提取是基于每个会话进行的。因此，通常需要将多个会话连接起来并从中提取时间序列，以进行连接性分析，例如 PPI 或 DCM。

步骤 (fMRI)

准备

连接中最棘手的部分是确保您的发病时间（跨会话整理）是正确的。如果您已在每个会话中获取了发病时间，为了将这些发病时间附加到新的连接的 GLM，您需要在这些发病时间中添加之前会话的会话长度。例如，假设您有 3 个会话以及相应的发病时间 1-3，每个会话分别有 60、55、44 个体积扫描

scans = [60 55 40];

% If onsets are in scans
onsets_concat = [onsets1,onsets2 + 60, onsets3 + 115];

% If onsets are in seconds
TR = 2; % Replace with the actual value according to your scanning protocol
onsets_concat = [onsets1,onsets2 + 60*TR, onsets3 + 115*TR];

其他回归量（例如头部运动回归量）应以类似的方式连接。

指定连接的 GLM

然后您可以指定您的单会话一级设计矩阵。您现在应该不再指定多个会话，而应该只指定一个会话，其中包含来自所有会话的图像体积以及附加的发病时间和其他回归量。运行批处理以指定模型，但现在不要进行估计。

按会话调整和估计

在进行估计之前，您仍然需要告诉 SPM 哪些扫描属于哪个原始会话，以便它可以调整其影响。这将为您的单会话 GLM 添加块效应回归量（替换设计矩阵中通常的平均值列），并更正高通滤波器和时间非球面度计算以说明原始会话长度。

在主 Matlab 窗口中运行以下代码。同样，这里的“扫描”表示原始时间序列中每个会话中的体积数量。

scans = [60 55 40];
spm_fmri_concatenate('SPM.mat', scans);

调整后的 GLM (SPM.mat) 将替换您之前指定的 GLM（它的副本将保存为 SPM_backup.mat）。现在您可以按照正常方式估计 GLM 并添加对比（最重要的是，用于时间序列提取的感兴趣效应 F 对比）。

现在您可以使用 SPM 感兴趣区域实用程序来提取单会话的时间序列。