Hi GPPI experts,

I have a question about the contrast I should use for the P.contrast field. I have 3 task conditions and 2 parametric modulators for one of the three conditions, and I am interested in the effects of these 2 parametric modulators. I've defined the 3 task conditions (see script below) but was wondering for the P.contrast field, should I create an effects of interest contrast which includes my 3 task conditions and the parametric modulators (i.e. 1 0 0 0 0; 0 1 0 0 0; 0 0 1 0 0; 0 0 0 1 0; 0 0 0 0 1)?

And if I exclude / have 0s in my F contrast for the motion regressors, does this mean that the GPPI analyses will remove the effects of motion regressors?

for i = 1:42

if i<10
precedingzero = '0';
else
precedingzero = '';
end

cd(horzcat('/Volumes/Seagate Backup Plus Drive/Bipolar EEG-fMRI/ARTREPAIR/', precedingzero, num2str(i), '/GLM/no_moody_EV'))

P.subject=horzcat(precedingzero, num2str(i));
P.directory= horzcat('/Volumes/Seagate Backup Plus Drive/Bipolar EEG-fMRI/ARTREPAIR/', precedingzero, num2str(i), '/GLM/no_moody_EV');
P.VOI= horzcat('/Volumes/Seagate Backup Plus Drive/Bipolar EEG-fMRI/ARTREPAIR/', precedingzero, num2str(i), '/GLM/no_moody_EV/LRstriatum.nii');
P.Region='VOI_LRstriatum';
P.Estimate=1;
P.contrast=9;
P.extract='eig';
P.Tasks={'1' 'Select' 'Antic' 'Outc'};
P.Weights=[];
P.analysis='psy';
P.method='cond';
P.CompContrasts=1;
P.Weighted=0;
P.ConcatR=0;
P.outdir = horzcat('/Volumes/Seagate Backup Plus Drive/Bipolar EEG-fMRI/ARTREPAIR/' , precedingzero, num2str(i), '/GLM/no_moody_EV/GPPI');
P.preservevarcorr=0;
P.wb=0;
P.equalroi=0;
P.FLmask=1;
%P.zipfiles=0;
%P.rWLS=0;
%P.FSFAST=0;
P.Contrasts(1).left={'Outc'};
P.Contrasts(1).right={'none'};
P.Contrasts(1).STAT='T';
P.Contrasts(1).Weighted=0;
P.Contrasts(1).MinEvents=1;
P.Contrasts(1).name='Outc';
P.Contrasts(2).left={'Outc'};
P.Contrasts(2).right={'none'};
P.Contrasts(2).Contrail={'xRPE^1'};
P.Contrasts(2).STAT='T';
P.Contrasts(2).Weighted=0;
P.Contrasts(2).MinEvents=1;
P.Contrasts(2).name='RPE';
P.Contrasts(3).left={'Outc'};
P.Contrasts(3).right={'none'};
P.Contrasts(3).Contrail={'xRPEdiff^1'};
P.Contrasts(3).STAT='T';
P.Contrasts(3).Weighted=0;
P.Contrasts(3).name='RPEdiff';
save(horzcat(precedingzero, num2str(i), '_LRstr_parameter.mat'), 'P')
PPPI(P)
clear

end