簡單來說,分為兩大塊。第一,PPI準(zhǔn)備活動(dòng);第二,正式的PPI操作。
【下面的例子是以我的實(shí)驗(yàn)數(shù)據(jù)參數(shù)為基準(zhǔn)的,TR=2s,總共5個(gè)run,每個(gè)run包含186個(gè)volume,并且默認(rèn)已經(jīng)經(jīng)過基礎(chǔ)的GLM建?!?/p>
PPI準(zhǔn)備活動(dòng)
準(zhǔn)備好預(yù)處理smooth后的功能像
合并頭動(dòng)
如果有5個(gè)run,那么需要將頭動(dòng)信息合并成一個(gè)run的。新的頭動(dòng)參數(shù)中第一個(gè)run的數(shù)據(jù)不變,第二個(gè)run的每一行數(shù)據(jù)加上第一個(gè)run的最后一行,第三個(gè)run的每一行數(shù)據(jù)加上新弄好的第二個(gè)run的最后一行,以此類推。合并時(shí)間序列
類似于上面的合并頭動(dòng),這一步的目的是將之前的每個(gè)run分別的timefile合并成一個(gè)大的timefile,因?yàn)樵赑PI操作中,只能輸入一個(gè)大run。我的timefile是以秒為單位的
具體操作:----> 第一個(gè)run的數(shù)據(jù)不變 ---->第二個(gè)run的所有數(shù)值加上186??2 ---->第三個(gè)run的所有數(shù)值加上372??2 ---->第4個(gè)run的所有數(shù)值加上558??2 ---->第5個(gè)run的所有數(shù)值加上744??2
正式的PPI操作
1.一個(gè)GLM
先做一個(gè)正常的GLM(Block regressor:1和0的矩陣,區(qū)分不同的run; 頭動(dòng)regressor放入頭動(dòng)參數(shù))。
一個(gè)PPI
點(diǎn)PPIs---->psychophysiologic interaction 生成一個(gè)變量PPI,然后用代碼循環(huán)生成不同名字的PPI。(生成對(duì)應(yīng)的PPI.mat文件)再做一個(gè)GLM(但是,不選conditions,不用填)
需要先load具體的PPI.mat,才能在后面value數(shù)值中直接寫PPI.ppi、PPI.Y等
Regressors: Regressor1, Name=PPI-interaction, value=PPI.ppi
Regressor2, Name=dmPFC_BOLD, value=PPI.Y
Regressor3, Name=Mental-Rnd, value=PPI.P
Regressor4, Name=Block1, value=[1,0]的矩陣
Regressor5,.......
....(回歸的Block參數(shù))
Multiple regressors:放入頭動(dòng)文件