Hello,
I'm running a cortical thickness analysis using Qdec. When I do Find Clusters and Go to Max, I obtain 9 clusters (min threshold 3, max threshold 5): ============================================================ Generating cluster stats using min threshold of 3... Found 9 clusters Contrast: 'lh-Diff-grupo1-grupo2-Intercept-thickness', 10fwhm, DOF: 30 ClusterNo Max VtxMax Size(mm2) TalX TalY TalZ NVtxs Annotation --------- --- ------ --------- ---- ---- ---- ----- ---------- 1 8.8371 89238 5.07 -28.2 12.0 -14.5 16 insula 2 -5.3729 64233 428.45 -35.4 50.5 -8.5 583 parsorbitalis 3 -5.1786 73273 192.98 -8.0 46.4 -13.0 292 medialorbitofrontal 4 -3.8763 121706 38.22 -19.6 61.0 1.5 53 rostralmiddlefrontal 5 -3.8612 93042 93.07 -23.8 40.2 -11.0 153 lateralorbitofrontal 6 -3.5732 14941 41.60 -7.0 41.1 -22.2 79 medialorbitofrontal 7 3.3511 160348 0.34 -30.3 10.7 -14.1 1 insula 8 3.3502 111085 31.96 -33.8 -13.2 -4.6 68 insula 9 -3.2412 99299 19.05 -34.6 47.2 18.4 26 rostralmiddlefrontal ============================================================
Why when using the command line mri_surfcluster with the same thresholds I do not have the same results (same clusters)?
# Cluster Growing Summary (mri_surfcluster) # $Id: mri_surfcluster.c,v 1.48 2010/04/12 17:36:58 greve Exp $ # $Id: mrisurf.c,v 1.678 2010/05/28 21:22:21 rpwang Exp $ # CreationTime 2011/05/02-19:07:52-GMT # cmdline mri_surfcluster --in qdec/lh_g1_g2/lh-Diff-grupo1-grupo2-Intercept-thickness/sig.mgh --subject fsaverage --hemi lh --thmin 3 --thmax 5 --sign abs --annot aparc --sum summary_3_5 --olab lh_g1_g2 # cwd /Applications/freesurfer/subjects/dados # sysname Darwin # hostname Fernanda-Palhanos-Mac-Pro.local # machine i386 # FixVertexAreaFlag 1 # FixSurfClusterArea 1 # # Input qdec/lh_g1_g2/lh-Diff-grupo1-grupo2-Intercept-thickness/sig.mgh # Frame Number 0 # srcsubj fsaverage # hemi lh # surface white # annot aparc # SUBJECTS_DIR /Applications/freesurfer/subjects/dados # Minimum Threshold 3 # Maximum Threshold 5 # Threshold Sign abs # AdjustThreshWhenOneTail 1 # Area Threshold 0 mm^2 # Overall max 8.83713 at vertex 89238 # Overall min -5.37286 at vertex 64233 # NClusters 12 # Total Cortical Surface Area 65416.6 (mm^2) # FixMNI = 1 # # ClusterNo Max VtxMax Size(mm^2) TalX TalY TalZ NVtxs Annot 1 -4.995 29209 398.78 -33.5 48.6 -10.8 548 rostralmiddlefrontal 2 -4.995 16744 185.39 -8.3 45.3 -12.5 283 medialorbitofrontal 3 4.914 28205 0.34 -28.4 11.8 -11.8 1 insula 4 4.850 22061 2.37 -26.8 9.7 -13.8 8 insula 5 4.512 52766 0.38 -29.8 11.0 -11.4 1 insula 6 -3.876 121706 38.22 -19.4 59.1 -1.6 53 rostralmiddlefrontal 7 -3.861 93042 93.07 -23.5 38.5 -11.2 153 lateralorbitofrontal 8 3.691 89239 0.68 -27.5 11.7 -12.5 2 insula 9 -3.573 14941 41.60 -7.0 38.8 -20.6 79 medialorbitofrontal 10 3.351 160348 0.34 -29.9 9.7 -12.3 1 insula 11 3.350 111085 31.96 -33.5 -12.9 -3.2 68 insula 12 -3.241 99299 19.05 -34.2 46.6 14.7 26 rostralmiddlefrontal
In the right hemisphere I have the same problem:
============================================================ Generating cluster stats using min threshold of 3... Found 6 clusters Contrast: 'rh-Diff-grupo1-grupo2-Intercept-thickness', 10fwhm, DOF: 30 ClusterNo Max VtxMax Size(mm2) TalX TalY TalZ NVtxs Annotation --------- --- ------ --------- ---- ---- ---- ----- ---------- 1 16.5175 22102 2.89 33.5 10.7 -13.0 7 insula 2 6.0843 160464 1.08 30.8 13.1 -13.1 3 insula 3 -5.9997 62430 1.08 36.1 3.3 -19.1 3 insula 4 -3.3483 60593 24.75 35.9 55.5 -7.6 36 rostralmiddlefrontal 5 3.2156 23631 10.78 19.8 -65.1 1.5 15 lingual 6 3.1968 120630 0.28 29.1 11.5 -15.8 1 insula ============================================================
cmdline mri_surfcluster --in qdec/rh_g1_g2/rh-Diff-grupo1-grupo2-Intercept-thickness/sig.mgh --subject fsaverage --hemi rh --thmin 3 --thmax 5 --sign abs --annot aparc --sum rh_summary_3_5 --olab rh_g1_g2 # cwd /Applications/freesurfer/subjects/dados # sysname Darwin # hostname Fernanda-Palhanos-Mac-Pro.local # machine i386 # FixVertexAreaFlag 1 # FixSurfClusterArea 1 # # Input qdec/rh_g1_g2/rh-Diff-grupo1-grupo2-Intercept-thickness/sig.mgh # Frame Number 0 # srcsubj fsaverage # hemi rh # surface white # annot aparc # SUBJECTS_DIR /Applications/freesurfer/subjects/dados # Minimum Threshold 3 # Maximum Threshold 5 # Threshold Sign abs # AdjustThreshWhenOneTail 1 # Area Threshold 0 mm^2 # Overall max 16.5175 at vertex 22102 # Overall min -5.99973 at vertex 62430 # NClusters 6 # Total Cortical Surface Area 65020.8 (mm^2) # FixMNI = 1 # # ClusterNo Max VtxMax Size(mm^2) TalX TalY TalZ NVtxs Annot 1 3.386 130519 0.70 29.7 12.4 -12.0 2 insula 2 -3.348 60593 24.75 35.5 53.5 -9.1 36 rostralmiddlefrontal 3 -3.245 28242 0.36 35.3 2.9 -15.6 1 insula 4 3.216 23631 10.78 19.6 -63.0 4.6 15 lingual 5 3.197 120630 0.28 28.8 10.4 -13.8 1 insula 6 3.067 89363 0.38 34.3 7.1 -12.2 1 insula
Am I doing something wrong? I want to use the command line to generate the labels for each cluster found.
Thanks in advance, Fernanda