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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://sccn.ucsd.edu/pipermail/eeglablist/attachments/20130612/c0b55bed/attachment.html>