14.02.2013 Views

Mathematics in Independent Component Analysis

Mathematics in Independent Component Analysis

Mathematics in Independent Component Analysis

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Chapter 20. Signal Process<strong>in</strong>g 86(3):603-623, 2006 277<br />

4<br />

3.5<br />

3<br />

2.5<br />

2<br />

1.5<br />

1<br />

0.5<br />

0<br />

SCA<br />

sNMF *<br />

sNMF<br />

NMF *<br />

NMF<br />

sNMF *<br />

sNMF<br />

(a) matrix comparison by Amari <strong>in</strong>dex<br />

NMF *<br />

NMF<br />

JADE<br />

Measure mean [a.u.]<br />

1.5<br />

1<br />

0.5<br />

0<br />

Renyi<br />

QMI<br />

KLD<br />

MuIn<br />

RSDSTW<br />

Measure<br />

Xcor<br />

JADE<br />

Source<br />

SCA NMFNMF Channels<br />

sNMF *<br />

sNMF<br />

*<br />

Signal<br />

(b) comparison of the recovered sources<br />

Fig. 6. Inter-component performance <strong>in</strong>dex comparisons. (a) Comparison of the<br />

Amari <strong>in</strong>dex of matrices acquired from two different methods each <strong>in</strong> the case of<br />

synthetic s-EMG decomposition. As the comparison matrix is symmetric, half of<br />

its values have been omitted for clarity. (b) Mean of the <strong>in</strong>ter-components mutual<br />

<strong>in</strong>formation for each method estimated by different measures. For comparison, the<br />

measures correspond<strong>in</strong>g to the source signals (m<strong>in</strong>imal <strong>in</strong>dices) and to the channels<br />

of the mixed s-EMG (maximal <strong>in</strong>dices) are also shown.<br />

F<strong>in</strong>ally, we perform SCA us<strong>in</strong>g a generalized Hough transform [37]. Note that<br />

there are also other algorithms possible for such extraction, and model generalizations<br />

are possible, see for example [45]. After whiten<strong>in</strong>g and dimension<br />

reduction, we perform Hough hyperplane identification with b<strong>in</strong>-size 180 and<br />

manually identified the maxima <strong>in</strong> the obta<strong>in</strong>ed the Hough accumulator. The<br />

recovered mix<strong>in</strong>g matrix is aga<strong>in</strong> visualized <strong>in</strong> Fig. 4(f). Similar to the previous<br />

results, the three components are roughly most active at locations 2 to 3,<br />

4 and 5 respectively. Multiplication with the pseudo<strong>in</strong>verse of the recovered<br />

matrices from the above experiments shows that the result co<strong>in</strong>cides quite well<br />

with the matrices from NMF, but differs slightly more from the ICA-recovery.<br />

We calculate and compare the Amari <strong>in</strong>dex for each method and as shown <strong>in</strong><br />

Fig. 6(a), no major differences can be detected. A comparison of the recovered<br />

sources us<strong>in</strong>g the <strong>in</strong>dices from section 2.3 is given <strong>in</strong> Fig. 6(b), where also<br />

the <strong>in</strong>dices correspond<strong>in</strong>g to the orig<strong>in</strong>al source signals (m<strong>in</strong>imum mutual<br />

<strong>in</strong>formation values) and to the channels of the s-EMG (maximum values) have<br />

been added. All the methods separate the mixtures (improvement <strong>in</strong> terms of<br />

<strong>in</strong>dices) but the methods yield somewhat different result, with JADE giv<strong>in</strong>g<br />

quite different sources than the rest of the methods, and NMF and NMF∗<br />

perform<strong>in</strong>g rather similar. In terms of source <strong>in</strong>dependence, of course JADE<br />

scores best rat<strong>in</strong>gs <strong>in</strong> the <strong>in</strong>dices, as can be confirmed by calculat<strong>in</strong>g the Amari<br />

<strong>in</strong>dex of the recovered matrix with the orig<strong>in</strong>al source matrix:<br />

16

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!