[Eeglablist] cluster ERP projected on scalp electrodes

[Marco Bellami]

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 
> <mailto: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
>     <mailto: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
>     <mailto:eeglablist-request at sccn.ucsd.edu>
>
>
>
>
> -- 
> 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/b0b5bee5/attachment.html>