When running 3dmvm with the -resid option having a large dataset (somehwere over 1500 subjects when having about 125 GB of RAM) it crashed due to memory error. I think it’s R that complains, I saw a smilar post here that said it was memory issues.
The stats are generated but not the residuals, so I am fairly sure that it’s the -resid part that crashes.
Would it be reasonable to make a separate directory called residuals and run 3dmvm twice with the table split in half, to generate residuals form the top and bottom half of the table and then average the two residuals volumes and feed it into 3dFWHMx?
Othwerwise I don’t think I can get the residuals. Running 3dFWHMx on each input file and average across subjects does in the end not yeild the same afc-parameters as if you do 3dMVM -resid. I confirmed this with a smaller sample. Since we are looking at structural maps (not the errts from fMRI data) it is important that we get the reisudals from mvm.