Charith Bhagya Karunarathna edited untitled.tex  almost 8 years ago

Commit id: 354248bb1588c17584e58dfa30b5233afc7d646e

deletions | additions      

       

\end{tabular}      \end{table} \section{Several popular methods}  Summary paragraph: We focus in this section on providing an overview of association tests, and how we used these association tests in our simulation study.  \subsection{Single-variant approach}  Under single-variant statistics, we evaluated Fisher's exact test of association. Fisher's exact test is a classical tool of studying association between genotype and disease traits with the use of contigency tables. We tested the null hypothesis of no association between rows (phenotypes), and columns (genotypes) of $2\times 3$ contigency table for each SNV. Each table contains the frequency of two homozygotes and the heterozygote among cases and controls. We then computed the P-value from each association.    \subsection{Pooled-variant method}  Variable threshold of Price et al.\cite{Price_2010} is based on the regression of phenotypes versus counts of mutations meeting the MAF threshold. Price et al. assume that the variants with MAF below some threshold, are more likely to be functional than variants with higher MAF. For each possible MAF threshold, a genotype score is computed based on given collapsing theme. The chosen MAF threshold maximizes the association signal and permutation testing is used to adjust for multiple thresholds. Price et al. found that VT approach had high power to detect the association between rare variant and disease trait, when effects are in one direction, in their simulations. Unlike VT-test, C-alpha test of Neale et al. is a non-burden single gene association method that helps improve power over burden test, especially when the effects of rare variants are in different directions. This C-alpha tests the variance of genetic effects under the assumption of rare variants observed in cases and controls are a mixture of deleterious, protective or neutral variants. We employed both VT-test and C-alpha test across the simulated region by using sliding windows of 20 SNVs overlapping by 5 SNVs.    \subsection{Joint-modeling method}   CAVIARBF\cite{Chen_2015} is a fine mapping method using marginal test statistics for the SNVs and their pairwise association in the Bayesian framework. This method approximates the Bayesian multivariate regression implemented in BIMBAM\cite{Servin_2005}. Chen et al. found that both CAVIARBF and BIMBAM have better performance than PAINTOR, and other methods. However, CAVIARBF is much faster than BIMBAM because it computes Bayse factors using only the SNVs in each causal model. These Bayes factors can be used to calculate the probability of SNVs being causal in the region.To compute the probability of SNVs being causal, set of models and their Bayes factors have to be considered. Let $p$ be the total number of SNVs in a candidate region, then the all possible number of causal models is $2^p$. Since it is difficult to compute all models for large $p$, this approach has a limitation on the number of causal variants in the model. So, this limitation reduces the number of models to evaluate in the model space, to $\sum_{i=0}^{l} \dbinom{p}{i}$, where $l$ is the number of causal SNVs in the model. Since there are 2630 SNVs in our data, to keep the computational load down, we considered $l=2$.