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@gmail.com