Dear Freesurfer,
I have a question during qdec.
I tried to compare the surface area between two groups using qdec program. Strange to say, if I control some combination of confounding factors, then the error develops and qdec is closed.
When I control age, sex, HAMD (hamilton depression) score, and mean surface area, there was no error. When I control age and sex, there was no error. When I control only surface area, there was no error. However, when I control age, sex, and mean surface area, then error happens.
Could you let me how to manage this problem?
I am looking forward to receive your answer. Thank you.
Best, Seung-Gul
---------------------------------------------------------------------------------below is the messages containing error
SUBJECTS_DIR is '/mnt/hgfs/share/suicide'
lh-Avg-Intercept-area -----------------------
Does the average area differ from zero?
Nuisance factors: age area
1.000 1.000 1.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000;
lh-Diff-NS-S-Intercept-area -----------------------
Does the average area, accounting for sex, differ between NS and S?
Nuisance factors: age area
1.000 -1.000 1.000 -1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000;
lh-Diff-M-F-Intercept-area -----------------------
Does the average area, accounting for diagnosis, differ between M and F?
Nuisance factors: age area
1.000 1.000 -1.000 -1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000;
lh-X-diagnosis-sex-Intercept-area -----------------------
Is there a diagnosis--sex interaction in the mean area?
Nuisance factors: age area
1.000 -1.000 -1.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000;
ninputs = 38
Checking inputs
nframestot = 38
Allocing output
Done allocing
nframes = 38
Writing to /mnt/hgfs/share/suicide/qdec/area_LH_AA/y.mgh
gdfReadHeader: reading /mnt/hgfs/share/suicide/qdec/area_LH_AA/qdec.fsgd
INFO: DeMeanFlag keyword not found, DeMeaning will NOT be done.
Continuous Variable Means (all subjects)
0 age 41.5263 12.8201
1 area 163656 15524.4
Class Means of each Continuous Variable
1 diagnosisNS-sexM 36.2000 176040.8391
2 diagnosisS-sexM 40.6250 170897.0234
3 diagnosisNS-sexF 46.5556 151977.3351
4 diagnosisS-sexF 42.9091 156686.6932
INFO: gd2mtx_method is dods
Reading source surface /mnt/hgfs/share/suicide/fsaverage/surf/lh.white
Number of vertices 163842
Number of faces 327680
Total area 65416.648438
AvgVtxArea 0.399267
AvgVtxDist 0.721953
StdVtxDist 0.195470
$Id: mri_glmfit.c,v 1.196.2.8 2012/11/01 18:51:41 greve Exp $
cwd /home/sgk
cmdline mri_glmfit --y /mnt/hgfs/share/suicide/qdec/area_LH_AA/y.mgh --fsgd /mnt/hgfs/share/suicide/qdec/area_LH_AA/qdec.fsgd dods --glmdir /mnt/hgfs/share/suicide/qdec/area_LH_AA --surf fsaverage lh --label /mnt/hgfs/share/suicide/fsaverage/label/lh.aparc.label --C /mnt/hgfs/share/suicide/qdec/area_LH_AA/contrasts/lh-Avg-Intercept-area.mat --C /mnt/hgfs/share/suicide/qdec/area_LH_AA/contrasts/lh-Diff-NS-S-Intercept-area.mat --C /mnt/hgfs/share/suicide/qdec/area_LH_AA/contrasts/lh-Diff-M-F-Intercept-area.mat --C /mnt/hgfs/share/suicide/qdec/area_LH_AA/contrasts/lh-X-diagnosis-sex-Intercept-area.mat
sysname Linux
hostname localhost.localdomain
machine x86_64
user sgk
FixVertexAreaFlag = 1
UseMaskWithSmoothing 1
OneSampleGroupMean 0
y /mnt/hgfs/share/suicide/qdec/area_LH_AA/y.mgh
logyflag 0
usedti 0
FSGD /mnt/hgfs/share/suicide/qdec/area_LH_AA/qdec.fsgd
labelmask /mnt/hgfs/share/suicide/fsaverage/label/lh.aparc.label
maskinv 0
glmdir /mnt/hgfs/share/suicide/qdec/area_LH_AA
IllCondOK 0
ReScaleX 1
DoFFx 0
Creating output directory /mnt/hgfs/share/suicide/qdec/area_LH_AA
Loading y from /mnt/hgfs/share/suicide/qdec/area_LH_AA/y.mgh
INFO: gd2mtx_method is dods
Saving design matrix to /mnt/hgfs/share/suicide/qdec/area_LH_AA/Xg.dat
Normalized matrix condition is 2083.95
Matrix condition is 1e+08
Found 148151 points in label.
Pruning voxels by thr: 0.000000
Found 148025 voxels in mask
Saving mask to /mnt/hgfs/share/suicide/qdec/area_LH_AA/mask.mgh
Reshaping mriglm->mask...
search space = 73649.347769
DOF = 26
Starting fit and test
Fit completed in 0.0923167 minutes
Computing spatial AR1 on surface
Residual: ar1mn=0.997068, ar1std=0.001719, gstd=7.468720, fwhm=17.587492
Writing results
lh-Avg-Intercept-area
maxvox sig=3.56625 F=17.7049 at index 25775 0 0 seed=1506511386
lh-Diff-NS-S-Intercept-area
maxvox sig=3.42944 F=16.7076 at index 10791 0 0 seed=1506511386
lh-Diff-M-F-Intercept-area
maxvox sig=-3.12066 F=14.5487 at index 103184 0 0 seed=1506511386
lh-X-diagnosis-sex-Intercept-area
maxvox sig=3.53629 F=17.4843 at index 158713 0 0 seed=1506511386
mri_glmfit done
ninputs = 4
Checking inputs
nframestot = 4
Allocing output
Done allocing
nframes = 4
Writing to /mnt/hgfs/share/suicide/qdec/area_LH_AA/contrasts.sig.mgh
reading colortable from annotation file...
colortable with 36 entries read (originally /autofs/space/terrier_001/users/nicks/freesurfer/average/colortable_desikan_killiany.txt)
reading colortable from annotation file...
colortable with 36 entries read (originally /autofs/space/terrier_001/users/nicks/freesurfer/average/colortable_desikan_killiany.txt)
gdfReadHeader: reading /mnt/hgfs/share/suicide/qdec/area_LH_AA/y.fsgd
INFO: ignoring tag Creator
INFO: ignoring tag SUBJECTS_DIR
INFO: ignoring tag SynthSeed
INFO: NOT demeaning continuous variables
Continuous Variable Means (all subjects)
0 age 41.5263 12.8201
1 area 163656 15524.3
Class Means of each Continuous Variable
1 diagnosisNS-sexM 36.2000 176040.6000
2 diagnosisS-sexM 40.6250 170897.1250
3 diagnosisNS-sexF 46.5556 151977.4444
4 diagnosisS-sexF 42.9091 156686.6364
MatrixMultiply: m1 is null!
No such file or directory
[sgk@localhost ~]$