Graham McVicker edited Estimating overdispersion parameters.tex  over 9 years ago

Commit id: dbdf49f3db28ccbb17ef0ff0cea08de9e66489ab

deletions | additions      

       

\subsubsection{Beta-Binomial}  We calculate the genome-wide likelihood of $\Upsilon_i$ by taking the product of likelihoods from all target region SNPs that are heterozygous in individual $i$. We again assume there is no genetic effect, so $p$ = 0.5, and we find the maximum likelihood estimate of $\Upsilon_i$ numerically. numerically:  \[  \textrm{L}\left(\Upsilon_i \left| D \right. \right) = \prod_k \Pr_{\mathrm{BB-mix}} \left( Y = y_{ik} \left| n_{ik}, p = 0.5, \Upsilon_i \right. \right) \\