I am trying to do the clusterwise correction for multiple comparisons and came upon a problem that maybe someone can help me with.
ri_volcluster --in brianweisinger/g1g2.intercept/sig.mgh --csd brianweisinger/csd/mc-z.negative.j001-g1g2.intercept.csd --mask brianweisinger/mask.mgh --cwsig brianweisinger/g1g2.intercept/mc-z.negative.sig.cluster.mgh --vwsig brianweisinger/g1g2.intercept/mc-z.negative.sig.voxel.mgh --sum brianweisinger/g1g2.intercept/mc-z.negative.sig.cluster.summary --ocn brianweisinger/g1g2.intercept/mc-z.negative.sig.ocn.mgh --csdpdf brianweisinger/g1g2.intercept/mc-z.negative.pdf.dat Adjusting threshold for 1-tailed test. If the input is not a -log10(p) volume, re-run with --no-adjust. INFO: loading mask volume: brianweisinger/mask.mgh INFO: MRIbinarize01(): nhits = 118078 Found 118078 voxels in mask Search Space = 118078 mm3 Computing voxel-wise significance CSDpvalMaxSigMap(): found 0/118078 above 0 ERROR: initializing hit map
It made some of the CSD files but not all of them. Any ideas?
Thank you for your help.
Looks like it did not find any vertices above the threshold that you used for the simulation (though it should give a better error msg). Before you run the simulation, you can look at the map at your given threshold. If there are no vertices above threshold, then nothing is going to show up in the simulation.
doug
Weisinger, Brian (NIH/OD) [E] wrote:
I am trying to do the clusterwise correction for multiple comparisons and came upon a problem that maybe someone can help me with.
ri_volcluster --in brianweisinger/g1g2.intercept/sig.mgh --csd brianweisinger/csd/mc-z.negative.j001-g1g2.intercept.csd --mask brianweisinger/mask.mgh --cwsig brianweisinger/g1g2.intercept/mc-z.negative.sig.cluster.mgh --vwsig brianweisinger/g1g2.intercept/mc-z.negative.sig.voxel.mgh --sum brianweisinger/g1g2.intercept/mc-z.negative.sig.cluster.summary --ocn brianweisinger/g1g2.intercept/mc-z.negative.sig.ocn.mgh --csdpdf brianweisinger/g1g2.intercept/mc-z.negative.pdf.dat Adjusting threshold for 1-tailed test. If the input is not a -log10(p) volume, re-run with --no-adjust. INFO: loading mask volume: brianweisinger/mask.mgh INFO: MRIbinarize01(): nhits = 118078 Found 118078 voxels in mask Search Space = 118078 mm3 Computing voxel-wise significance CSDpvalMaxSigMap(): found 0/118078 above 0 ERROR: initializing hit map
It made some of the CSD files but not all of them. Any ideas?
Thank you for your help.
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
freesurfer@nmr.mgh.harvard.edu