[Marsbar-users] fitted event responses

[dsridhar <sri...@gm...>]

Hi,

I am trying to create fitted event responses using the following code: 

--------------------------------------------------------------------
% extract data
Y = get_marsy(roi, D, 'mean');

% MarsBaR estimation
E = estimate(D, Y);

% create fitted time courses
[this_tc dt] = event_fitted(E, event_spec, event_duration);
--------------------------------------------------------------------

This works great. However, when I extract raw time-courses with MarsBar and
average the raw response across events and subjects, the raw and fitted
waveforms seem to differ by quite a bit. 

I was wondering if it would it be possible for me to access/see the "raw"
waveform/timeseries that is being fitted (using the event_fitted function)?
I tried the following: 

------------------------
K>> SPM=des_struct(E) 

SPM = 

         nscan: [289 262]
          Sess: [1x2 struct]
           xBF: [1x1 struct]
           xGX: [1x1 struct]
            xX: [1x1 struct]
           xVi: [1x1 struct]
            xY: [1x1 struct]
         SPMid: 'SPM2: MarsBaR estimation. mardo_2 version 1.14'
            xM: [1x1 struct]
         xsDes: [1x1 struct]
         xMars: [1x1 struct]
         betas: [16x1 double]
    ResidualMS: 30.6604
          xCon: [1x1 struct]
         marsY: [1x1 marsy]

------------------------

But I couldn't find the raw waveform/timeseries (to-be-fitted) anywhere
inside this structure. Any help would be appreciated.

Thanks,
Sridhar
-- 
View this message in context: http://www.nabble.com/fitted-event-responses-tf4218667.html#a12001993
Sent from the MarsBaR mailing list archive at Nabble.com.