Re: [Freesurfer] Matrix ill-conditioned error for dods thickness mri_glmfit analysis

2017-08-16 Thread Douglas N Greve
posting the response to the list: The problem is that femalezero and 
femaletwo only have 4 and 3 subjects, respectively. As you are trying to 
estimate 4 slopes + an intercept, you will need at least 5 subjects per 
group.

doug




On 08/16/2017 03:27 PM, Sullivan, Danielle Renee wrote:
>
> Hi Doug,
>
>
> I am just following up to this email. If you prefer I can repost on 
> the FreeSurfer mail archive.
>
>
> Thanks,
>
> Danielle R. Sullivan (Miller), Ph.D.
>
> National Center for PTSD
>
> VA Boston Healthcare System Jamaica Plain
>
> Boston University School of Medicine
>
> OFFICE: (857) 364-4022
>
>
> 
> *From:* Sullivan, Danielle Renee
> *Sent:* Friday, August 11, 2017 1:33:03 PM
> *To:* Douglas N Greve
> *Subject:* Re: [Freesurfer] Matrix ill-conditioned error for dods 
> thickness mri_glmfit analysis
>
> I am still getting the same error. My apologies for the back and forth.
>
>
> Here are updated files (now there are 5 classes and contrast mtx has 
> been changed to reflect that).
>
>
> Danielle R. Sullivan (Miller), Ph.D.
>
> National Center for PTSD
>
> VA Boston Healthcare System Jamaica Plain
>
> Boston University School of Medicine
>
> OFFICE: (857) 364-4022
>
>
> 
> *From:* Douglas N Greve <gr...@nmr.mgh.harvard.edu>
> *Sent:* Friday, August 11, 2017 1:10:53 PM
> *To:* Sullivan, Danielle Renee
> *Subject:* Re: [Freesurfer] Matrix ill-conditioned error for dods 
> thickness mri_glmfit analysis
> remove the definition too
>
>
> On 08/11/2017 12:46 PM, Sullivan, Danielle Renee wrote:
> >
> > Hi,
> >
> >
> > Thank you for the tips. I tried doing both and still run into the
> > error (although now I get a warning that the class femaleone is
> > defined but not used). Attached is my new FSGD (now all variables
> > normalized via a Z-score), XG.mtx, Contrast, and Terminal output.
> >
> >
> > Thanks again for the help!
> >
> >
> > Danielle R. Sullivan (Miller), Ph.D.
> >
> > National Center for PTSD
> >
> > VA Boston Healthcare System Jamaica Plain
> >
> > Boston University School of Medicine
> >
> > OFFICE: (857) 364-4022
> >
> >
> > 
> > *From:* Douglas N Greve <gr...@nmr.mgh.harvard.edu>
> > *Sent:* Friday, August 11, 2017 11:12:21 AM
> > *To:* Sullivan, Danielle Renee
> > *Subject:* Re: [Freesurfer] Matrix ill-conditioned error for dods
> > thickness mri_glmfit analysis
> > there appear to be two problems. One is that you only have one instance
> > of femaleone. This might be ok if you had no covariates, but you can't
> > estimate a covariate slope from one data point. Second, the scale of
> > your covariates is much different, so you may need to normalize them. I
> > would try removing the lone femaleone and rerun
> >
> >
> > On 08/11/2017 10:35 AM, Sullivan, Danielle Renee wrote:
> > >
> > > Hi Doug,
> > >
> > >
> > > Thanks for your response and help.
> > >
> > >
> > > Attached is the matrix and the terminal output. Please let me know if
> > > you need anything else.
> > >
> > >
> > >
> > > Danielle R. Sullivan (Miller), Ph.D.
> > >
> > > National Center for PTSD
> > >
> > > VA Boston Healthcare System Jamaica Plain
> > >
> > > Boston University School of Medicine
> > >
> > > OFFICE: (857) 364-4022
> > >
> > >
> > > can you send the design matrix Xg.mtx and the full terminal output?
> > > On 08/08/2017 09:13 AM, Sullivan, Danielle Renee wrote:
> > > >
> > > > Hi FreeSurfer Experts,
> > > >
> > > >
> > > > I am trying to run a dods cortical thickness analysis examining the
> > > > interaction of group on a continuous variable while controlling for
> > > > other variables.
> > > >
> > > >
> > > > Specifically, I have 6 classes (gender by three groups [zero, one,
> > > > two]). I am interested in the group interaction with a continuous
> > > > variable while controlling for gender (among 4 other continuous
> > > > variables). I set up my fsgd file (see attached) and contrast file
> > (36
> > > > dimensions- F test):
> > > >
> >

