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