Dear Experts,
during an fsfast analysis an error occured after selxavg3-sess was executed. The corresponding commands were:
mkanalysis-sess -fsd bold -surface self lh -fwhm 0 -event-related -paradigm sensory.par -nconditions 1 -spmhrf 0 -TR 2 -refeventdur 10 -nskip 8 -polyfit 2 -analysis analysis.lh -delay 0 -per-run -force -mcextreg
mkcontrast-sess -analysis analysis.lh -contrast stim-v-base -a 1
selxavg3-sess -s ./sessdir -analysis analysis.lh -no-preproc -force
The output of selxavg3-sess was:
>> >> >> >> >> >> >>
sxa3pwd =
/.../evaluation
>>
sxa3cmd =
/usr/local/freesurfer/fsfast/bin/selxavg3-sess -s ./sessdir -analysis analysis.lh -no-preproc -force
>> /usr/local/freesurfer/fsfast/toolbox/fast_selxavg3.m
>> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >>
#@# sessdir ###############################
/.../evaluation/sessdir
-------------------------
$Id: fast_selxavg3.m,v 1.100.2.2 2012/11/30 18:40:38 greve Exp $
/usr/local/freesurfer/fsfast/toolbox/fast_selxavg3.m
/usr/local/freesurfer/fsfast/toolbox/fast_ldanaflac.m
/usr/local/freesurfer/matlab/MRIread.m
-------------------------
outtop = /.../evaluation
Extension format = nii.gz
1 stim-v-base.mat
Excluding 8 points
nruns = 1
autostimdur =
outanadir = /.../evaluation/sessdir/bold/analysis.lh
Excluding 8 points
Found 145573/155188 (93.8) voxels in mask
Creating Design Matrix
... creation time = 0.008 sec
DoMCFit = 1
ntptot = 108, nX = 15, DOF = 93
Saving X matrix to /.../evaluation/sessdir/bold/analysis.lh/Xtmp.mat
XCond = 64.7067 (normalized)
Computing compensation for resdual AR1 bias
1 -0.5 -0.468819 (t=0.005249)
2 -0.25 -0.252194 (t=0.007053)
3 0 -0.0449419 (t=0.008126)
4 0.25 0.149516 (t=0.00935)
5 0.5 0.319324 (t=0.010795)
AR1 Correction M: 0.0749452 1.26122
Computing contrast matrices
OLS Beta Pass
run 1 t= 0.0
Global Mean
Global In-Mask Mean = (0.664951)
??? Error using ==> mrdivide
Matrix dimensions must agree.
Error in ==> fast_selxavg3 at 434
RescaleFactor = flac0.inorm/gmean;
>> ------------------------------------------
ERROR: fast_selxavg3() failed\n
Do you have an idea what went wrong?
Sincerely yours pfannmoe