Re: [Freesurfer] Matrix ill-conditioned error for dods thickness mri_glmfit analysis

2017-08-09 Thread Douglas N Greve
can you send the design matrix Xg.mtx and the full terminal output?


On 08/08/2017 09:13 AM, Sullivan, Danielle Renee wrote:
>
> Hi FreeSurfer Experts,
>
>
> I am trying to run a dods cortical thickness analysis examining the 
> interaction of group on a continuous variable while controlling for 
> other variables.
>
>
> Specifically, I have 6 classes (gender by three groups [zero, one, 
> two]). I am interested in the group interaction with a continuous 
> variable while controlling for gender (among 4 other continuous 
> variables). I set up my fsgd file (see attached) and contrast file (36 
> dimensions- F test):
>
> *0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 [1 -1 0 1 
> -1 0]/2*
>
> *0 0 0 0 0 0 **0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 [1 -0 
> -1 1  0 -1]/2*
>
>
> and ran this command:
>
>
> *mri_glmfit --y 
> **lh.age.PC1_forFSGD.PC2_forFSGD.PC3_forFSGD.GGT_forFSGD.rs8042149_A_bysex_forfsgd.thickness.10.mgh
>  
> --fsgd 
> fsgd.age.PC1_forFSGD.PC2_forFSGD.PC3_forFSGD.GGT_forFSGD.rs8042149_A_bysex_forFSGD.152.thickness.08082017.txt
>  
> dods --C 
> /GLM_CONTRAST.MAT/Contrast.Male.Female.Age.PC1.PC2.PC3.GGT.RORAbySex.mtx 
> --surf fsaverage lh --cortex --glmdir 
> lh.age.PC1_forFSGD.PC2_forFSGD.PC3_forFSGD.GGT_forFSGD.rs8042149_A_bysex_forfsgd.thickness.10.glmdir
>  
> *
>
>
> But keep getting this error:
>
>
> *ERROR: matrix is ill-conditioned or badly scaled, condno = 1e+08*
>
> **
>
> *Possible problem with experimental design:*
>
> *Check for duplicate entries and/or lack of range of*
>
> *continuous variables within a class.*
>
> *If you seek help with this problem, make sure to send:*
>
> *  1. Your command line:*
>
> *mri_glmfit --y 
> lh.age.PC1_forFSGD.PC2_forFSGD.PC3_forFSGD.GGT_forFSGD.rs8042149_A_bysex_forfsgd.thickness.10.mgh
>  
> --fsgd 
> fsgd.age.PC1_forFSGD.PC2_forFSGD.PC3_forFSGD.GGT_forFSGD.rs8042149_A_bysex_forFSGD.152.thickness.08082017.txt
>  
> dods --C 
> /Volumes/VA_Imaging/Projects/salat/2389_TRT/dmiller/GGT_CRP/GLM_CONTRAST.MAT/Contrast.Male.Female.Age.PC1.PC2.PC3.GGT.RORAbySex.mtx
>  
> --surf fsaverage lh --cortex --glmdir 
> lh.age.PC1_forFSGD.PC2_forFSGD.PC3_forFSGD.GGT_forFSGD.rs8042149_A_bysex_forfsgd.thickness.10.glmdir
>  
> *
>
> *  2. The FSGD file (if using one)*
>
> *  3. And the design matrix above*
>
>
>
> Any idea on what I am doing wrong?
>
>
> Danielle R. Miller, Ph.D.
>
> National Center for PTSD
>
> VA Boston Healthcare System Jamaica Plain
>
> Boston University School of Medicine
>
> OFFICE: (857) 364-4022
>
>
>
>
> ___
> Freesurfer mailing list
> Freesurfer@nmr.mgh.harvard.edu
> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

-- 
Douglas N. Greve, Ph.D.
MGH-NMR Center
gr...@nmr.mgh.harvard.edu
Phone Number: 617-724-2358
Fax: 617-726-7422

Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2
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


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.