[Eeglablist] cluster ERP projected on scalp electrodes

[Makoto Miyakoshi]

Dear Marco,

Try my backprojection toolbox. Unzip it and locate the folder under
/eeglab/plugins and relaunch eeglab. Let me know if it is useful. Currently
this is not distributed, but if you say it is useful I'll consider
distributing it.

Makoto


2013/6/13 Marco Bellami <mrcbellami at gmail.com>

>  Hi Makoto
>
> thank you so much for your suggestion.
>
> I took a look to  std_envtopo. However, can you help me in understanding
> how to use this plugin to extract a matrix with size  [number_of_channels,
> times],  having in each row the erp of the selected cluster ( e.g. cluster
> 2) projected  on one channel (e.g. on Fp1)?
>
> Thank you again
>
> Marco
>
>
>
> Il 12/06/2013 21:34, Makoto Miyakoshi ha scritto:
>
> Dear Marco,
>
>  Without reading your code, let me first ask if you have tried
> std_envtopo(). That plugin is available since version 11. After clustering,
> go to 'edit/plot cluster' and you'll find the button in the bottom in the
> GUI. It seems that is the function you need. If not, please let me know the
> difference.
>
>  Makoto
>
>
> 2013/6/10 Marco Bellami <mrcbellami at gmail.com>
>
>> Hi All,
>>
>> I performed an ICA decomposition on my raw EEG data (for each subject
>> and condition). Then, I created a STUDY and a design and I clustered ICs.
>> Now, I want to calculate, for each cluster and condition of my design,
>> the projection of the ERP of the cluster on the channels and I'd like to
>> know if my approach is correct.
>> First, I calculated the projection on the scalp of each  component of
>> each cluster (using the icaproj function). Then, I averaged the
>> projected ERPs between the components of the same cluster:
>>
>> for example,  let clust=1
>>
>>   for cond = 1:total_conditions % choose a condition (from
>> STUDY.condition)
>>
>>          % compute average cluster ERP projected on EEG channels
>>
>>          for ic = 1:length(STUDY.cluster(clust).allinds{cond}) % for
>> each component in the cluster
>>
>>              design_idx = STUDY.cluster(clust).setinds{cond}(1,ic);
>>              setidx =
>> STUDY.design(design_number).cell(design_idx).dataset;
>>              %select the component
>>              comp = STUDY.cluster(clust).allinds{cond}(1,ic);
>>              %load the dataset with the selected component
>>              [ALLEEG EEG CURRENTSET] = pop_newset(ALLEEG, EEG,
>> CURRENTSET,...
>>              'retrieve',setidx,'study',CURRENTSTUDY);
>>
>>              winvs =
>> ALLEEG(setidx).icawinv(:,:)*STUDY.cluster(clust).topopol(ic);
>>              % calculate the projection of the IC on all channels
>>
>> prjics(:,:,ic)=icaproj2(mean(ALLEEG(setidx).data(:,:,:),3),pinv(winvs),comp);
>>
>>          end;
>>          %average projection over trials
>>          prjclus=mean(prjics,3);
>>
>> Thanks,
>>
>> Marco
>>
>> _______________________________________________
>> Eeglablist page: http://sccn.ucsd.edu/eeglab/eeglabmail.html
>> To unsubscribe, send an empty email to
>> eeglablist-unsubscribe at sccn.ucsd.edu
>> For digest mode, send an email with the subject "set digest mime" to
>> eeglablist-request at sccn.ucsd.edu
>>
>
>
>
>  --
> Makoto Miyakoshi
> Swartz Center for Computational Neuroscience
> Institute for Neural Computation, University of California San Diego
>
>
>


-- 
Makoto Miyakoshi
Swartz Center for Computational Neuroscience
Institute for Neural Computation, University of California San Diego
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://sccn.ucsd.edu/pipermail/eeglablist/attachments/20130613/2587deee/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: backproj.zip
Type: application/zip
Size: 22545 bytes
Desc: not available
URL: <http://sccn.ucsd.edu/pipermail/eeglablist/attachments/20130613/2587deee/attachment.zip>