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.