Hello FreeSurfers,
I'm trying to analyze some BOLD data where the subject moves through a
dozen or so states over the course of an hour. None of these repeat
but they are the same for all subjects. I have a .par file that
specifies the duration of each segment with an integer index for each
one. I don't need an HRF; a square wave is preferable. But that
doesn't seem to be an option in mkanalysis-sess so I used a gamma
function with an exponent of 0 and a long dispersion (~300 sec.) which
seems to approximate a square wave well enough. There is no "resting
state" to speak of, each state is meant to be contrasted against all
the others. It's a very simple design and yet difficult to fit into
an analysis. I thought I did it right (or well enough) but when I try
to run it I get a series of errors that Google is of no help in
elucidating:
??? Attempted to access a(1); index out of bounds because numel(a)=0.
Error in ==> fast_st2fir at 121
b(1) = a(1);
Error in ==> flac_desmat at 56
flacnew.ev(nthev).Xfir = fast_st2fir(st,flac.ntp,flac.TR,ev.psdwin,1);
Error in ==> flac_customize at 272
flacnew = flac_desmat(flacnew);
Error in ==> fast_selxavg3 at 97
flac0 = flac_customize(flac0);
>> ------------------------------------------
ERROR: fast_selxavg3() failed\n
Any help figuring this out would certainly be appreciated!
Thank you,
Alessandro D. Gagliardi
Ph.D. Candidate
Integrative Neuroscience Program
Rutgers University Mind Brain Analysis
adfgagliardi(a)gmail.com