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.
