Removing scans in SPM
is difficult. But if you want to give it a go... First have a look at the message by Tom Nichols on the SPM list. You will need his fix to the spm_bch_tsampl.m file, attached to the mail; it is already installed at the CBU.Find the bad scans; lets say they are scans 50, 51 52
Do the slice timing as if there are two sessions for this one session: session 1 is scans 1 to 49, session 2 is 53 to the end.
Delete images 50, 51, 52. Replace them with another scan from the same series, maybe scan 49, with something like:
rm scan50.* scan51.* scan52.* cp scan49.img scan50.img cp scan49.hdr scan50.hdr cp scan49.mat scan50.mat cp scan49.img scan51.img cp scan49.hdr scan51.hdr cp scan49.mat scan51.mat cp scan49.img scan52.img cp scan49.hdr scan52.hdr cp scan49.mat scan52.mat
Do the realignment, normalization etc.Specify the design as if the scans were not missing. Then, in your batch file, do something like the following. Let's say that the scans were missing from session 2. Fill up your model batch structure as normal, then add something to specify the scans that are actually going to used:
% session scans that are missing misscans = { [], [50:52], [], [] }; % specify scans remaining in analysis for i = 1:4 remain = 1:scannos(i); remain(misscans{i}) = []; model(i).remain{i} = remain; end
Maybe that will work, who knows?Matthew 6/2/02 NRNR