14.02.2013 Views

Mathematics in Independent Component Analysis

Mathematics in Independent Component Analysis

Mathematics in Independent Component Analysis

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

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

In � this case we often omit the other variables and write f(x1, . . . , xn) =<br />

f1(x1), . . . , fn(xn) �<br />

or f = f1 × · · · × fn.<br />

Consider now the postnonl<strong>in</strong>ear Bl<strong>in</strong>d Source Separation model:<br />

X = f(AS)<br />

where aga<strong>in</strong> S is an <strong>in</strong>dependent random vector, A ∈ Gl(n) and f is a diagonal<br />

nonl<strong>in</strong>earity. We assume the components fi of f to be <strong>in</strong>jective analytic<br />

are analytic.<br />

functions with non-vanish<strong>in</strong>g derivatives. Then also the f −1<br />

i<br />

Def<strong>in</strong>ition 6 Let A ∈ Gl(n) be an <strong>in</strong>vertible matrix. Then A is said to<br />

be mix<strong>in</strong>g if A has at least two nonzero entries <strong>in</strong> each row 2 . And A =<br />

(aij)i,j=1...n is said to be absolutely degenerate if there are two columns l �= m<br />

such that a 2 il = λa 2 im for a λ �= 0 i.e. the normalized columns differ only by the<br />

signs of the entries.<br />

Postnonl<strong>in</strong>ear BSS is a generalization of l<strong>in</strong>ear BSS, so the <strong>in</strong>determ<strong>in</strong>acies<br />

of postnonl<strong>in</strong>ear ICA conta<strong>in</strong> at least the <strong>in</strong>determ<strong>in</strong>acies of l<strong>in</strong>ear BSS: A<br />

can only be reconstructed up to scal<strong>in</strong>g and permutation. Here of course additional<br />

<strong>in</strong>determ<strong>in</strong>acies come <strong>in</strong>to play because of translation: fi can only<br />

be recovered up to a constant. Also, if L ∈ Gl(n) is a scal<strong>in</strong>g matrix, then<br />

f(AS) = (f ◦L) �<br />

(L−1A)S �<br />

, so f and A can <strong>in</strong>terchange scal<strong>in</strong>g factors <strong>in</strong> each<br />

component. Another <strong>in</strong>determ<strong>in</strong>acy could occur if A is not mix<strong>in</strong>g, i.e. at least<br />

one observation xi conta<strong>in</strong>s only one source; <strong>in</strong> this case fi can obviously not<br />

be recovered. For example if A = I, then f(S) is already aga<strong>in</strong> <strong>in</strong>dependent,<br />

because <strong>in</strong>dependence is <strong>in</strong>variant under component-wise nonl<strong>in</strong>ear transformation;<br />

so f cannot be found us<strong>in</strong>g this method.<br />

A not so obvious <strong>in</strong>determ<strong>in</strong>acy occurs if A is absolutely degenerate. Then<br />

only the matrix A but not the nonl<strong>in</strong>earities can be recovered by look<strong>in</strong>g at<br />

the edges of the support of the fully-bounded random vector. For example<br />

1 1<br />

consider<br />

�<br />

the case n = 2, A = ( 2 −2 ) and the analytic function f(x1, x2) =<br />

x1 + 1<br />

2π s<strong>in</strong>(πx1), x2 + 1 x2 s<strong>in</strong>(π π 2 )�.<br />

Then A−1 ◦ f ◦ A maps [0, 1] 2 onto [0, 1] 2 .<br />

S<strong>in</strong>ce both components of f are <strong>in</strong>jective, we can verify this by look<strong>in</strong>g at the<br />

edges:<br />

2 A slightly more general def<strong>in</strong>ition of ’mix<strong>in</strong>g’ can be def<strong>in</strong>ed still guarantee<strong>in</strong>g<br />

identifiability of the sources; it is however omitted for the sake of simplicity.<br />

5

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

Saved successfully!

Ooh no, something went wrong!