Question: association analysis with stratification
gravatar for alessandrotestori7
8 days ago by
alessandrotestori7100 wrote:

Hello! I want to run a classical genome-wide association analysis (cases vs controls) in order to understand whether particular polymorphisms are relevant in a disease. The problem I have encountered is that cases come from a different population than controls. Is it possible to run an association analysis anyway (with a correction) and get reliable statistics or should I discard a lot of samples? Are there any detailed pipelines or packages I can use? Please let me know. Thanks!

ADD COMMENTlink modified 7 days ago by theobroma22710 • written 8 days ago by alessandrotestori7100
gravatar for theobroma22
7 days ago by
theobroma22710 wrote:

You can use the mantelhaen.test function in the stats R package, which conveniently is loaded in the library directory when you installed R. Alternatively, look at the epiR package, although I've never used it.

ADD COMMENTlink written 7 days ago by theobroma22710

Interesting, can you elaborate how it deals with population stratification? I considered studies with different control population vs patients dead upon arrival.

ADD REPLYlink written 7 days ago by WouterDeCoster19k

It would've been correct to consider populations with vastly different polymorphisms as dead upon arrival...but, in statistical modelling you could have also considered 1) having two controls, or 2) contrasting population cases with non-population controls, or 3) stratifying them separately in the contingency table which can help with the analysis like detecting any confounding variables after you obtain the point estimates and chi-squared values.

ADD REPLYlink written 7 days ago by theobroma22710

Thanks for your reply, though I'm not sure that is exactly what I need. However, please do let me know if I'm wrong about it. Indeed, I was rather looking for a pipeline or a package which can adjust p-values from an association study in case the controls and the cases have different allele distributions. I think several softwares/tools are available (e.g.: EIGENSTRAT, Price el al 2006, GenABEL?), but I was unable to follow a pipeline or find a recent package.

ADD REPLYlink written 6 days ago by alessandrotestori7100

If you need to adjust p-values you can get the q-values and local FDR values, using the qvalue package available on Bioconductor.

ADD REPLYlink written 6 days ago by theobroma22710

Now I see your study more clearly, and think the non-population controls were added to the data set as variation, or noise.

ADD REPLYlink written 6 days ago by theobroma22710
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: 811 users visited in the last hour