Graham McVicker edited Correcting read depth and GC.tex  almost 10 years ago

Commit id: 3990193bea1e28f2fceda0cf72119e69c61e79ad

deletions | additions      

       

For each target region $j$, we count the total number of reads $v_j$ and calculate the GC content $w_j$. Then, for each individual $i$, we fit coefficients $a_{0i}, a_{1i}, \ldots, b_{4i} $ that maximize the likelihood of the observed counts given the adjusted expected counts, $T^{*}_{ij}$:  \[  \Like\left( \textrm{L}\left(  D_i \left| a_{0i}, a_{1i}, \ldots, b_{4i} \right. \right) = \prod_j \Pr_{\mathrm{Pois}} \left(X_{ij} = x_{ij} \left| T^{*}_{ij} \right. \right)\\ \]  \[  T^{*}_{ij} = \exp\left(a_{0i} + a_{1i} w_j + a_{2i} w_j^2 + a_{3i} w_j^3 + a_{4i} w_j^4 \right) \left(b_{1i} v_j + b_{2i} v_j^2 + b_{3i} v_j^3 + b_{4i} v_j^4 \right)