Re: [Freesurfer] glmfit with aseg table: view surface overlay
Not easily. You could do something in matlab, eg
aseg = MRIread('aseg.mgz');
sig = aseg;
sig.vol = zeros(seg.volsize);
ind = find(aseg.vol == 17); % left hippocampus
sig.vol(ind) = SigValueForLeftHippo; % from your table
% Do the same with the rest of your segmentations
MRIwrite(sig,'mysig.mgz')
On 4/16/2021 4:09 PM, Griem, Julia wrote:
External Email - Use Caution
External Email - Use Caution
Hi Doug,
Ok, that’s fine. Is it possible to visualize subcortical volume
differences in any way?
Thank you!
Julia
/NIHR Maudsley Biomedical Research Centre PhD Student/
/Forensic & Neurodevelopmental Sciences/
/Institute of Psychiatry, Psychology & Neuroscience/
/King’s College London/
/16 De Crespigny Park, Camberwell/
/London, SE5 8AB/
*From:*freesurfer-boun...@nmr.mgh.harvard.edu
*On Behalf Of *Douglas N. Greve
*Sent:* 16 April 2021 16:11
*To:* freesurfer@nmr.mgh.harvard.edu
*Subject:* Re: [Freesurfer] glmfit with aseg table: view surface overlay
Are you trying to correct the table? You can't do that with
mri_gtlmfit-sim. That looks for clusters, which does not make sense
when you are using table data. There are a couple of ways you can do
it. One is with bonerfonni correction (ie, multiply the p-value by the
number of ROIs) or with FDR.
On 4/15/2021 11:13 AM, Griem, Julia wrote:
*External Email - Use Caution *
*External Email - Use Caution *
Hi Doug,
Thank you. Makes sense. Is there any way to visualize significance
overlays for the subcortical volumes?
I’ve also moved on to the correction for multiple comparisons in
my subcortical volumes, using:
mri_glmfit-sim -- glmdir roi.binarygroup.glmdir --perm 1000 1.3
abs --cwp 0.05 --2spaces --bg 1
This did not create any new files in the Contrast subdirectories,
but the log file does say that it was run and finished. Can I find
a table with corrected p-values anywhere? Or would it make more
sense to move the aseg_table into a different stats package and
calculate my GLM there, correcting for multiple comparisons
post-hoc (manually)?
Thank you so much
Julia
/NIHR Maudsley Biomedical Research Centre PhD Student/
/Forensic & Neurodevelopmental Sciences/
/Institute of Psychiatry, Psychology & Neuroscience/
/King’s College London/
/16 De Crespigny Park, Camberwell/
/London, SE5 8AB/
*From:*freesurfer-boun...@nmr.mgh.harvard.edu
<mailto:freesurfer-boun...@nmr.mgh.harvard.edu>
<mailto:freesurfer-boun...@nmr.mgh.harvard.edu> *On Behalf Of
*Douglas N. Greve
*Sent:* 15 April 2021 15:14
*To:* freesurfer@nmr.mgh.harvard.edu
<mailto:freesurfer@nmr.mgh.harvard.edu>
*Subject:* Re: [Freesurfer] glmfit with aseg table: view surface
overlay
The input is a table, so the output is a table (the sig.mgh file
will just have the number of values that are in the table; this is
not a "volume" you can look at in FV). Look in the sig.table.dat
file for the sigs.
On 4/13/2021 8:41 AM, Griem, Julia wrote:
*External Email - Use Caution *
*External Email - Use Caution *
Dear Freesurfer team,
I am comparing subcortical volumes (aseg table) in 2 groups
with no other variables. Based on the GLM lecture slides from
the course I attended, I ran the following glmfit command:
mri_glmfit
--table aseg_table.txt
--fsgd FSGD/DemographicsSubcortBinaryGroup.fsgd
--C Contrasts/HC-ASPD.mtx
--C Contrasts/ASPD-HC.mtx
--glmdir roi.binarygroup.glmdir
It created all the correct output and subdirectories.
I then used
freeview -f
$SUBJECTS_DIR/fsaverage/surf/lh.inflated:overlay=ASPD-HC/sig.mgh
but an error saying "failed to load surface overlay" comes up
and the terminal error message is "number of vertices in
overlay data does not match with surface". I assume this is
because I am trying to load a 'whole-brain' sig.mgh onto a lh
surface? Do I have to change the glmfit command to specify
hemispheres, and if so, how do I do that since I do not want
to look at cortical but subcortical data? Or alternatively,
can I load both hemispheres into freeview so the overlay
sig.mgh fits?
Thank you so much,
Julia
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu <mailto:Freesurfer@nmr.mgh.harvard.edu>
https://secure-web.cisco.com/1lCp4yPi_vnD1_jDHjqMFFlcPKYlBFqg-Dg5nn2IbJrmH-ru9ATLtvP78StdglWKXz9Aoe2wkuZl3qWXtUlkzaPRyFPRyN-OfCI2-hTiv-3a80hsVv3ZU0pu_dlM-mpTTSgILcHfolVQX-BFcoHymMYXlXOQy8cuIBiKE6WSwV49wP1b003OX5yqYIduz3OPXVBpXTw-rkyt1naV1kGogJEbRs3QNS8dOuep2scIGZ
Re: [Freesurfer] glmfit with aseg table: view surface overlay
External Email - Use Caution Hi Doug, Ok, that's fine. Is it possible to visualize subcortical volume differences in any way? Thank you! Julia NIHR Maudsley Biomedical Research Centre PhD Student Forensic & Neurodevelopmental Sciences Institute of Psychiatry, Psychology & Neuroscience King's College London 16 De Crespigny Park, Camberwell London, SE5 8AB From: freesurfer-boun...@nmr.mgh.harvard.edu On Behalf Of Douglas N. Greve Sent: 16 April 2021 16:11 To: freesurfer@nmr.mgh.harvard.edu Subject: Re: [Freesurfer] glmfit with aseg table: view surface overlay Are you trying to correct the table? You can't do that with mri_gtlmfit-sim. That looks for clusters, which does not make sense when you are using table data. There are a couple of ways you can do it. One is with bonerfonni correction (ie, multiply the p-value by the number of ROIs) or with FDR. On 4/15/2021 11:13 AM, Griem, Julia wrote: External Email - Use Caution External Email - Use Caution Hi Doug, Thank you. Makes sense. Is there any way to visualize significance overlays for the subcortical volumes? I've also moved on to the correction for multiple comparisons in my subcortical volumes, using: mri_glmfit-sim -- glmdir roi.binarygroup.glmdir --perm 1000 1.3 abs --cwp 0.05 --2spaces --bg 1 This did not create any new files in the Contrast subdirectories, but the log file does say that it was run and finished. Can I find a table with corrected p-values anywhere? Or would it make more sense to move the aseg_table into a different stats package and calculate my GLM there, correcting for multiple comparisons post-hoc (manually)? Thank you so much Julia NIHR Maudsley Biomedical Research Centre PhD Student Forensic & Neurodevelopmental Sciences Institute of Psychiatry, Psychology & Neuroscience King's College London 16 De Crespigny Park, Camberwell London, SE5 8AB From: freesurfer-boun...@nmr.mgh.harvard.edu<mailto:freesurfer-boun...@nmr.mgh.harvard.edu> <mailto:freesurfer-boun...@nmr.mgh.harvard.edu> On Behalf Of Douglas N. Greve Sent: 15 April 2021 15:14 To: freesurfer@nmr.mgh.harvard.edu<mailto:freesurfer@nmr.mgh.harvard.edu> Subject: Re: [Freesurfer] glmfit with aseg table: view surface overlay The input is a table, so the output is a table (the sig.mgh file will just have the number of values that are in the table; this is not a "volume" you can look at in FV). Look in the sig.table.dat file for the sigs. On 4/13/2021 8:41 AM, Griem, Julia wrote: External Email - Use Caution External Email - Use Caution Dear Freesurfer team, I am comparing subcortical volumes (aseg table) in 2 groups with no other variables. Based on the GLM lecture slides from the course I attended, I ran the following glmfit command: mri_glmfit --table aseg_table.txt --fsgd FSGD/DemographicsSubcortBinaryGroup.fsgd --C Contrasts/HC-ASPD.mtx --C Contrasts/ASPD-HC.mtx --glmdir roi.binarygroup.glmdir It created all the correct output and subdirectories. I then used freeview -f $SUBJECTS_DIR/fsaverage/surf/lh.inflated:overlay=ASPD-HC/sig.mgh but an error saying "failed to load surface overlay" comes up and the terminal error message is "number of vertices in overlay data does not match with surface". I assume this is because I am trying to load a 'whole-brain' sig.mgh onto a lh surface? Do I have to change the glmfit command to specify hemispheres, and if so, how do I do that since I do not want to look at cortical but subcortical data? Or alternatively, can I load both hemispheres into freeview so the overlay sig.mgh fits? Thank you so much, Julia ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu<mailto:Freesurfer@nmr.mgh.harvard.edu> https://secure-web.cisco.com/1lCp4yPi_vnD1_jDHjqMFFlcPKYlBFqg-Dg5nn2IbJrmH-ru9ATLtvP78StdglWKXz9Aoe2wkuZl3qWXtUlkzaPRyFPRyN-OfCI2-hTiv-3a80hsVv3ZU0pu_dlM-mpTTSgILcHfolVQX-BFcoHymMYXlXOQy8cuIBiKE6WSwV49wP1b003OX5yqYIduz3OPXVBpXTw-rkyt1naV1kGogJEbRs3QNS8dOuep2scIGZBsiejUTWQ70ctg1Q1aHGxp0cDDmqUiw5JiB4-w-xtjucQ/https%3A%2F%2Fmail.nmr.mgh.harvard.edu%2Fmailman%2Flistinfo%2Ffreesurfer<https://secure-web.cisco.com/15TPPgYUIAbq_mfehNPq_gFJ_A4L6U-A3i04XgnxgmRj4TXZHT7RWtUEzj2_mNIXg7bUpa2Wu_ZHYaZtDt9OlhcBMTKEt-f5caWg2qrBP4NAnkSCCrSF6xkh3tdlNnoO4TWHCYFBatZRbnE5f1my2rp3jyKmNVeKI1_FtMalmWxsq01gC-UyJiqJ-cc7zEIOYupXJlSaS1OIQvD8uH5qQwiUH7mlb-M8e1yXuYb6PUEE/https%3A%2F%2Feur03.safelinks.protection.outlook.com%2F%3Furl%3Dhttps%253A%252F%252Fsecure-web.cisco.com%252F1p4hRmT0vowWhGIfUeLnb_VJVgx9U6ZITRXr63s2lm8_k7kAez1k6R-KysBh1p63Yb2XHYAg8vuquiA9J6e89ahoYjRji0z1mc9YRR_I_qzsg0h-M5P_pG46-KieVhwAW8f4CV0uxI0YGWj3l8ztd7Lf_noyMpRRe1cPMBnuATLarbyVP_qjP4CkSx6Q_0-NjnKLXM1EK1rlqJg-MFoBrxRM3c59ypxnMuMekTfqQpTPwXMtlGEVxSHbgM2hGwnr6WsSuCmNuOer7HDHskpvukA%252Fhttps%25253A%25252F%25252Feur03.
Re: [Freesurfer] glmfit with aseg table: view surface overlay
Are you trying to correct the table? You can't do that with mri_gtlmfit-sim. That looks for clusters, which does not make sense when you are using table data. There are a couple of ways you can do it. One is with bonerfonni correction (ie, multiply the p-value by the number of ROIs) or with FDR. On 4/15/2021 11:13 AM, Griem, Julia wrote: External Email - Use Caution External Email - Use Caution Hi Doug, Thank you. Makes sense. Is there any way to visualize significance overlays for the subcortical volumes? I’ve also moved on to the correction for multiple comparisons in my subcortical volumes, using: mri_glmfit-sim -- glmdir roi.binarygroup.glmdir --perm 1000 1.3 abs --cwp 0.05 --2spaces --bg 1 This did not create any new files in the Contrast subdirectories, but the log file does say that it was run and finished. Can I find a table with corrected p-values anywhere? Or would it make more sense to move the aseg_table into a different stats package and calculate my GLM there, correcting for multiple comparisons post-hoc (manually)? Thank you so much Julia /NIHR Maudsley Biomedical Research Centre PhD Student/ /Forensic & Neurodevelopmental Sciences/ /Institute of Psychiatry, Psychology & Neuroscience/ /King’s College London/ /16 De Crespigny Park, Camberwell/ /London, SE5 8AB/ *From:*freesurfer-boun...@nmr.mgh.harvard.edu *On Behalf Of *Douglas N. Greve *Sent:* 15 April 2021 15:14 *To:* freesurfer@nmr.mgh.harvard.edu *Subject:* Re: [Freesurfer] glmfit with aseg table: view surface overlay The input is a table, so the output is a table (the sig.mgh file will just have the number of values that are in the table; this is not a "volume" you can look at in FV). Look in the sig.table.dat file for the sigs. On 4/13/2021 8:41 AM, Griem, Julia wrote: *External Email - Use Caution * *External Email - Use Caution * Dear Freesurfer team, I am comparing subcortical volumes (aseg table) in 2 groups with no other variables. Based on the GLM lecture slides from the course I attended, I ran the following glmfit command: mri_glmfit --table aseg_table.txt --fsgd FSGD/DemographicsSubcortBinaryGroup.fsgd --C Contrasts/HC-ASPD.mtx --C Contrasts/ASPD-HC.mtx --glmdir roi.binarygroup.glmdir It created all the correct output and subdirectories. I then used freeview -f $SUBJECTS_DIR/fsaverage/surf/lh.inflated:overlay=ASPD-HC/sig.mgh but an error saying "failed to load surface overlay" comes up and the terminal error message is "number of vertices in overlay data does not match with surface". I assume this is because I am trying to load a 'whole-brain' sig.mgh onto a lh surface? Do I have to change the glmfit command to specify hemispheres, and if so, how do I do that since I do not want to look at cortical but subcortical data? Or alternatively, can I load both hemispheres into freeview so the overlay sig.mgh fits? Thank you so much, Julia ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu <mailto:Freesurfer@nmr.mgh.harvard.edu> https://secure-web.cisco.com/1lCp4yPi_vnD1_jDHjqMFFlcPKYlBFqg-Dg5nn2IbJrmH-ru9ATLtvP78StdglWKXz9Aoe2wkuZl3qWXtUlkzaPRyFPRyN-OfCI2-hTiv-3a80hsVv3ZU0pu_dlM-mpTTSgILcHfolVQX-BFcoHymMYXlXOQy8cuIBiKE6WSwV49wP1b003OX5yqYIduz3OPXVBpXTw-rkyt1naV1kGogJEbRs3QNS8dOuep2scIGZBsiejUTWQ70ctg1Q1aHGxp0cDDmqUiw5JiB4-w-xtjucQ/https%3A%2F%2Fmail.nmr.mgh.harvard.edu%2Fmailman%2Flistinfo%2Ffreesurfer <https://secure-web.cisco.com/1p4hRmT0vowWhGIfUeLnb_VJVgx9U6ZITRXr63s2lm8_k7kAez1k6R-KysBh1p63Yb2XHYAg8vuquiA9J6e89ahoYjRji0z1mc9YRR_I_qzsg0h-M5P_pG46-KieVhwAW8f4CV0uxI0YGWj3l8ztd7Lf_noyMpRRe1cPMBnuATLarbyVP_qjP4CkSx6Q_0-NjnKLXM1EK1rlqJg-MFoBrxRM3c59ypxnMuMekTfqQpTPwXMtlGEVxSHbgM2hGwnr6WsSuCmNuOer7HDHskpvukA/https%3A%2F%2Feur03.safelinks.protection.outlook.com%2F%3Furl%3Dhttps%253A%252F%252Fsecure-web.cisco.com%252F1lCp4yPi_vnD1_jDHjqMFFlcPKYlBFqg-Dg5nn2IbJrmH-ru9ATLtvP78StdglWKXz9Aoe2wkuZl3qWXtUlkzaPRyFPRyN-OfCI2-hTiv-3a80hsVv3ZU0pu_dlM-mpTTSgILcHfolVQX-BFcoHymMYXlXOQy8cuIBiKE6WSwV49wP1b003OX5yqYIduz3OPXVBpXTw-rkyt1naV1kGogJEbRs3QNS8dOuep2scIGZBsiejUTWQ70ctg1Q1aHGxp0cDDmqUiw5JiB4-w-xtjucQ%252Fhttps%25253A%25252F%25252Fmail.nmr.mgh.harvard.edu%25252Fmailman%25252Flistinfo%25252Ffreesurfer%26data%3D04%257C01%257Cjulia.griem%2540kcl.ac.uk%257C3f5e76a919d749738d0808d90018e9b9%257C8370cf1416f34c16b83c724071654356%257C0%257C0%257C637540929284868052%257CUnknown%257CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%253D%257C1000%26sdata%3DsjiRm%252BXt2HfAkzyBe4V8fixcdeM0VPaJDFm11VM2g%252BE%253D%26reserved%3D0> ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://secur
Re: [Freesurfer] glmfit with aseg table: view surface overlay
External Email - Use Caution Hi Doug, Thank you. Makes sense. Is there any way to visualize significance overlays for the subcortical volumes? I've also moved on to the correction for multiple comparisons in my subcortical volumes, using: mri_glmfit-sim -- glmdir roi.binarygroup.glmdir --perm 1000 1.3 abs --cwp 0.05 --2spaces --bg 1 This did not create any new files in the Contrast subdirectories, but the log file does say that it was run and finished. Can I find a table with corrected p-values anywhere? Or would it make more sense to move the aseg_table into a different stats package and calculate my GLM there, correcting for multiple comparisons post-hoc (manually)? Thank you so much Julia NIHR Maudsley Biomedical Research Centre PhD Student Forensic & Neurodevelopmental Sciences Institute of Psychiatry, Psychology & Neuroscience King's College London 16 De Crespigny Park, Camberwell London, SE5 8AB From: freesurfer-boun...@nmr.mgh.harvard.edu On Behalf Of Douglas N. Greve Sent: 15 April 2021 15:14 To: freesurfer@nmr.mgh.harvard.edu Subject: Re: [Freesurfer] glmfit with aseg table: view surface overlay The input is a table, so the output is a table (the sig.mgh file will just have the number of values that are in the table; this is not a "volume" you can look at in FV). Look in the sig.table.dat file for the sigs. On 4/13/2021 8:41 AM, Griem, Julia wrote: External Email - Use Caution External Email - Use Caution Dear Freesurfer team, I am comparing subcortical volumes (aseg table) in 2 groups with no other variables. Based on the GLM lecture slides from the course I attended, I ran the following glmfit command: mri_glmfit --table aseg_table.txt --fsgd FSGD/DemographicsSubcortBinaryGroup.fsgd --C Contrasts/HC-ASPD.mtx --C Contrasts/ASPD-HC.mtx --glmdir roi.binarygroup.glmdir It created all the correct output and subdirectories. I then used freeview -f $SUBJECTS_DIR/fsaverage/surf/lh.inflated:overlay=ASPD-HC/sig.mgh but an error saying "failed to load surface overlay" comes up and the terminal error message is "number of vertices in overlay data does not match with surface". I assume this is because I am trying to load a 'whole-brain' sig.mgh onto a lh surface? Do I have to change the glmfit command to specify hemispheres, and if so, how do I do that since I do not want to look at cortical but subcortical data? Or alternatively, can I load both hemispheres into freeview so the overlay sig.mgh fits? Thank you so much, Julia ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu<mailto:Freesurfer@nmr.mgh.harvard.edu> https://secure-web.cisco.com/1lCp4yPi_vnD1_jDHjqMFFlcPKYlBFqg-Dg5nn2IbJrmH-ru9ATLtvP78StdglWKXz9Aoe2wkuZl3qWXtUlkzaPRyFPRyN-OfCI2-hTiv-3a80hsVv3ZU0pu_dlM-mpTTSgILcHfolVQX-BFcoHymMYXlXOQy8cuIBiKE6WSwV49wP1b003OX5yqYIduz3OPXVBpXTw-rkyt1naV1kGogJEbRs3QNS8dOuep2scIGZBsiejUTWQ70ctg1Q1aHGxp0cDDmqUiw5JiB4-w-xtjucQ/https%3A%2F%2Fmail.nmr.mgh.harvard.edu%2Fmailman%2Flistinfo%2Ffreesurfer<https://secure-web.cisco.com/1p4hRmT0vowWhGIfUeLnb_VJVgx9U6ZITRXr63s2lm8_k7kAez1k6R-KysBh1p63Yb2XHYAg8vuquiA9J6e89ahoYjRji0z1mc9YRR_I_qzsg0h-M5P_pG46-KieVhwAW8f4CV0uxI0YGWj3l8ztd7Lf_noyMpRRe1cPMBnuATLarbyVP_qjP4CkSx6Q_0-NjnKLXM1EK1rlqJg-MFoBrxRM3c59ypxnMuMekTfqQpTPwXMtlGEVxSHbgM2hGwnr6WsSuCmNuOer7HDHskpvukA/https%3A%2F%2Feur03.safelinks.protection.outlook.com%2F%3Furl%3Dhttps%253A%252F%252Fsecure-web.cisco.com%252F1lCp4yPi_vnD1_jDHjqMFFlcPKYlBFqg-Dg5nn2IbJrmH-ru9ATLtvP78StdglWKXz9Aoe2wkuZl3qWXtUlkzaPRyFPRyN-OfCI2-hTiv-3a80hsVv3ZU0pu_dlM-mpTTSgILcHfolVQX-BFcoHymMYXlXOQy8cuIBiKE6WSwV49wP1b003OX5yqYIduz3OPXVBpXTw-rkyt1naV1kGogJEbRs3QNS8dOuep2scIGZBsiejUTWQ70ctg1Q1aHGxp0cDDmqUiw5JiB4-w-xtjucQ%252Fhttps%25253A%25252F%25252Fmail.nmr.mgh.harvard.edu%25252Fmailman%25252Flistinfo%25252Ffreesurfer%26data%3D04%257C01%257Cjulia.griem%2540kcl.ac.uk%257C3f5e76a919d749738d0808d90018e9b9%257C8370cf1416f34c16b83c724071654356%257C0%257C0%257C637540929284868052%257CUnknown%257CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%253D%257C1000%26sdata%3DsjiRm%252BXt2HfAkzyBe4V8fixcdeM0VPaJDFm11VM2g%252BE%253D%26reserved%3D0> ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
Re: [Freesurfer] glmfit with aseg table: view surface overlay
The input is a table, so the output is a table (the sig.mgh file will just have the number of values that are in the table; this is not a "volume" you can look at in FV). Look in the sig.table.dat file for the sigs. On 4/13/2021 8:41 AM, Griem, Julia wrote: External Email - Use Caution External Email - Use Caution Dear Freesurfer team, I am comparing subcortical volumes (aseg table) in 2 groups with no other variables. Based on the GLM lecture slides from the course I attended, I ran the following glmfit command: mri_glmfit --table aseg_table.txt --fsgd FSGD/DemographicsSubcortBinaryGroup.fsgd --C Contrasts/HC-ASPD.mtx --C Contrasts/ASPD-HC.mtx --glmdir roi.binarygroup.glmdir It created all the correct output and subdirectories. I then used freeview -f $SUBJECTS_DIR/fsaverage/surf/lh.inflated:overlay=ASPD-HC/sig.mgh but an error saying "failed to load surface overlay" comes up and the terminal error message is "number of vertices in overlay data does not match with surface". I assume this is because I am trying to load a 'whole-brain' sig.mgh onto a lh surface? Do I have to change the glmfit command to specify hemispheres, and if so, how do I do that since I do not want to look at cortical but subcortical data? Or alternatively, can I load both hemispheres into freeview so the overlay sig.mgh fits? Thank you so much, Julia ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://secure-web.cisco.com/1lCp4yPi_vnD1_jDHjqMFFlcPKYlBFqg-Dg5nn2IbJrmH-ru9ATLtvP78StdglWKXz9Aoe2wkuZl3qWXtUlkzaPRyFPRyN-OfCI2-hTiv-3a80hsVv3ZU0pu_dlM-mpTTSgILcHfolVQX-BFcoHymMYXlXOQy8cuIBiKE6WSwV49wP1b003OX5yqYIduz3OPXVBpXTw-rkyt1naV1kGogJEbRs3QNS8dOuep2scIGZBsiejUTWQ70ctg1Q1aHGxp0cDDmqUiw5JiB4-w-xtjucQ/https%3A%2F%2Fmail.nmr.mgh.harvard.edu%2Fmailman%2Flistinfo%2Ffreesurfer ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
[Freesurfer] glmfit with aseg table: view surface overlay
External Email - Use Caution Dear Freesurfer team, I am comparing subcortical volumes (aseg table) in 2 groups with no other variables. Based on the GLM lecture slides from the course I attended, I ran the following glmfit command: mri_glmfit --table aseg_table.txt --fsgd FSGD/DemographicsSubcortBinaryGroup.fsgd --C Contrasts/HC-ASPD.mtx --C Contrasts/ASPD-HC.mtx --glmdir roi.binarygroup.glmdir It created all the correct output and subdirectories. I then used freeview -f $SUBJECTS_DIR/fsaverage/surf/lh.inflated:overlay=ASPD-HC/sig.mgh but an error saying "failed to load surface overlay" comes up and the terminal error message is "number of vertices in overlay data does not match with surface". I assume this is because I am trying to load a 'whole-brain' sig.mgh onto a lh surface? Do I have to change the glmfit command to specify hemispheres, and if so, how do I do that since I do not want to look at cortical but subcortical data? Or alternatively, can I load both hemispheres into freeview so the overlay sig.mgh fits? Thank you so much, Julia ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
