Skip to main content

Table 2 Initial parameter vector (α1, α2, α3) and their estimation after several iterations of the EM algorithm

From: PSE-HMM: genome-wide CNV detection from NGS data using an HMM with Position-Specific Emission probabilities

Initial

Estimated

\( \begin{array}{cc}& \begin{array}{ccc}{\alpha}_1& {\alpha}_2& {\alpha}_3\end{array}\\ {}\begin{array}{c} diploid\\ {}\hfill 1\kern0.5em \mathrm{copy}\hfill \\ {}\hfill 0\kern0.5em \mathrm{copy}\hfill \\ {}\hfill \mathrm{copy}>2\hfill \end{array}& \left[\begin{array}{ccc}\hfill 1\hfill & \hfill 0\hfill & \hfill 0\hfill \\ {}\hfill 0.5\hfill & \hfill 0.5\hfill & \hfill 0\hfill \\ {}\hfill 0\hfill & \hfill 1\hfill & \hfill 0\hfill \\ {}\hfill 0.5\hfill & \hfill 0\hfill & \hfill 0.5\hfill \end{array}\right]\end{array} \)

\( \begin{array}{cc}& \begin{array}{ccc}{\alpha}_1& {\alpha}_2& {\alpha}_3\end{array}\\ {}\begin{array}{c} diploid\\ {}\hfill 1\kern0.5em \mathrm{copy}\hfill \\ {}\hfill 0\kern0.5em \mathrm{copy}\hfill \\ {}\hfill \mathrm{copy}>2\hfill \end{array}& \left[\begin{array}{ccc}\hfill 0.998\hfill & \hfill \approx 0\hfill & \hfill \approx 0\hfill \\ {}\hfill 0.5\hfill & \hfill 0.5\hfill & \hfill 0\hfill \\ {}\hfill 0\hfill & \hfill 1\hfill & \hfill 0\hfill \\ {}\hfill 0.47\hfill & \hfill 0\hfill & \hfill 0.53\hfill \end{array}\right]\end{array} \)