08.02.2013 Views

New Statistical Algorithms for the Analysis of Mass - FU Berlin, FB MI ...

New Statistical Algorithms for the Analysis of Mass - FU Berlin, FB MI ...

New Statistical Algorithms for the Analysis of Mass - FU Berlin, FB MI ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

48 CHAPTER 3. MATHEMATICAL MODELING AND ALGORITHMS<br />

Figure 3.6.15: This shows a spectrum alignment using key reference key peaks.<br />

The top left picture show <strong>the</strong> reference key peaks and <strong>the</strong>ir respective distances in xdirection.<br />

The top right picture show a new spectrum to be aligned with <strong>the</strong> reference.<br />

The result <strong>of</strong> this alignment process is shown in <strong>the</strong> bottom picture.<br />

peaks. As can be seen, <strong>the</strong> part in range a1 is compressed while range b1 is<br />

stretched and partly filled with zero elements.<br />

O<strong>the</strong>r Alignment <strong>Algorithms</strong><br />

The problem <strong>of</strong> spectra alignment appears not only in mass spectrometry but<br />

<strong>for</strong> many o<strong>the</strong>r data, such as in NMR spectra alignment (see (Kim et al., 2006)<br />

<strong>for</strong> an Bayesian approach) or electrophoretic lane alignment (see (Aittokallio<br />

et al., 2001) <strong>for</strong> a dynamic time warping approach). In all <strong>of</strong> <strong>the</strong>se domains<br />

(and many more) different and similar algorithms have been developed (see<br />

e.g. (Mäkinen, 2007) <strong>for</strong> an example derived from sequence alignment). One<br />

large family <strong>of</strong> algorithms we have analyzed during this work have <strong>the</strong> following<br />

drawbacks:<br />

� They are operating globally, that is, local changes (during alignment)<br />

<strong>of</strong> an ordered spectrum x0 . . . xn at a particular position xi affects all<br />

positions xi+1 . . . xn.<br />

� They align exactly two spectra, that is, a new spectrum to a reference<br />

spectrum.<br />

� They can not incorporate any statistical in<strong>for</strong>mation.<br />

As we have shown above, our approach allows to circumvent <strong>the</strong>se shortcomings<br />

while using statistical in<strong>for</strong>mation which are incorporated naturally by<br />

aligning groups <strong>of</strong> peaks (masterpeaks) opposed to single peaks from two single<br />

spectra.

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

Saved successfully!

Ooh no, something went wrong!