Hi Cathy
You should put a comma between X and [1 2]
lhstats = lme_mass_fit_vw(X, [1 2],ni,lhcortex);
The model with a single random effect for the intercept term must always converge:
lhstats = lme_mass_fit_vw(X, [1],ni,lhcortex);
Can you tell me with words what the columns of your design matrix are? How many subjects and how many repeated measures for each subject do you have?
Best -Jorge
Sent from my iPad
On Jun 21, 2013, at 4:53, Catherine Bois C.Bois@sms.ed.ac.uk wrote:
Dear Jorge/freesurfer group,
I am using the scripts you sent me that do not require the newer version of matlab, and now using a computer that has the Statistics toolbox. I can get the scripts (with the suffix 1) to run with; lhstats = lme_mass_fit_vw(X[1 2],ni,lhcortex);) however it "fails to reach convergence" at most vertices. We have tried simplifying the model to include only one random effect at a time, however the problem persists. Our design matrix is as follows;
X = [ones(length(M),1) M(:,1) Mat(:,1) Mat(:,1).*M(:,1) Mat(:,2) Mat(:,2).*M(:,1) Mat(:,3) Mat(:,3).*M(:,1) M(:,3)-1 M(:,4)];
I have read on the mailing list that non-convergence at some vertices is normal, however what we are getting far exceeds 10%. Any help on this matter would be greatly appreciated!
Best Wishes,
Cathy
-- The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336.
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.