FreeSurfer group analysis can either be carried out with the graphical interface, Qdec, or by creating your own GLM and contrasts and running mri_glmfit on the command line. Instructions for both these methods can be found on the FreeSurfer wiki pages linked above. In both cases you will need to smooth the data first. This can be done with the following script on Artemis.
Before running the script you will need to add in the project name, SUBJECTS_DIR, and a subject list.
#PBS -P
#PBS -N QDECproc
#PBS -l select=1:ncpus=2:mem=32GB
#PBS -q=main
#PBS -l walltime=10:00:00
#PBS -q defaultQ
module load mrtrix3/3.0.1
module load freesurfer/6.0.0
source $FREESURFER_HOME/SetUpFreeSurfer.sh
export SUBJECTS_DIR=/scratch/Project/Folder
cd $SUBJECTS_DIR
sub_list='FTDXXXX
FTDXXXX
FTDXXXX
'
cp -r /usr/local/freesurfer/6.0.0/subjects/fsaverage ${SUBJECTS_DIR}
for i in $sub_list ; do
mv ${SUBJECTS_DIR}/${i}/${i}_FreeSurferAnat/* ${SUBJECTS_DIR}/${i}/;
mris_preproc --s ${i} --hemi lh --meas thickness --fwhm 15 --target fsaverage --out ${i}/surf/lh.thickness.fwhm15.fsaverage.mgh;
mris_preproc --s ${i} --hemi rh --meas thickness --fwhm 15 --target fsaverage --out ${i}/surf/rh.thickness.fwhm15.fsaverage.mgh;
mris_preproc --s ${i} --hemi lh --meas thickness --fwhm 20 --target fsaverage --out ${i}/surf/lh.thickness.fwhm20.fsaverage.mgh;
mris_preproc --s ${i} --hemi rh --meas thickness --fwhm 20 --target fsaverage --out ${i}/surf/rh.thickness.fwhm20.fsaverage.mgh;
done
0 Comments