<div dir="ltr">thanks a lot. it clears the matter<div>Yuval<br><br><div class="gmail_quote">On 19 March 2012 18:59, Gio Piantoni <span dir="ltr"><<a href="mailto:g.piantoni@nin.knaw.nl">g.piantoni@nin.knaw.nl</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Yuval,<br>
<br>
If you pass data as third argument, the function will substract the<br>
weights of rejected component(s) from your data. So, your data will<br>
have the initial rank minus the number of rejected component. No<br>
reconstruction occurs.<br>
<br>
See line 121 of ft_rejectcomponent or try the code below.<br>
<br>
nelec = 50;<br>
fs = 1000;<br>
ntrl = 5;<br>
<br>
data = [];<br>
data.label = cellfun(@(x)['E' num2str(x)], num2cell(1:nelec), 'uni', 0);<br>
data.fsample = fs;<br>
for i =1:ntrl<br>
  data.trial{i} = randn(nelec, fs);<br>
  data.time{i} = 1/fs:1/fs:1;<br>
end<br>
<br>
cfg = [];<br>
cfg.method = 'fastica';<br>
cfg.numcomponent = 10;<br>
comp = ft_componentanalysis(cfg, data);<br>
<br>
cfg = [];<br>
cfg.component = 2;<br>
data2 = ft_rejectcomponent(cfg, comp, data); % subtract comp from data<br>
data3 = ft_rejectcomponent(cfg, comp); % reconstruct from comp<br>
<br>
rank(data.trial{1}) % full rank<br>
rank(data2.trial{1}) % full rank - 1 rejected<br>
rank(data3.trial{1})  % rank of 'comp' - 1 rejected<br>
<div><div></div><div class="h5"><br>
On Mon, Mar 19, 2012 at 16:46, Yuval Harpaz <<a href="mailto:yuvharpaz@gmail.com">yuvharpaz@gmail.com</a>> wrote:<br>
> Dear ft_emaillist<br>
><br>
> I want to reject components obtained with with fastica.<br>
><br>
> if I only run ft_componentanalysis with cfg.numcomponent=20, creating only<br>
> 20 components, can I still reject a specific component by<br>
> ft_rejectcomponent(cfg,comp,data) or will it only reconstruct the date based<br>
> on the 20 components I calculated?<br>
><br>
> thanks,Yuval<br>
><br>
> --<br>
><br>
> Y.Harpaz<br>
><br>
> a link to the BIU MEG lab:<br>
> <a href="http://faculty.biu.ac.il/~goldsa/index.html" target="_blank">http://faculty.biu.ac.il/~goldsa/index.html</a><br>
><br>
><br>
</div></div>> _______________________________________________<br>
> fieldtrip mailing list<br>
> <a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
> <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
_______________________________________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
<a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div dir="ltr"><br>Y.Harpaz<br><br>a link to the BIU MEG lab:<br><a href="http://faculty.biu.ac.il/%7Egoldsa/index.html" target="_blank">http://faculty.biu.ac.il/~goldsa/index.html</a></div>

<br>
</div></div>