Eqtl Analysis In Plink
1
2
Entering edit mode
13.5 years ago
Rossella ▴ 370

Hi, I am trying to perform an eQTL analysis using Plink. I consider the gene expression of a particular gene as a quantitative trait and run Plink with the following options --linear --genotypic --mperm 100 --mperm-save-all

According to Plink the first line of the assoc.linear.perm file should contain the stats for the original dataset while the following lines contain stats for permuted phenotypes. When I compare the assoc.linear file that contains p-values and statistics and the first line of the perm file I get different statistics and I don't understand why.

Does somebody know why there is a difference between the stats?

Thanks a lot in advance

plink snp • 7.4k views
ADD COMMENT
0
Entering edit mode

I noticed this as well; the P values from association without permutation are from the Wald test (so say the docs). In principle the corresponding values in the permutation result should be the result of the same point-wise test, but in my results they differ slightly.

ADD REPLY
1
Entering edit mode
13.5 years ago

I would not use PLINK - from what I hear from my colleagues. Take a look at http://www.plosgenetics.org/article/info%3Adoi%2F10.1371%2Fjournal.pgen.1002078 for a recent article on identifying eQTL associations in primary human liver samples.

ADD COMMENT
1
Entering edit mode

My reason to use Plink resides in the fact that I want to reproduce results found in a paper and then apply the same technique (which is considered state of the art) to my dataset. The problem is that I cannot easily understand the output given by Plink. I could code everything myself but that would take a lot of time.

ADD REPLY

Login before adding your answer.

Traffic: 2439 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6