spm-批量数据统计建模

吡哆

已于 2024-03-17 18:44:31 修改

阅读量93

点赞数 2

文章标签： matlab

于 2024-03-17 18:37:52 首次发布

本文链接：https://blog.csdn.net/2201_75816423/article/details/136787036

版权

该代码片段展示了如何在Matlab中使用SPM进行功能性连接（FC）数据的三组方差分析，分别针对NSA、NSI和HC脑区。通过批量处理，提升了分析效率，适用于大规模可循环的数据集。

摘要由CSDN通过智能技术生成

clc
clear
raw_path='D:\2023_12_NSSI_YLL\2024_01_15\results_nsa_nsi_hc\fmri\FC_FunImgARWDCF\'
for j=3:44
    roi_num=['zROI',num2str(j),'FCMap_']
    output_dir=[raw_path,roi_num,'\result_f']
    nsa_dir=dir([raw_path,roi_num,'\nsa33\','*.nii'])
    nsi_dir=dir([raw_path,roi_num,'\nsi23\','*.nii'])
    hc_dir=dir([raw_path,roi_num,'\hc23\','*.nii'])
    for i=1:length(nsa_dir)
        nsa_data{i,1}=[raw_path,roi_num,'\nsa33\',nsa_dir(i).name];
    end
    for i=1:length(nsi_dir)
        nsi_data{i,1}=[raw_path,roi_num,'\nsi23\',nsi_dir(i).name];
    end
    for i=1:length(hc_dir)
        hc_data{i,1}=[raw_path,roi_num,'\hc23\',hc_dir(i).name];
    end
    matlabbatch{1}.spm.stats.factorial_design.dir = {output_dir};
%%
matlabbatch{1}.spm.stats.factorial_design.des.anova.icell(1).scans = nsa_data;
%%
%%
matlabbatch{1}.spm.stats.factorial_design.des.anova.icell(2).scans =nsi_data;
%%
%%
matlabbatch{1}.spm.stats.factorial_design.des.anova.icell(3).scans =hc_data;
%%
matlabbatch{1}.spm.stats.factorial_design.des.anova.dept = 0;
matlabbatch{1}.spm.stats.factorial_design.des.anova.variance = 1;
matlabbatch{1}.spm.stats.factorial_design.des.anova.gmsca = 0;
matlabbatch{1}.spm.stats.factorial_design.des.anova.ancova = 0;
matlabbatch{1}.spm.stats.factorial_design.cov = struct('c', {}, 'cname', {}, 'iCFI', {}, 'iCC', {});
matlabbatch{1}.spm.stats.factorial_design.multi_cov = struct('files', {}, 'iCFI', {}, 'iCC', {});
matlabbatch{1}.spm.stats.factorial_design.masking.tm.tm_none = 1;
matlabbatch{1}.spm.stats.factorial_design.masking.im = 1;
matlabbatch{1}.spm.stats.factorial_design.masking.em = {'D:\2023_12_NSSI_YLL\2024_01_15\RSFC_2024_02_28\Masks\AllResampled_GreyMask_02_91x109x91.nii,1'};
matlabbatch{1}.spm.stats.factorial_design.globalc.g_omit = 1;
matlabbatch{1}.spm.stats.factorial_design.globalm.gmsca.gmsca_no = 1;
matlabbatch{1}.spm.stats.factorial_design.globalm.glonorm = 1;
matlabbatch{2}.spm.stats.fmri_est.spmmat(1) = cfg_dep('Factorial design specification: SPM.mat File', substruct('.','val', '{}',{1}, '.','val', '{}',{1}, '.','val', '{}',{1}), substruct('.','spmmat'));
matlabbatch{2}.spm.stats.fmri_est.write_residuals = 0;
matlabbatch{2}.spm.stats.fmri_est.method.Classical = 1;
matlabbatch{3}.spm.stats.con.spmmat(1) = cfg_dep('Model estimation: SPM.mat File', substruct('.','val', '{}',{2}, '.','val', '{}',{1}, '.','val', '{}',{1}), substruct('.','spmmat'));
matlabbatch{3}.spm.stats.con.consess{1}.fcon.name = 'nsa-nsi-hc';
matlabbatch{3}.spm.stats.con.consess{1}.fcon.weights = [1 -1 0
                                                        1 0 -1];
matlabbatch{3}.spm.stats.con.consess{1}.fcon.sessrep = 'none';
matlabbatch{3}.spm.stats.con.delete = 0;
spm_jobman('run',matlabbatch)
end

由于功能连接基于先验生成了很多种子点，为了节约时间，因此批量用spm做三组方差分析，同理也可以应用于其他的大量可循环的统计