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.

Chapter 15. Neurocomput<strong>in</strong>g, 69:1485-1501, 2006 217<br />

1<br />

0.5<br />

0<br />

−0.5<br />

Sources<br />

−1<br />

0 20 40<br />

4<br />

2<br />

0<br />

−2<br />

−4<br />

0 20 40<br />

2<br />

1<br />

0<br />

−1<br />

−2<br />

0 20 40<br />

(n)<br />

2<br />

1<br />

0<br />

−1<br />

1st step<br />

−2<br />

0 20 40<br />

4<br />

2<br />

0<br />

−2<br />

−4<br />

0 20 40<br />

2<br />

0<br />

−2<br />

−4<br />

0 20 40<br />

(n)<br />

2<br />

1<br />

0<br />

−1<br />

2st step<br />

−2<br />

0 20 40<br />

4<br />

2<br />

0<br />

−2<br />

−4<br />

0 20 40<br />

2<br />

1<br />

0<br />

−1<br />

−2<br />

0 20 40<br />

Fig. 5. Comparison of output signals result<strong>in</strong>g after the first step (second column)<br />

and the second step (last column) of dAMUSE.<br />

proton signal of the water solvent contam<strong>in</strong>at<strong>in</strong>g the prote<strong>in</strong> spectrum. This<br />

water artifact cannot be suppressed completely with technical means, hence<br />

it would be <strong>in</strong>terest<strong>in</strong>g to remove it dur<strong>in</strong>g the analysis of the spectra.<br />

BSS techniques have been shown to solve this separation problem [27,28]. BSS<br />

algorithms are based on an ICA [2] which extracts a set of underly<strong>in</strong>g <strong>in</strong>dependent<br />

source signals out of a set of measured signals without know<strong>in</strong>g how<br />

the mix<strong>in</strong>g process is carried out. We have used an algebraic algorithm [35,36]<br />

based on second order statistics us<strong>in</strong>g the time structure of the signals to<br />

separate this and related artifacts from the rema<strong>in</strong><strong>in</strong>g prote<strong>in</strong> spectrum. Unfortunately<br />

due to the statistical nature of the algorithm unwanted noise is<br />

<strong>in</strong>troduced <strong>in</strong>to the reconstructed spectrum as can be seen <strong>in</strong> figure 6. The water<br />

artifact removal is effected by a decomposition of a series of NMR spectra<br />

<strong>in</strong>to their uncorrelated spectral components apply<strong>in</strong>g a generalized eigendecomposition<br />

of a congruent matrix pencil [37]. The latter is formed with a<br />

correlation matrix of the signals and a correlation matrix with delayed or<br />

filtered signals [32]. Then we can detect and remove the components which<br />

conta<strong>in</strong> only a signal generated by the water and reconstruct the rema<strong>in</strong><strong>in</strong>g<br />

prote<strong>in</strong> spectrum from its ICs. But the latter now conta<strong>in</strong>s additional noise<br />

20<br />

(n)

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

Saved successfully!

Ooh no, something went wrong!