Hi, I am working on R1 maps. I ran a very simple glm on smoothed R1 maps, FWHM = 8 mm, population 1, population 2 , and age as nuisance factor.
After POS MC simulation CW correction, I had two significant clusters resulting from this in the planum temporale and Heschl’s sulcus. The contrast being population 1 - population 2, my interpretation was higher R1 values in these two clusters in Pop 1 as compared to Pop 2. After that I wanted to run some post-hoc tests let say, in a statistical software similar to SPSS, JMP, on the mean R1 values from each of these clusters. The strategy I picked was to stay in the common space, fsaverage, and apply mri_segstats: #Extract mean R1 in montecarlo's significantcluster for i in `cat subjects` do mri_segstats --mul 1000 --annot fsaverage lh PT_mc-z.pos.sig.ocn --i $SUBJECTS_DIR/${i}/surf/lh.R1.fsaverage.fwhm8.mgh --sum lh.PT.R1.${i}.stats done for i in `cat subjects` do mri_segstats --mul 1000 --annot fsaverage lh HS_mc-z.pos.sig.ocn --i $SUBJECTS_DIR/${i}/surf/lh.R1.fsaverage.fwhm8.mgh --sum lh.HS.R1.${i}.stats done Everything when fine but as compared to my glm analysis when looking at mean group differences in my stat software I have the opposite : larger R1 values in Pop 2 as compared to Pop 1 (after removing the effect of age or not, and removing one outlier or not). I went back to my glm and checked everything, fsgd file, matrix and so on, everything seems ok I am not mistaken on the direction of my glm contrast (pop 1 - pop 2) so I do not get it at all. Any advices? For now I am thinking about going back in the individual space, registering my clusters in the individual surfaces (label2surf something like that) and then extract my mean R1 values, or extract values in non-smoothed maps, but I do not think it will change anything. Damien
_______________________________________________ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Partners Compliance HelpLine at http://www.partners.org/complianceline . If the e-mail was sent to you in error but does not contain patient information, please contact the sender and properly dispose of the e-mail.