External Email - Use Caution        

Hi Sophie,

it's a bit hard to tell, but maybe one option would be to try fitting the model 
with just a random-intercept effect instead of both random-intercept and 
random-slope effects.

I.e. using [1] instead of [1 2] for the Zcols argument to lme_mass_fit_EMinit.

In my experience the random-intercept models are a bit easier to fit.

Please write if the error still persists.

Best regards,

Kersten

On Fri, 2023-11-24 at 01:42 +0000, Sophie Matis wrote:

CAUTION: This email originated from outside of DZNE. Do not click links or open 
attachments unless you recognize the sender and know the content is safe.
ACHTUNG: Dies ist eine externe E-Mail, bitte seien Sie vorsichtig beim 
Anklicken von Links oder Öffnen von Anhängen


        External Email - Use Caution
Hello FreeSurfer experts!

I have carried out some analysis of changes in cortical thickness over time 
using the LME Matlab (R23b) Toolbox distributed within Freesurfer (7.4.1) in 
line with the instructions on this page:  MailScanner has detected a possible 
fraud attempt from "secure-web.cisco.com" claiming to be 
https://secure-web.cisco.com/1vfqsNpqmJmN2xLvXPQjtL3uRg_7ssJ4vs2AWTjpgborp6wn-Y8UpD_80Jy1Uknaqn3X0QTcBVkjig91xtXifI2EfCt8L9h1JrxvBukl-VZ_CVZuFjf4FjZdcPACl_56MAixZDzSRinExy1pkDaKwNP-_kW969LdD8iuxQCjEGs99j32nInU8c_jeSpzkq0jqA28D3Ha39KMwm9FpWDRZI5pQWW-iIuHOWiYN__JGXTqtGljHGUc79ZjpT2tSmPZtHeadhtn0W6jNclzksxJLHY8-9JclXCWIkNaug8xwsWnUiuU3PmT8FRaRKU1YcSeN/https%3A%2F%2Fsurfer.nmr.mgh.harvard.edu%2Ffswiki%2FLinearMixedEffectsModels<https://secure-web.cisco.com/1T0jceJqjr8o0YltVDmmB4EAlEhKr0pMp8gIavJ6vaUX1HtK-Q7_-yis2I7b9vCr3VJvdjHFpm7iOSm-HwZDzLYf9Nkj3tT_o8O99hjM6rAdq1g-WORzpm0WPWMgQJgFAhXft6qOrg2xj11ey4u2NsJXZ3jTQi3YP8bPpvvJB7IqfGoNub-6dGd_s6cilB6iqLKwORhHOIfRB0ZG-bCXTKNUbxYdeLONr_Yrr1o35sEMNfnZLfwT9iBoaqpJ007DIJRYohd5R1rxIU-WwpMzx_LQw7DlgyJUcyfbsyDixY8O15o-b4J-7Z3i5E1XSMaoz/https%3A%2F%2Fsurfer.nmr.mgh.harvard.edu%2Ffswiki%2FLinearMixedEffectsModels>

At this stage I just want to look at changes over time within each of three 
groups. For two groups I have been able to proceed through the steps fine and 
have my significance files which look normal. For the last group I am running 
into an error on the “[lhTh0,lhRe] = lme_mass_fit_EMinit(X,[1 
2],Y,ni,lhcortex,3); “ step. MATLAB returns “Warning: Matrix is singular, close 
to singular or badly scaled” and lhTh0 and lhRe are created as matrices with 
only NaN values.

Has anyone run into this error before or have any tips for troubleshooting it?

I haven’t had any errors on the longitudinal processing steps and visual QC on 
my files looks normal. No errors on mris_preproc or mri_surf2surf.  An 
aparcstats2table for thickness doesn’t reveal any NaNs.

My design matrix is created with X = [ones(length(M),1) M] after following the 
suggested sortData code from the wiki
Contrast is [0 1]
This is my qdec file:

