Question: GWASTools: assocTestRegression, forcing the analysis to use reference allele
0
gravatar for Guest User
5.2 years ago by
Guest User12k
Guest User12k wrote:
In the package GWASTools, the function assocTestRegression performs an association test based on the minor allele. However this makes comparison across different datasets complicated because the minor allele in one population may be different from another. Is there a way to perform this test using the reference allele instead? Thanks, Danica -- output of sessionInfo(): -- Sent via the guest posting facility at bioconductor.org.
gwastools • 554 views
ADD COMMENTlink modified 5.2 years ago by Stephanie M. Gogarten720 • written 5.2 years ago by Guest User12k
Answer: GWASTools: assocTestRegression, forcing the analysis to use reference allele
0
gravatar for Stephanie M. Gogarten
5.2 years ago by
University of Washington
Stephanie M. Gogarten720 wrote:
Hi Danica, There is not a way to change which allele is being tested in the code, but it is fairly straightforward to adjust the output to compare with a different effect allele. The output contains a "minor.allele" column which will tell you whether allele A or B was the minor (tested) allele in your sample. If the minor allele was B but the reference allele was A, then you would reverse the sign of "Stat" (beta) and take the inverse of the odds ratio to get the effect in terms of allele A. Defining which allele to test is a somewhat arbitrary choice, and since the definition of A and B alleles for GWAS data is usually determined by the array manufacturer, we chose to use the minor allele for greater consistency. (This is the same convention followed by PLINK, a popular tool for association testing.) Testing the reference allele would require additional information not contained in the dataset itself (which organism? which version of the reference sequence?). Stephanie On 8/12/14 2:37 AM, Danica [guest] wrote: > In the package GWASTools, the function assocTestRegression performs an association test based on the minor allele. However this makes comparison across different datasets complicated because the minor allele in one population may be different from another. > > Is there a way to perform this test using the reference allele instead? > > Thanks, > Danica > > -- output of sessionInfo(): > > > > -- > Sent via the guest posting facility at bioconductor.org. >
ADD COMMENTlink written 5.2 years ago by Stephanie M. Gogarten720
Hi Stephanie, Thank you for your advice. Indeed as you have said it is easy enough to alter the results to suit the allele of interest. However I was asking the question mostly in the context of a dataset involving thousands of SNPs where doing this alteration would require a bit of extra coding. Regards, Danica [[alternative HTML version deleted]]
ADD REPLYlink written 5.2 years ago by Danica Fabrigar10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 247 users visited in the last hour