Question: Differences In Significant Snps Between Plink And Pasw Using Logistic Regression
2
gravatar for Philipp Bayer
8.6 years ago by
Philipp Bayer6.9k
Australia/Perth/UWA
Philipp Bayer6.9k wrote:

Hello everyone,

for fun I checked whether there are any differences between PLINK and PASW/SPSS when doing a logistic regression (and standard parameters) on a binary trait. Using the same dataset, I get quite different significant SNPs - they sometimes overlap slightly, but sometimes the most significant SNPs in PLINK don't even appear in the model that PASW builds.

For testing I don't use any SNP-specific cleaning-steps etc., just wanted to check how both programs compare.

Does anybody have any explanation?

Thanks!

gwas plink • 2.6k views
ADD COMMENTlink written 8.6 years ago by Philipp Bayer6.9k
3
gravatar for ff.cc.cc
8.6 years ago by
ff.cc.cc1.3k
European Union
ff.cc.cc1.3k wrote:

Probably, PLINK and SPSS have two different default model approach. From the plink online help:

For the additive effects of SNPs, the direction of the regression coefficient represents the effect of each extra minor allele

Plink usually works on SNPs as a continous predictor, maybe SPSS manage your SNPs as a factor.

ADD COMMENTlink written 8.6 years ago by ff.cc.cc1.3k

So in doing a GWAS, I should "trust" Plink more than SPSS?

ADD REPLYlink written 8.6 years ago by Philipp Bayer6.9k
1

You could trust the plink logistic reg. approach if you believe in additive effects of the rare allele (more alleles -> more effect in almost linear way). If you believe there could be other possible explanation for your outcome (e.g. recessive, dominant, co-dominant...) you could treat SNPs as factors (each allele has its own effect on the outcome).

p.s. if you do not worry about covariates like age sex etc... but only about genetic component consider the --model option of plink

ADD REPLYlink modified 8.6 years ago • written 8.6 years ago by ff.cc.cc1.3k
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 2.3.0
Traffic: 952 users visited in the last hour
_