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

Commit id: 6721d632053ce0bd30bc3076b56584d554d7cf8d

deletions | additions      

       

\subsection{Correcting for GC content and other effects on expected read depth}  Since the number of mapped reads can differ between sequencing lanes and runs, we initially modelled modeled  the expected number of counts, $\lambda_{h,i}$, as a linear function of the total number of mapped reads for each individual, $T_i$. However, technical variation between experiments can change this relationship and reduce the power to detect true differences in read depths between samples or cause spurious associations. We directly model some known sources of technical variation and estimate adjusted total read depths, $T^{*}_{i,j}$, for each individual, $i$, and target region, $j$. This gives us a more accurate estimate of the expected number of reads and boosts our ability to detect true QTLs. \subsubsection{Adjusting total read depth}