Question: How Can I Calculate The Interaction Between Two Snp
gravatar for Liyf
9.3 years ago by
Liyf290 wrote:

I wonder that if I have two snps genotype data of 1000 cases and control, how can I calculate the p-value of interaction between two snp. I think there are so many models, so I am confused.

interaction association • 4.8k views
ADD COMMENTlink modified 9.3 years ago by • written 9.3 years ago by Liyf290
gravatar for
9.3 years ago by
European Union wrote:

The number of candidate tools to check epistatic patterns in GWAS is large, and only a few of them are able to really manage huge amount of data.

A couple of should-read articles by H.J. Cordell are:

i) Epistasis: what it means, what it doesn't mean, and statistical methods to detect it in humans. ii) Genome-wide association studies: Detecting gene-gene interactions that underlie human diseases.

A practical starting point is the Plink tool (Purcell). It performs a couple of epistatic test for snp pairs:

i) --logistic command will produce a p-value for the fitting of the beta parameter in a logistic model that check for allelic interaction
ii) --fast-epistasis command will generate a p-value related to a z-score test that in many situations reproduce the results of the logistic model, but performs much faster

Both options are well documented in the online manual.

There are many other methods and tools you could use or implement, but I suggest you to play with them in a second phase.

I cite only a few:

  • MDR & C. (Moore & others)
  • SNP-Harvester (Yang)
  • Mega SNP-Hunter (Wan)
  • BEAM (Zhang & Liu)
  • PIA (Goodman)
  • Boost (Wan)

Still there is not a standard for epistasis detection, you will face logistic models, classifiers, multi-stage filters, mutual information and so on, and the matter is very challenging

ADD COMMENTlink written 9.3 years ago by

Thank you very much. I have another question, if I detect some epistasis using plink and Beam, how can I know whether they are real or not? If I use another cohort to do replication, I still use plink and BEAM?

ADD REPLYlink written 9.3 years ago by Liyf290

Yes, if you have a good result, one of the best things to do (sometime difficult in GWAS) is to replicate it with the same method into an independent dataset

ADD REPLYlink written 9.3 years ago by
gravatar for Larry_Parnell
9.3 years ago by
Boston, MA USA
Larry_Parnell16k wrote:

You do need to make some assumptions - such as the recessive vs dominant action of the minor alleles involved in the epistasis. When someone invokes or applies epistasis, there is not a formulaic approach that dictates only one model to use to test for the association and epistasis. There is a bit of trial and error involved to determine which model to use.

ADD COMMENTlink written 9.3 years ago by Larry_Parnell16k

en, I think you are right. Thanks very much.

ADD REPLYlink written 9.3 years ago by Liyf290
Please log in to add an answer.


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