<div dir="ltr">Hi Peter,<div><br></div><div style>Regarding your first question: after  maxfilter, you need to do a pca, simply because your rank of data has been reduced from 306 to about 64. You can do this by:</div><div style>
<br></div><div style><div>cfg        = [];</div><div>cfg.method = 'runica'; </div><div>n_comp = rank(squeeze(data.trial{1}) * squeeze(data.trial{1})');</div><div>cfg.runica.pca = n_comp;</div><div>cfg.runica.stop = 1e-7;</div>
<div>comp = ft_componentanalysis(cfg, data);<br></div><div><br></div><div style>Here, n_comp is the number of component and it is equal to the rank of data set. If you don't do this, the algorithm takes a long time and normally does not converge. Note that the data should be continuous and it is highly recommended to remove big and jump artefacts before ICA (for example using ft_databrowser). I found ICA using infomax 'runica' quite robust with high accuracy fro removing EOG and ECG.</div>
<div style><br></div><div style><br></div><div style>Regarding your second question: Using only 'MAG' to investigate the topograph of component is fine, and better than 'GRAD', but have a look at the components time-series as well. If you remove that component it will be removed from the whole data set.<br>
</div><div><br></div><div><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On 9 March 2013 01:56, Peter Goodin <span dir="ltr"><<a href="mailto:pgoodin@swin.edu.au" target="_blank">pgoodin@swin.edu.au</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div>
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">Hi Fieldtrip list, 
<div><br>
</div>
<div>Firstly I'd like to thank those who have asked questions using ICA on neuromag 306 data and the answers given. It's been extremely informative and taught me a lot about the general principles of ICA (it's no longer completely black magic, only partially...).<br>

<br>
I have two and a bit questions about using ICA on neuromag data, the first on the full rank assumption. After maxfiltering, the output data is rank deficient due to removal of the b-out components. I see that fastica is able to detect this and makes adjustments
 to the amounts of components extracted accordingly, but is doing this justified? What actually happens to the data if ICA is applied to a rank deficient data set? Are there any papers that examine this? </div>
<div><br>
</div>
<div>Secondly, I've been including all the meg channels (mags + grads) when running ICA (due to my interpretation of the conversation in <a href="http://mailman.science.ru.nl/pipermail/fieldtrip/2012-April/005016.html" target="_blank">http://mailman.science.ru.nl/pipermail/fieldtrip/2012-April/005016.html</a>) then examining the time series components for those
 related to artifact (specifically ECG and EOG), while using a magnetometer only layout to view component topographies as a backup. My question (well, really a confirmation) is, by removing the artifact components calculated using the mixed sensors, does this
 remove the associated components from all the sensor types? Viewing pre and post ICA data appears to confirm this, but expert opinion is always good. The bit of a question is - is it acceptable to use a magnetometer only layout to interpret topoplots of the
 components if the components are based on mixed sensors? I've been using it because the readily identifiable artifact components have a distribution I'd expect from ECG and EOG as does the associated time series but I'm wanting to make sure.  </div>

<div><br>
</div>
<div>Also, if this is the wrong forum to ask these questions can someone suggest one that might be a better fit? </div>
<div><br>
</div>
<div>Thanks again, </div>
<div><br>
</div>
<div>Peter.<br>
<div><br>
<div style="font-family:Tahoma;font-size:13px">
<div style="font-family:Tahoma;font-size:13px">
<div style="font-family:Tahoma;font-size:13px">
<div><font>__________________________</font></div>
<div><font>Peter Goodin,  </font>
<div>
<div>
<div>BSc (Hons), Ph.D Candidate.</div>
<div><br>
</div>
<div><font>Brain and Psychological Sciences Research Centre (BPsych)</font></div>
<div>Swinburne University, </div>
<div>Hawthorn, Vic, 3122</div>
</div>
<div><br>
</div>
<div><font>Monash Alfred Psychiatry Research Centre (MAPrc)</font></div>
<div><font>Level 1, Old Baker Building</font></div>
<div><font>Commercial Road</font></div>
<div><font>Melbourne, Vic, 3004</font><span style="font-size:small"> </span></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>

<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>Hamid R. Mohseni, PhD<br>
Post-Doctoral Research Assistant<br>Institute of Biomedical Engineering<br>University of Oxford, OX3 7DQ, UK<br>Tel: +44 (0) 1865 2 83826<br>
</div>