Dear Freesurfer,
We have a recurrent problem with an ill-conditioned matrix. Could someone shed some light on our problem.
Here is the command line:
mri_glmfit --y lh.ELS-CON.thickness.10.mgh --fsgd CAT-SP_June26.FSGF dods --C ELSvsCON_COV-Gen-Age-ICV.mtx --surf fsaverage lh --cortex --glmdir ELS_CON_June26_contrast_one
Below is the Design matrix and attached is the fsgd
Thanks very much,
Creating output directory ELS_CON_June26_contrast_one
Loading y from /Applications/freesurfer/subjects/SP-CAT_data/lh.ELS-CON.thickness.10.mgh
INFO: gd2mtx_method is dods
Saving design matrix to ELS_CON_June26_contrast_one/Xg.dat
Normalized matrix condition is 44113.9
Design matrix ------------------
0.000 0.000 0.000 1.000 0.000 0.000 0.000 29.600 0.000 0.000 0.000 1302835.625;
0.000 1.000 0.000 0.000 0.000 30.000 0.000 0.000 0.000 1785397.750 0.000 0.000;
0.000 0.000 0.000 1.000 0.000 0.000 0.000 29.930 0.000 0.000 0.000 1350884.250;
0.000 0.000 0.000 1.000 0.000 0.000 0.000 29.610 0.000 0.000 0.000 1492907.500;
0.000 0.000 0.000 1.000 0.000 0.000 0.000 29.380 0.000 0.000 0.000 1479810.250;
0.000 1.000 0.000 0.000 0.000 30.220 0.000 0.000 0.000 1408085.125 0.000 0.000;
0.000 0.000 0.000 1.000 0.000 0.000 0.000 30.470 0.000 0.000 0.000 1257495.875;
0.000 0.000 0.000 1.000 0.000 0.000 0.000 30.810 0.000 0.000 0.000 1221350.500;
0.000 0.000 0.000 1.000 0.000 0.000 0.000 29.740 0.000 0.000 0.000 1325784.500;
0.000 0.000 0.000 1.000 0.000 0.000 0.000 29.770 0.000 0.000 0.000 1428926.125;
0.000 0.000 0.000 1.000 0.000 0.000 0.000 30.310 0.000 0.000 0.000 1099130.750;
0.000 1.000 0.000 0.000 0.000 29.000 0.000 0.000 0.000 1766575.750 0.000 0.000;
0.000 0.000 0.000 1.000 0.000 0.000 0.000 29.050 0.000 0.000 0.000 1311269.250;
0.000 1.000 0.000 0.000 0.000 29.550 0.000 0.000 0.000 1749889.250 0.000 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 22.440 0.000 0.000 0.000 1617660.625 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 25.560 0.000 0.000 0.000 1749327.375 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 22.190 0.000 0.000 0.000 1532579.000 0.000;
1.000 0.000 0.000 0.000 25.530 0.000 0.000 0.000 2144050.750 0.000 0.000 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 25.950 0.000 0.000 0.000 1664001.125 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 22.790 0.000 0.000 0.000 1544209.625 0.000;
1.000 0.000 0.000 0.000 25.800 0.000 0.000 0.000 1718912.750 0.000 0.000 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 23.810 0.000 0.000 0.000 1624601.625 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 22.800 0.000 0.000 0.000 1706050.125 0.000;
1.000 0.000 0.000 0.000 21.250 0.000 0.000 0.000 1508681.750 0.000 0.000 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 25.720 0.000 0.000 0.000 1526436.125 0.000;
1.000 0.000 0.000 0.000 25.600 0.000 0.000 0.000 1558595.125 0.000 0.000 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 23.950 0.000 0.000 0.000 1523555.500 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 22.280 0.000 0.000 0.000 1436572.375 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 24.840 0.000 0.000 0.000 1391633.500 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 23.790 0.000 0.000 0.000 1484742.375 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 23.770 0.000 0.000 0.000 1523672.500 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 24.130 0.000 0.000 0.000 1616389.875 0.000;
0.000 0.000 1.000 0.000 0.000 0.000 25.790 0.000 0.000 0.000 1620330.625 0.000;
1.000 0.000 0.000 0.000 23.830 0.000 0.000 0.000 1785148.125 0.000 0.000 0.000;
1.000 0.000 0.000 0.000 23.790 0.000 0.000 0.000 1685555.500 0.000 0.000 0.000;
1.000 0.000 0.000 0.000 24.960 0.000 0.000 0.000 1791343.625 0.000 0.000 0.000;
--------------------------------
ERROR: matrix is ill-conditioned or badly scaled, condno = 44113.9
--------------------------------
Possible problem with experimental design:
Check for duplicate entries and/or lack of range of
continuous variables within a class.
If you seek help with this problem, make sure to send:
1. Your command line:
mri_glmfit --y lh.ELS-CON.thickness.10.mgh --fsgd CAT-SP_June26.FSGF dods --C ELSvsCON_COV-Gen-Age-ICV.mtx --surf fsaverage lh --cortex --glmdir ELS_CON_June26_contrast_one
2. The FSGD file (if using one)
3. And the design matrix above
--
<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>
Carl M. Anderson, Ph.D. (mail stop 204)
Assistant Professor of Psychiatry,
Harvard Medical School &
The Neuroimaging Center,
McLean Hospital 115 Mill St. Belmont, MA 02478
Ph: 617-855-2972; Fax: 617-855-2770
<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>