Significant in group analysis of what? The commands look ok
On 4/27/2021 9:18 AM, LIN, Yezhe wrote:
External Email - Use Caution
Dear Freesurfer experts,
My name is Yezhe Lin. I am trying to extract the thickness in an area we found significant in group analysis. However, when we did a t-test on the extracted values, the significance was not obtained. I wonder something could go wrong somewhere and what possibly caused the difference. I attached my scripts for the adjusted pipeline FYI.
steps #make a 3x3x3 ROI using the peak coordinate obtained in group analysis in standard space (MNI152) by fslview
#check ROI location tkmedit -f /usr/local/fsl/data/standard/MNI152_T1_2mm_brain.nii.gz -overlay MNI152_T1_2mm_rh.thickness.CBPHCagesexeTIVmask.nii -fthresh 0.5 # register to fsaverage of the dataset cd ../fsaverage/surf/ (Freesurfer_Home fsaverage) fslregister --s fsaverage --mov /mnt/backup2/restore-2021/CBPFreesurfer/FreesurferResults/TP1/FS/mask/MNI152_T1_2mm_brain.nii --reg MNI2mm_to_fsaverage.dat #To check results, run: tkregister2 --mov /mnt/backup2/restore-2021/CBPFreesurfer/FreesurferResults/TP1/FS/mask/MNI152_T1_2mm_brain.nii --reg MNI2mm_to_fsaverage.dat --surf orig #view your ROI in the fsaverage subject volume to ensure that it is located where you expect it to be: (two pials clear) cd $SUBJECTS_DIR/fsaverage/surf tkmedit fsaverage T1.mgz \ -overlay /mnt/backup2/restore-2021/CBPFreesurfer/FreesurferResults/TP1/FS/mask/MNI152_T1_2mm_rh.thickness.CBPHCagesexeTIVmask.nii \ -overlay-reg MNI2mm_to_fsaverage.dat \ -fthresh 0.5 \ -surface lh.white -aux-surface rh.white
#register ROI to fsaverage space cd $SUBJECTS_DIR/fsaverage/surf mri_vol2surf \ --mov /mnt/backup2/restore-2021/CBPFreesurfer/FreesurferResults/TP1/FS/mask/MNI152_T1_2mm_rh.thickness.CBPHCagesexeTIVmask.nii \ --reg MNI2mm_to_fsaverage.dat \ --projdist-max 0 1 0.1 \ --interp nearest \ --hemi rh \ --out rh.fsaverage.thickness.CBPHCagesexeTIVmask.mgh #view the ROI registration (on the peak center of significant cluster)
tksurfer fsaverage rh inflated -overlay rh.fsaverage.thickness.CBPHCagesexeTIVmask.mgh -fthresh 0.5
# extract values from rh.XXX.CBPvsHC_age_sex_eTIV.10.mgh mri_segstats --i ../rh.thickness.CBPvsHC_age_sex_eTIV.10.mgh --seg rh.fsaverage.thickness.CBPHCagesexeTIVmask.mgh --excludeid 0 --avgwf CBPHC_th_results.dat
Hope my questions is clear! Appreciate your kind help and time!
Best regards,
Yezhe
Yezhe Lin, MD MMed (she/her/hers)
Postdoctoral Associate, Department of Psychiatry
University of Rochester Medical Center
yezhe.lin.bee@gmail.com
Twitter https://secure-web.cisco.com/17T5lubn6nlIZ8y4sRQ9atwBlwoYN79xOySRhIHB7tYJ6OfM1OaOZAb7h0m5Gt1wqHzwRZiVxaLaolTIa0H34_obqegyt6jLRmWDX3RzMpeg1Bzj8zxIKcfvxodjYYc5j3LWMchkyWZx7w_nNy9yaA7_cy7LKrOadv-r1F3Cmc40fl7ZZAJgtP6cyMueZ5tw7iA4PAHIG-oy0tuxipyGIVSJABontdzFkTcRiLP42AAtF72Aeohb-tCKi99DzzI91X0-FnrgV6RKlVvh8JBWUYg/https%3A%2F%2Ftwitter.com%2FLinYezhe: @LinYezhe | LinkedIn https://secure-web.cisco.com/1NXc-SzQPbOIPzAqdQ4nxLbbQPvlXgxeaK7h2P-GMT2mOxEpslSv2Aq_oZZQNXHk_7WJCyWV8XVlZNBdvrOowUMj11c2f-ruIKsSFxw9bTS_R3GMqko5UnckotE8uZ8j2qeGvWBRnZBM8dYGqrvbX_FPxNCnBXzJcuSPXfzcPHQl1I-xgQiyryL_lxyP15CBmIkTb__ZQ-xnAXehOz81NKr8uP2nyeGaGWe2poxbnQfHaE-4Z7R5y5Togx_6i-RonP4tCQKFk1mRmj5r6ONjiRg/https%3A%2F%2Fwww.linkedin.com%2Fin%2Fyezhe-lin-md-mm-7063b2a1%2F | Researchgate https://secure-web.cisco.com/1h-lNaTYyp9GHwWyGnmRYHdr1jD9GARdf4R_WiOS0EG2qhl76derArpHYKcYGvgPZmloIjFzM7pazBATXCQ0_J2VQB_MytBAaLvK1LqPfyF-D8timX5eRK4UALiSq4PDKMvTzttgyxLhyTuAelo9fctNySNcEOsc5J7v1Vznc6Kzm2zm26bvJUchfZR3AULA-ZKDo5KlfYpg5C2A5V5Cjm7NvlEbibyGyrCRkI6AuFj5mCQoJ1mJFyO3eHUHwS5DgSw1CcC2nO-DZk87qzXQvVQ/https%3A%2F%2Fwww.researchgate.net%2Fprofile%2FYezhe_Lin
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer