SPM/連線
一個會話(也稱為“執行”)是在掃描器中收集資料的一段時間。在會話之間,資料採集會暫停,通常是為了讓參與者休息並能夠與實驗者交流。SPM 預設情況下會為單變數推斷指定多會話 GLM。但是,您可能希望將所有會話連線成一個,以便進行後續(特別是連線性)分析。
對於基於 GLM 的單變數推斷,連線不是必需的。此外,它可能會導致與多會話“完整模型”相比產生不同的估計值,因為 1) 如果一個條件太靠近會話結束,在與 HRF 卷積後它可能會錯誤地對來自後續會話的資料進行建模(在每個會話結束之前獲取額外的掃描可以有效地減輕這種影響)和 2) 它透過在整個連線的會話中為每個條件提供一個 beta 來停用條件 × 會話互動。
然而,SPM 時間序列提取 是基於每個會話進行的。因此,通常需要將多個會話連線起來並從中提取時間序列,以進行連線性分析,例如 PPI 或 DCM。
連線中最棘手的部分是確保您的發病時間(跨會話整理)是正確的。如果您已在每個會話中獲取了發病時間,為了將這些發病時間附加到新的連線的 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];
其他迴歸量(例如頭部運動迴歸量)應以類似的方式連線。
然後您可以指定您的單會話一級設計矩陣。您現在應該不再指定多個會話,而應該只指定一個會話,其中包含來自所有會話的影像體積以及附加的發病時間和其他迴歸量。執行批處理以指定模型,但現在不要進行估計。
在進行估計之前,您仍然需要告訴 SPM 哪些掃描屬於哪個原始會話,以便它可以調整其影響。這將為您的單會話 GLM 新增塊效應迴歸量(替換設計矩陣中通常的平均值列),並更正高通濾波器和時間非球面度計算以說明原始會話長度。
在主 Matlab 視窗中執行以下程式碼。同樣,這裡的“掃描”表示原始時間序列中每個會話中的體積數量。
scans = [60 55 40];
spm_fmri_concatenate('SPM.mat', scans);
調整後的 GLM (SPM.mat) 將替換您之前指定的 GLM(它的副本將儲存為 SPM_backup.mat)。現在您可以按照正常方式估計 GLM 並新增對比(最重要的是,用於時間序列提取的感興趣效應 F 對比)。
現在您可以使用 SPM 感興趣區域實用程式來提取單會話的時間序列。