Hi all,
I have some data that consists of 5 trials, each of which has 5 slices of functional data (265 time points per trial, TR=300ms). Since the experiment is quite basic (2 sec. visual flash, 3 times per trial), I assume that I would begin to see a bold response after averaging the 5 trials together (3dmean file_1+orig file_2+orig etc). However, I seem to be getting into trouble and think it may be due to the way in which I'm comparing my ideal.1D file to the real data. The ideal.1D file I created consists of 265 pts with '0's (rest) and '1's (activation), and then transformed it using:
waver -dt 0.3 -GAM -input ideal.1D > ideal_hrf.1D
Could I now baseline correct/detrend the data in each slice of each trial independently, and THEN average overall trials? If so, how would I do this?