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 ...
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.