...
| Code Block |
|---|
#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
|
After carrying out statistical analysis with QDEC or the glmfit command, you can correct the results for multiple comparisons using the below code. Please refer to the link for descriptions of each option, how to change the vertex-wise and cluster-wise thresholds, and how to do a one-sided test.
| Code Block |
|---|
|
#PBS -P
#PBS -N FScorrection
#PBS -l select=1:ncpus=4:mem=32GB
#PBS -l walltime=2:00:00
#PBS -q defaultQ
module load freesurfer/6.0.0
source $FREESURFER_HOME/SetUpFreeSurfer.sh
##Add the subject directory location and the location of all the glm directories below##
export SUBJECTS_DIR=
glmdir=/scratch/long_fixel_FTD/CorticalThickness/Data/qdec/
cd $SUBJECTS_DIR
list=$(ls ${glmdir})
for i in $list_ ; do
mri_glmfit-sim \
--glmdir ${glmdir}${i} \
--perm 1000 3.0 abs \
--cwp 0.05 \
--2spaces \
--bg 4 ;
done
|