<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
</head>
<body dir="auto">
<div>You could always do an svd in matlab?</div>
<div><br>
</div>
<div>Then you could do an accumulative sum of the eigenvalues on the diagonal of the output (S I think). The components are usually in rank order, so the first PCA component would explain the most variance.</div>
<div><br>
</div>
<div>Tyler</div>
<div><br>
<br>
</div>
<div><br>
On 13 Dec 2014, at 9:42 am, Joëlle Albrecht <<a href="mailto:joelle.albrecht@hotmail.com">joelle.albrecht@hotmail.com</a>> wrote:<br>
<br>
</div>
<blockquote type="cite">
<div><style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style>
<div dir="ltr">
<p class="MsoNormal">Dear EEGLAB-ers<o:p></o:p></p>
<p class="MsoNormal"><br>
</p>
<p class="MsoNormal">I’m trying to find out, how much variance of the original data the components (computed by ICA) explain, but I’m not sure how to do this. (My goal is something like: Component 1 explains 60% of the variance in the original data. etc.)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal" style="margin-bottom: 0.0001pt;">1) I used the following function:<o:p></o:p></p>
<p class="MsoNormal" style="margin-bottom: 0.0001pt;"><o:p> </o:p></p>
<p class="MsoNormal" style="margin-bottom: 0.0001pt;"><span style="font-size: 10pt; font-family: 'Courier New';">[pc,eigvec,sv] = runpca(EEG.data);</span><span style="font-size:12.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal" style="margin-bottom: 0.0001pt;"><span style="font-size: 10pt; font-family: 'Courier New';">[proj pvaf] = compvar(EEG.data, pc, eigvec, components);</span><span style="font-size:12.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">What exactly does the „percent variance accounted for (pvaf)“ mean? For example if I got a pvaf of about 60, does this mean that this component explains 60% of the original data (or is it more complicated)? These pvaf’s sum up to about
80, does this mean that 20% of the original data cannot be explained by all the components?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">2) I used another function, too (after completing the ICA):<o:p></o:p></p>
<p class="MsoNormal" style="margin-bottom: 0.0001pt;"><span style="font-size:
10.0pt;font-family:"Courier New"">[pvaf,pvafs,vars]=eeg_pvaf(EEG,component);</span><span style="font-size:12.0pt;font-family:"Courier New""><o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Again: What exactly do the pvaf’s mean? Plus, these pvaf‘s are very different from the pvaf’s I obtained with the function above. What is the difference? These pvaf’s add to over 400 (whereas the ones from above sum up to only 80 which
seems to be more reasonable to me), does this mean that the independent components overlap? The pvaf’s resulting from this function (but not the ones from the one above) are not monotonously decreasing – how are they sorted?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thank you very much for the support! <o:p></o:p></p>
<p class="MsoNormal">Joëlle </p>
</div>
</div>
</blockquote>
<blockquote type="cite">
<div><span>_______________________________________________</span><br>
<span>Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a></span><br>
<span>To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu">
eeglablist-unsubscribe@sccn.ucsd.edu</a></span><br>
<span>For digest mode, send an email with the subject "set digest mime" to <a href="mailto:eeglablist-request@sccn.ucsd.edu">
eeglablist-request@sccn.ucsd.edu</a></span></div>
</blockquote>
</body>
</html>