Thank you, you are correct, I figured that my analysis was actually using a subset of the runs which doesn't contain all 69 conditions. I would still like to use this subset of runs so I think I will make new para files that assign less conditions (some conditions are stimuli of similar types that can probably be put together, like multiple different faces, multiple different bodies, etc.).
Cheers,
Gabe
Within a given session you must have all conditions represented. They all do not need to be in the same run, but they must all appear at least once across all runs. It looks like only 33 of the 69 are there. BTW, 69 is probably not a realistic number of conditions, way too big.
doug
Gabriel A. Tobon wrote:
Sorry about that, I attached the wrong error msg from while I was trying "selxavg3-sess". They are numbered 1-69.
Wed Jul 14 18:24:59 EDT 2010
INFO (kan_ecev_207ab): RunList = 008 009 010 011 012 013 014 015 016 017 018 019 020 021
ERROR: there are 33 non-null conditions in the
paradigm file. However, 69 non-null conditions were
specified when mkanalysis was run. These must be made
consistent in order to run selxavg-sess.
What is the difference between "selxavg-sess" and "selxavg3-sess"? The error from "selxavg3-sess" says:
??? Error using ==> svd
Input to SVD must not contain NaN or Inf.
Error in ==> cond at 40
s = svd(A);
Error in ==> fast_selxavg3 at 248
XCond = cond(XtX);
>> ------------------------------------------
ERROR: fast_selxavg3() failed\n
Thank you,
Gabe
<mailto:Freesurfer@nmr.mgh.harvard.edu>On Wed, Jul 14, 2010 at 6:21 PM, Douglas N Greve <greve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>> wrote:
69 conditions? Are they numbered 1-69? That output does not have
the right error msg in it.
Gabriel A. Tobon wrote:
Hello,
I have event related data with 69 different conditions, and
every condition does not appear in every run. This causes
selxavg-sess to throw an error. Part of the relevant debug
output is below. Why does every condition need to appear in
every run?
Thank you,
Gabe
echo "$sess " | & tee -a $LF
echo
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/data/kan_ecev_207ab
tee -a
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/analysis/log/selxavg3-sess-bold-mc_analyze-100714181048.log
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/data/kan_ecev_207ab
date | & tee -a $LF
date
tee -a
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/analysis/log/selxavg3-sess-bold-mc_analyze-100714181048.log
Wed Jul 14 18:10:48 EDT 2010
if ( ! -d $sess ) then
if ( ! -d
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/data/kan_ecev_207ab
) then
set funcdir = $sess/$fsd
set funcdir =
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/data/kan_ecev_207ab/bold
set anadir = $sess/$fsd/$analysisname
set anadir =
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/data/kan_ecev_207ab/bold/mc_analyze
if ( $#OutParent == 0 ) then
if ( 0 == 0 ) then
set anadir = $sess/$fsd/$analysisname
set anadir =
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/data/kan_ecev_207ab/bold/mc_analyze
else
else
echo anadir = $anadir | tee -a $LF
echo anadir =
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/data/kan_ecev_207ab/bold/mc_analyze
tee -a
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/analysis/log/selxavg3-sess-bold-mc_analyze-100714181048.log
anadir =
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/data/kan_ecev_207ab/bold/mc_analyze
if ( -e $anadir && $DoGLMFit ) then
if ( -e
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/data/kan_ecev_207ab/bold/mc_analyze
&& 1 ) then
echo "analysis $analysis alread exists for $sessbase" | tee -a $LF
echo analysis mc_analyze alread exists for kan_ecev_207ab
tee -a
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/analysis/log/selxavg3-sess-bold-mc_analyze-100714181048.log
analysis mc_analyze alread exists for kan_ecev_207ab
if ( $Skip == 0 && $Overwrite == 0 && ! $monly && ! $perrun &&
! $jkrun ) then
if ( 0 == 0 && 0 == 0 && ! 0 && ! 0 && ! 0 ) then
echo "ERROR: analysis $analysis already exists. You must " |
tee -a $LF
echo ERROR: analysis mc_analyze already exists. You must
tee -a
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/analysis/log/selxavg3-sess-bold-mc_analyze-100714181048.log
ERROR: analysis mc_analyze already exists. You must
echo " specify whether to -skip or -overwrite" | tee -a $LF
echo specify whether to -skip or -overwrite
tee -a
/mindhive/nklab/projects/eightcat-cluster/gtobon_analysis/8cat-ev2-final/analysis/log/selxavg3-sess-bold-mc_analyze-100714181048.log
specify whether to -skip or -overwrite
exit 1 ;
exit 1
------------------------------------------------------------------------
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edugreve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html>
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
-- Douglas N. Greve, Ph.D.
MGH-NMR Center
--
Douglas N. Greve, Ph.D.
MGH-NMR Center
greve@nmr.mgh.harvard.edu
Phone Number: 617-724-2358 Fax: 617-726-7422
Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
FileDrop: www.nmr.mgh.harvard.edu/facility/filedrop/index.html