123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103 |
- % function ContratAnalysis:
- % this function creates contrast of SP vs fix. It uses te SPM.mat file within the provided directory
- %
- % input:
- % directory - directory to search for SPM.mat file
- function ContrastAnalysis(directory)
- % create batch file
- matlabbatch = [];
- matlabbatch{1}.spm.stats.con.spmmat = cellstr([directory '/SPM.mat']);
- matlabbatch{1}.spm.stats.con.consess{1}.tcon.name = 'SP > sacc';
- contrast = zeros(8,24);
- contrast(:,4) = 1;
- contrast(:,16) = -1;
- contrast = contrast';
- matlabbatch{1}.spm.stats.con.consess{1}.tcon.convec = contrast(:);
- % equivalent to
- % matlabbatch{1}.spm.stats.con.consess{1}.tcon.convec = [0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 -1 0 0 0 0 0 0 0 0 ...
- % 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 -1 0 0 0 0 0 0 0 0 ... % session 2
- % 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 -1 0 0 0 0 0 0 0 0 ... % session 3
- % 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 -1 0 0 0 0 0 0 0 0 ... % session 4
- % 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 -1 0 0 0 0 0 0 0 0 ... % session 5
- % 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 -1 0 0 0 0 0 0 0 0 ... % session 6
- % 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 -1 0 0 0 0 0 0 0 0 ... % session 7
- % 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 -1 0 0 0 0 0 0 0 0 ... % session 8
- % ];
- matlabbatch{1}.spm.stats.con.consess{1}.tcon.sessrep = 'none';
- matlabbatch{1}.spm.stats.con.consess{2}.tcon.name = 'sacc > SP';
- contrast = zeros(8,24);
- contrast(:,4) = -1;
- contrast(:,16) = 1;
- contrast = contrast';
- matlabbatch{1}.spm.stats.con.consess{2}.tcon.convec = contrast(:);
- matlabbatch{1}.spm.stats.con.consess{2}.tcon.sessrep = 'none';
- matlabbatch{1}.spm.stats.con.consess{3}.tcon.name = 'motion > sacc';
- contrast = zeros(8,24);
- contrast(:,10) = 1;
- contrast(:,16) = -1;
- contrast = contrast';
- matlabbatch{1}.spm.stats.con.consess{3}.tcon.convec = contrast(:);
- matlabbatch{1}.spm.stats.con.consess{3}.tcon.sessrep = 'none';
- matlabbatch{1}.spm.stats.con.consess{4}.tcon.name = 'sacc > motion';
- contrast = zeros(8,24);
- contrast(:,10) = -1;
- contrast(:,16) = 1;
- contrast = contrast';
- matlabbatch{1}.spm.stats.con.consess{4}.tcon.convec = contrast(:);
- matlabbatch{1}.spm.stats.con.consess{4}.tcon.sessrep = 'none';
- matlabbatch{1}.spm.stats.con.consess{5}.tcon.name = 'sp > motion';
- contrast = zeros(8,24);
- contrast(:,4) = 1;
- contrast(:,10) = -1;
- contrast = contrast';
- matlabbatch{1}.spm.stats.con.consess{5}.tcon.convec = contrast(:);
- matlabbatch{1}.spm.stats.con.consess{5}.tcon.sessrep = 'none';
- matlabbatch{1}.spm.stats.con.consess{6}.tcon.name = 'motion > sp';
- contrast = zeros(8,24);
- contrast(:,4) = -1;
- contrast(:,10) = 1;
- contrast = contrast';
- matlabbatch{1}.spm.stats.con.consess{6}.tcon.convec = contrast(:);
- matlabbatch{1}.spm.stats.con.consess{6}.tcon.sessrep = 'none';
- matlabbatch{1}.spm.stats.con.consess{7}.tcon.name = 'sp + motion > sacc';
- contrast = zeros(8,24);
- contrast(:,4) = 1;
- contrast(:,10) = 1;
- contrast(:,16) = -2;
- contrast = contrast';
- matlabbatch{1}.spm.stats.con.consess{7}.tcon.convec = contrast(:);
- matlabbatch{1}.spm.stats.con.consess{7}.tcon.sessrep = 'none';
- matlabbatch{1}.spm.stats.con.consess{8}.tcon.name = 'SP > baseline';
- contrast = zeros(8,24);
- contrast(:,4) = 1;
- contrast = contrast';
- matlabbatch{1}.spm.stats.con.consess{8}.tcon.convec = contrast(:);
- matlabbatch{1}.spm.stats.con.consess{8}.tcon.sessrep = 'none';
- matlabbatch{1}.spm.stats.con.consess{9}.tcon.name = 'sacc > baseline';
- contrast = zeros(8,24);
- contrast(:,16) = 1;
- contrast = contrast';
- matlabbatch{1}.spm.stats.con.consess{9}.tcon.convec = contrast(:);
- matlabbatch{1}.spm.stats.con.consess{9}.tcon.sessrep = 'none';
- matlabbatch{1}.spm.stats.con.consess{10}.tcon.name = 'motion > baseline';
- contrast = zeros(8,24);
- contrast(:,10) = 1;
- contrast = contrast';
- matlabbatch{1}.spm.stats.con.consess{10}.tcon.convec = contrast(:);
- matlabbatch{1}.spm.stats.con.consess{10}.tcon.sessrep = 'none';
- % execute job
- spm_jobman('run',matlabbatch);
- end
|