Dear Freesurfer experts,
We are using mri_glmfit and wanted to compare our results to qdec to
decide which method we prefer for group analysis. The results are
consistent between the two methods, except for the monte carlo results.
When we run mri_glmfit-sim to do the monte carlo simulations there are
no significant clusters. We are not surprised by this as when we apply a
FDR corrections all of our clusters go away. However, when we are using
the monte carlo method in qdec we get a large significant cluster (in
the csd summary file there are two large significant clusters, but using
the image in qdec we only see one of them). The model is the same for
both methods. We are (or think) we are using the exact same parameters
as qdec. Is there a reason why the simulation results are so different?
Which result should we use?
The following are two of the simulation commands we have tried using v5.1:
mri_glmfit-sim --glmdir glm/lh.asrs_sex_age.thickness.glmdir/ --sim mc-z
100 4 mc-z.unsigned --sim-sign abs --cwpvalthresh 0.01 --overwrite
mri_glmfit-sim --glmdir glm/lh.asrs_sex_age.thickness.glmdir/ --sim mc-z
100 2 mc-z.unsigned --sim-sign abs --cwpvalthresh 0.01 --overwrite
For the following command, we saw three clusters (they are inconsistent
with the qdec ones), but none of them significant (p-values=.2):
mri_glmfit-sim --glmdir glm/lh.asrs_sex_age.thickness.glmdir/ --sim mc-z
100 4 mc-z.unsigned --sim-sign abs --cwpvalthresh 0.995 --overwrite
We haven't run 10 000 simulations yet, because we wanted to see if there
were any significant clusters there at all before running that many
simulations.
Any advice would be very helpful.
Thank you for your time,
Marie
_______________________________________________
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.