fsid    fsid-base       yearsfromBL
108_MR1       108   0.00
199_MR1       199   0.00
199_MR2       199   1.10
199_MR3       199   2.25
202_MR1       202   0.00
204_MR1       204   0.00
204_MR2       204   1.16
204_MR3       204   2.30
223_MR1       223   0.00
0821_MR1      0821  0.00
0931_MR1      0931  0.00
1010_MR1      1010  0.00
1010_MR2      1010  1.84
1031_MR1      1031  0.00
1031_MR2      1031  1.00
1031_MR3      1031  2.13
1064_MR1      1064  0.00
1064_MR2      1064  2.01
1125_MR1      1125  0.00
1134_MR1      1134  0.00
1134_MR2      1134  1.17
1134_MR3      1134  2.28
1233_MR1      1233  0.00
1233_MR2      1233  1.05
1373_MR1      1373  0.00
1439_MR1      1439  0.00
1439_MR2      1439  1.14
1439_MR3      1439  2.28
1475_MR1      1475  0.00
1475_MR2      1475  1.04
1475_MR3      1475  2.03
1522_MR1      1522  0.00
1522_MR2      1522  2.07
1541_MR1      1541  0.00
1541_MR2      1541  1.47
1551_MR1      1551  0.00
1551_MR2      1551  0.98
1551_MR3      1551  2.05
1572_MR1      1572  0.00
1572_MR2      1572  1.03
1625_MR1      1625  0.00
1625_MR2      1625  1.06
1637_MR1      1637  0.00
1641_MR1      1641  0.00
1656_MR1      1656  0.00
1656_MR2      1656  0.94
1656_MR3      1656  2.26
1755_MR1      1755  0.00
1755_MR2      1755  1.04
1755_MR3      1755  2.01
1829_MR1      1829  0.00
09154_MR1     09154 0.00
09154_MR2     09154 1.01
09154_MR3     09154 2.01
12110_MR1     12110 0.00
12110_MR2     12110 1.04


Any help appreciated!

Thank you,

Sophie Matis | Data Management & Research Officer
School of Health Sciences & Brain and Mind Centre
Faculty of Medicine and Health

The University of Sydney
Brain and Mind Centre, 94-100 Mallett Street, Camperdown (M02G) | The 
University of Sydney | NSW | 06
sophie.ma...@sydney.edu.au<mailto:shona.brin...@sydney.edu.au>  | 
sydney.edu.au<http://secure-web.cisco.com/1TuWK5WUev2Fb_3Y4GMAFPyrgRjl5K2o4Isb-heH5rrBUMMXjrcSBv64tyDOqZgQRYa26qVsbjDGLjAnnXaDJXUQdbW0LVvM_Y_Sd6Fda41zdGFyuVrM5eSK_pZ-24SCZ2gnso-c6sS_rF7WhUD1QQKKugvYO7mxXTOv2xrLMvIcbQaqJTphk2bzWXI753rXjMuj-06CGcUR9jT9D_6uhEPw5QNl_vnNJ8wlP3J4YCMwG0erQNjpVhM3fBFtW0a922mBpG6D88G1Y88s-AcNu7XMpIsZEV7DcDAiwr7ekJHGqy95lNBfIW2SNfaAsx5H1/http%3A%2F%2Fsydney.edu.au%2F>

CRICOS 00026A
This email plus any attachments to it are confidential. Any unauthorised use is 
strictly prohibited. If you receive this email in error, please delete it and 
any attachments.

Please think of our environment and only print this email if necessary.

_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu<mailto:Freesurfer@nmr.mgh.harvard.edu>
https://secure-web.cisco.com/1ydqXvTO5osDgovBFUJnQ5tUTYxdqDqwqjTMp3F_K1uxaHRM0YnNrANZx8GTtsd8_Djqt2ObvAZj7-LPywaA5don48aNL5TSzFYl3yhEYceczVAAhli-y7Dahxj2SGo6fzvY9Mh-jql__8Ia4L3dVPW-ZpMZlomFWyF_YZPwD9Q8AFgwCjbYbErcVEe5e5OMg0NGRk2t0CGllwpAktbHDrC2ravLyRSivl8q05NTvgNm5TkuBGgdpDUQT9LlFaQCelrm-HEMYUVkbX4bRZF6dhVymV-llZxi3lS4B2wNyjbADMoJnJLOiLZvqR4N8OSDf/https%3A%2F%2Fmail.nmr.mgh.harvard.edu%2Fmailman%2Flistinfo%2Ffreesurfer
Please note that this e-mail is not secure (encrypted).  If you do not wish to 
continue communication over unencrypted e-mail, please notify the sender of 
this message immediately.  Continuing to send or respond to e-mail after 
receiving this message means you understand and accept this risk and wish to 
continue to communicate over unencrypted e-mail.

_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
Please note that this e-mail is not secure (encrypted).  If you do not wish to 
continue communication over unencrypted e-mail, please notify the sender of 
this message immediately.  Continuing to send or respond to e-mail after 
receiving this message means you understand and accept this risk and wish to 
continue to communicate over unencrypted e-mail. 

Reply via email to