A simple pipeline that takes a freebayes vcf file of biallelic SNPs and an expression matrix in DESeq format and uses two different tests to detect eQTLs. First, we test for allelic imbalance in heterozygous individuals by fitting a beta binomial distribution for imbalance at each locus and using a likelihood ratio test for deviation from 5050. Next, we use limma/voom to log transform and calculate precision weights for expression values and use lm to test for associations between genotype and gene expression. P values are calculated using Stouffer's method. Fst and Fst outliers (OutFLANK), Differential expression between populations (DESeq2), and the proportion of population difference in expression explained by eQTLs can be calculated. In addition, functions for plotting expression associations, allelic imbalance, and allele distributions between populations are provided.
Package details 


Maintainer  
License  MIT License 
Version  1.0 
Package repository  View on GitHub 
Installation 
Install the latest version of this package by entering the following in R:

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.