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.

100 Chapter 4. Neurocomput<strong>in</strong>g 64:223-234, 2005<br />

p<br />

1<br />

0<br />

-1<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.6<br />

0.2<br />

0.4<br />

0.2 0.4 0.6 0.8 1<br />

1<br />

1<br />

2<br />

3<br />

0.4 0.6<br />

4<br />

5<br />

0.6<br />

0.4<br />

0.4 0.6<br />

0.2 0.4 0.6 0.8 1<br />

Orig<strong>in</strong>al<br />

-1 0 1<br />

Relative volume error<br />

q<br />

2<br />

0<br />

-2<br />

0<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

6<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0<br />

�<br />

�<br />

− log vol(∆p,q)<br />

Unmix<strong>in</strong>g with p = 0.65, q = 0.5<br />

2<br />

0<br />

-2<br />

1<br />

0<br />

-1<br />

Unmix<strong>in</strong>g with p = 0.65, q = 0.65<br />

0<br />

Unmix<strong>in</strong>g with p = 0.5, q = 0.65<br />

-1 0 1<br />

Fig. 3. The top-left image graphs the separation error ∆p,q by measur<strong>in</strong>g the negative<br />

logarithm of its volume. The error ∆p,q denotes the difference set of po<strong>in</strong>ts from<br />

either the support of the recovered source distribution or a quadrangle with the<br />

same vertices. The other plots show the distributions of the recovered sources at<br />

some comb<strong>in</strong>ations of the parameters p, q of the nonl<strong>in</strong>earities. At p = 0.5, q = 0.5<br />

this is the orig<strong>in</strong>al source distribution. Here light grey areas represent ∆p,q and dark<br />

grey areas the <strong>in</strong>tersection of the support and the quadrangle.<br />

A simple estimator for mutual <strong>in</strong>formation based on histogram estimation of<br />

the entropy (with 10 b<strong>in</strong>s <strong>in</strong> each dimension) is used to check the <strong>in</strong>dependence<br />

of the recovered sources. A more elaborate histogram-based estimator<br />

by Moddemeijer [17] yields similar results. As shown <strong>in</strong> figure 2 the mutual<br />

<strong>in</strong>formation of the recovered sources is m<strong>in</strong>imal at the parameters p = q = 0.5,<br />

which correspond to the mix<strong>in</strong>g model. It can also be noticed that the m<strong>in</strong>imum<br />

is much less dist<strong>in</strong>ct <strong>in</strong> the second component. This <strong>in</strong>dicates that <strong>in</strong><br />

numerical application it should be easier to detect nonl<strong>in</strong>ear functions which<br />

are bounded.<br />

The second graph (figure 3) further illustrates that the criterion for the borders<br />

11

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

Saved successfully!

Ooh no, something went wrong!