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

Commit id: 344da8e470f5794daf3a41468445ad5c62a85dd7

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 model the expected number of counts, $\lambda_{hi}$, 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}$, $T^{*}_{ij}$,  for each individual and target region. 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}