Hi to All,
I have come across a problem using qdec, specified below. What I am trying to do is to look for an interaction between two discrete variables (diagnostik and PTSD). It is my understanding that this appears to be a common problem. However, is there a way around it, such as conceptualizing the design matrix differently? And how could that be done?
Thanks,
Hannah
[LOE5145MAC:~/Borderline_2/subjects] bruehh01% lh-Avg-Intercept-thickness -----------------------
Does the average thickness differ from zero?
1.000 1.000 1.000 1.000;
lh-Diff-NonPTSD-PTSD-Intercept-thickness -----------------------
Does the average thickness, accounting for diagnostik, differ between NonPTSD and PTSD?
1.000 -1.000 1.000 -1.000;
lh-Diff-BPD-KG-Intercept-thickness -----------------------
Does the average thickness, accounting for PTSD, differ between BPD and KG?
1.000 1.000 -1.000 -1.000;
lh-X-PTSD-diagnostik-Intercept-thickness -----------------------
Is there a PTSD--diagnostik interaction in the mean thickness?
1.000 -1.000 -1.000 1.000;
ninputs = 49
Checking inputs
nframestot = 49
Allocing output
nframes = 49
Writing to /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/y.mgh
gdfReadHeader: reading /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/qdec.fsgd
WARNING: gdfReadV1: class PTSDPTSD-diagnostikKG is defined but not used.
INFO: demeaning continous variables
Continuous Variable Means (all subjects)
Class Means of each Continuous Variable
1 PTSDNonPTSD-diagnostikBPD
2 PTSDPTSD-diagnostikBPD
3 PTSDNonPTSD-diagnostikKG
4 PTSDPTSD-diagnostikKG
INFO: gd2mtx_method is dods
Reading source surface /Users/bruehh01/Borderline_2/subjects/fsaverage/surf/lh.white
reading group avg surface area 822 cm^2 from file
Reading in average area /Users/bruehh01/Borderline_2/subjects/fsaverage/surf/lh.white.avg.area.mgh
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.138.2.15 2009/07/27 17:43:10 greve Exp $
cwd /Users/bruehh01/Borderline_2/subjects
cmdline mri_glmfit --y /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/y.mgh --fsgd /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/qdec.fsgd dods --glmdir /Users/bruehh01/Borderline_2/subjects/qdec/Untitled --surf fsaverage lh --label /Users/bruehh01/Borderline_2/subjects/fsaverage/label/lh.cortex.label --C /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/contrasts/lh-Avg-Intercept-thickness.mat --C /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/contrasts/lh-Diff-NonPTSD-PTSD-Intercept-thickness.mat --C /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/contrasts/lh-Diff-BPD-KG-Intercept-thickness.mat --C /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/contrasts/lh-X-PTSD-diagnostik-Intercept-thickness.mat
sysname Darwin
hostname LOE5145MAC.campus.fu-berlin.de
machine i386
user bruehh01
FixVertexAreaFlag = 1
UseMaskWithSmoothing 1
OneSampleGroupMean 0
y /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/y.mgh
logyflag 0
usedti 0
FSGD /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/qdec.fsgd
labelmask /Users/bruehh01/Borderline_2/subjects/fsaverage/label/lh.cortex.label
maskinv 0
glmdir /Users/bruehh01/Borderline_2/subjects/qdec/Untitled
IllCondOK 0
DoFFx 0
Creating output directory /Users/bruehh01/Borderline_2/subjects/qdec/Untitled
Loading y from /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/y.mgh
INFO: gd2mtx_method is dods
Saving design matrix to /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/Xg.dat
Matrix condition is 1e+08
ERROR: matrix is ill-conditioned or badly scaled, condno = 1e+08
Possible problem with experimental design:
Check for duplicate entries and/or lack of range of
continuous variables within a class.
1.000 0.000 0.000 0.000;
1.000 0.000 0.000 0.000;
0.000 1.000 0.000 0.000;
1.000 0.000 0.000 0.000;
1.000 0.000 0.000 0.000;
1.000 0.000 0.000 0.000;
0.000 1.000 0.000 0.000;
0.000 1.000 0.000 0.000;
1.000 0.000 0.000 0.000;
0.000 1.000 0.000 0.000;
0.000 1.000 0.000 0.000;
1.000 0.000 0.000 0.000;
1.000 0.000 0.000 0.000;
1.000 0.000 0.000 0.000;
0.000 1.000 0.000 0.000;
1.000 0.000 0.000 0.000;
0.000 1.000 0.000 0.000;
0.000 1.000 0.000 0.000;
1.000 0.000 0.000 0.000;
0.000 1.000 0.000 0.000;
0.000 1.000 0.000 0.000;
1.000 0.000 0.000 0.000;
1.000 0.000 0.000 0.000;
0.000 1.000 0.000 0.000;
1.000 0.000 0.000 0.000;
0.000 1.000 0.000 0.000;
1.000 0.000 0.000 0.000;
0.000 1.000 0.000 0.000;
0.000 1.000 0.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
0.000 0.000 1.000 0.000;
Error in Analyze: command failed: mri_glmfit --y /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/y.mgh --fsgd /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/qdec.fsgd dods --glmdir /Users/bruehh01/Borderline_2/subjects/qdec/Untitled --surf fsaverage lh --label /Users/bruehh01/Borderline_2/subjects/fsaverage/label/lh.cortex.label --C /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/contrasts/lh-Avg-Intercept-thickness.mat --C /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/contrasts/lh-Diff-NonPTSD-PTSD-Intercept-thickness.mat --C /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/contrasts/lh-Diff-BPD-KG-Intercept-thickness.mat --C /Users/bruehh01/Borderline_2/subjects/qdec/Untitled/contrasts/lh-X-PTSD-diagnostik-Intercept-thickness.mat
