Dear all,
I am running an analysis with 2 groups with different stimulus types, levels and presentation modes. With a few exceptions, each subject had 4 stimulus-evoked BOLD runs, each with a particular set of stimulus type that resulted in a total of unique 10 EVs per subject (1st run=2 EVs, 2nd run= 2 EVs, 3rd run = 4 EVs, 4th run = 2 EVs).
Therefore, 1st level results outputs a total of 602 copes (296 for the 1st group, 306 copes for the 2nd group). Since the FEAT GUI becomes slow with such a huge amount of input copes for higher-level analysis, I ran a higher-level analysis to calculate group means for each group separately which outputs 10 copes per group (total of 20 copes).
I am now trying to run a higher level analysis using as input the 20 copes from the group means high-level analysis to perform F-tests, however FEAT is not happy with it. I am running the analysis using FSL 5.0.5 in a computing cluster using Red Hat Enterprise Linux Server release 5.6 (Tikanga), requesting 24 GB of RAM to be allocated to the FEAT command issued. When I run the analysis I get an error in the log (see below). I also attached here the design.fsf file I used to call FEAT in the command line, if that helps. Can you please help me finding the reasons for such error? Or maybe there is an alternative way to run this analysis that you might suggest?
Thank you.
Estephan Moana-Filho, D.D.S., M.S.
Ph.D. Candidate
Clinical Fellow, Oral & Maxillofacial Pain Program
Regional Center for Neurosensory Disorders
UNC School of Dentistry
[log in to unmask]
ln -s ../.files .files
Higher-level stats
cat ../design.lcon | awk '{ print }' > design.lcon
/nas02/apps/fsl-5.0.5/fsl/bin/flameo --cope=filtered_func_data --vc=var_filtered_func_data --dvc=tdof_filtered_func_data --mask=mask --ld=stats --dm=design.mat --cs=design.grp --tc=design.con --fc=design.fts --runmode=flame1
An exception has been thrown
nevs >= ntpts, Singular matrix.Trace: Gsmanager::flame_stage1; Gsmanager::run.
Log directory is: stats
Setting up:
ntptsing=20.000000
evs_group=1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
1.000000
WARNING: The passed in varcope file, var_filtered_func_data, contains voxels inside the mask with zero (or negative) values. These voxels will be excluded from the analysis.
nevs=20
ntpts=20
ngs=1
nvoxels=1732269
Running:
/bin/rm -f stats/zem* stats/zols* stats/mask* ; /bin/mv dof stats
Post-stats
child process exited abnormally
while executing
"if { [ catch {
for { set argindex 1 } { $argindex < $argc } { incr argindex 1 } {
switch -- [ lindex $argv $argindex ] {
-I {
incr argin..."
(file "/nas02/apps/fsl-5.0.5/fsl/bin/feat" line 119)
|