10 Jul
2013
10 Jul
'13
11:57 p.m.
The --cache is completely independent of --qcache (it was an unfortunate name for me to give it). The --cache means that the simulations have already been run. Think about a statistics text book that has a table of t-values, DOF, and corresponding p-values. I created a similar table that we distribute with FreeSurfer so that you don't have to run the simulation yourself. doug
On 7/10/13 11:52 PM, jh kim wrote:
Dear Prof. Greve
Thanks for the quick responses. My questions are now resolved.
Here is my last question regarding permutation testing. I cannot figure out what the '--cache' flag stands for. I've already done my data using recon-all with '-qcache' flag. Should I use '--cache' flag option in mri_glmfit-sim command?
Thank you in advance for your help.
Kim
2013/7/11 Douglas N Greve <greve@nmr.mgh.harvard.edu mailto:greve@nmr.mgh.harvard.edu>
On 07/09/2013 10:07 PM, jh kim wrote: Dear Prof. Greve I have additional (and probably silly) questions. 1) Is there much difference between result of monte-carlo null-z simulation (QDEC GUI) analysis and that of clusterwise correction for multiple comparisons using 10000 permutations? I'm guessing similar results between both methods; however I'm wondering whether significant clusters from the former processing might disappeared when using the latter process, or vice versa? When you say permutation, do you mean the permutation method where you permute your own data? It is hard to say whether they will give the same result. I usually find permutation to give less significant results (which may just be more accurate). 2) Does QDEC analysis automatically produce *.glmdir files for each contrast? If not, how can I make *.glmdir from the QDEC (not using the mri_glmfit command)? It creates a glmdir for the whole analysis. Within the glmdir, it creates a folder for each contrast doug Thank you for your kindness Kim 2013/7/9 Douglas N Greve <greve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> <mailto:greve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>>> On 07/07/2013 01:36 PM, jh kim wrote: > Dear FreeSurfer experts... > > Please be patient with a newbie's questions. > For the first time, I'm trying to analyze between-group differences in > cortical thickness (controls vs. patients) by using FS 5.1.0 and QDEC. > I've read through the online manual and mailing lists; however, I > cannot figure out the following queries. > > 1) I selected 'group' (con, pat) as a fixed factor (gender is not > applicable because all are females), and demeaned age as a continuous > covariate. After analyzing using DODS option, several questions are > shown in the Display GUI panel. Does the question "Does the average > thickness differ between con and pat?" mean between-group difference > of cortical thickness after *controlling for the effect of age*? Yes, it controls for all your other factors. > > 2) Should I include average cortical thickness of the left or right > cortex (calculated from the command, mris_anatomical_stats -l > lh.cortex.label -f subid.txt subid lh)as a nuisance variable? I know > this issue has been addressed many times before, and it is generally > recommended not to include as a nuisance factor. However, in case of > statistically significant difference in mean cortical thickness of one > hemisphere (unpaired two sample t-test by SPSS), is it still not > necessary to include it as a nuisance factor? From the similar points > of view, how about volume and surface area analyses (include or not > total cortical volume and total surface area as a nuisance factor, > respectively) ? I don't have a strong opinion about removing the global mean thickness, but I think others (Mike Harms?) recommend to do so. For volume and surface area, you probably want to normalize by the ICV. > > 3) I found several significant clusters that are corrected for > multiple comparisons using Monte Carlo simulation (corrected P < 0.01) > in the bottom of Display panel. I assume that this procedures take > 10,000 permutations by default. Am I right? It is 10k iterations of generating white noise (not a permutation test unless you specified that explicitly). > > 4) Last question, how can I extract individual value (cortical > thickness per subject) from the significant clusters? I want to do > this 'all-at-once'. Please let me know the command. There should already be a file in the output folder called something.y.ocn.dat that has a row for each input subject and a column for each cluster. doug > > Apology for bothering you. > Any comments would be greatly appreciated. > > Thanks. > Kim > > > _______________________________________________ > Freesurfer mailing list > Freesurfer@nmr.mgh.harvard.edu <mailto:Freesurfer@nmr.mgh.harvard.edu> <mailto:Freesurfer@nmr.mgh.harvard.edu <mailto:Freesurfer@nmr.mgh.harvard.edu>> > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer -- Douglas N. Greve, Ph.D. MGH-NMR Center greve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> <mailto:greve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>> Phone Number: 617-724-2358 <tel:617-724-2358> <tel:617-724-2358 <tel:617-724-2358>> Fax: 617-726-7422 <tel:617-726-7422> <tel:617-726-7422 <tel:617-726-7422>> Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2 www.nmr.mgh.harvard.edu/facility/filedrop/index.html <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/ _______________________________________________ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu <mailto:Freesurfer@nmr.mgh.harvard.edu> <mailto:Freesurfer@nmr.mgh.harvard.edu <mailto: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. -- Douglas N. Greve, Ph.D. MGH-NMR Center greve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> Phone Number: 617-724-2358 <tel:617-724-2358> Fax: 617-726-7422 <tel:617-726-7422> Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2 www.nmr.mgh.harvard.edu/facility/filedrop/index.html <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer