dipole time course

[Floris de Lange]

Dear Fieldtrippers,

I have a fairly basic Fieldtrip question; I hope you can point me in the
right direction.
I'm running a dipole model on the early auditory response (60-90 ms after
the click). The dipoloe fit works very well. Just for completeness, these
are the lines of code I used:

cfg = [];
cfg.numdipoles = 2;
cfg.symmetry = 'y';
cfg.model = 'regional';
cfg.latency = [-0.44 -0.41];
cfg.vol = vol;
cfg.inwardshift  = 0;
source = dipolefitting(cfg,erf{1});

% visualize where the dipoles are
figure(1); clf; headmodelplot(cfg,erf{1})
hold on;
plot3([source.dip.pos(1,1)],[source.dip.pos(1,2)],[source.dip.pos(1,3)],'ro');
plot3([source.dip.pos(2,1)],[source.dip.pos(2,2)],[source.dip.pos(2,3)],'bo');

I have fitted the dipole over a 30 msec window. But now, I would like to
plot the response of these dipoles for the whole time window sampled.
In other words, I would like to use the coordinates obtained from the dipole
fitting procedure, and use them as 'virtual sensors' to reduce the
dimensionality of my data from 275 channel time courses to 2. I don't want
to use the dipole moments, because they're only defined within the 30 ms
that I used to calculate my dipole.
I suppose I have to project my data through the dipoles, but I'm not sure
how. Any help is greatly appreciated!

Best wishes,
Floris

--
Floris de Lange
http://www.florisdelange.com

----------------------------------
The aim of this list is to facilitate the discussion between users of the FieldTrip  toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis. See also http://listserv.surfnet.nl/archives/fieldtrip.html and http://www.ru.nl/neuroimaging/fieldtrip.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20090706/43bb8b8f/attachment.html>