Question: association analysis with stratification
1
gravatar for alessandrotestori7
6 months ago by
alessandrotestori7250 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 6 months ago by theobroma221.0k • written 6 months ago by alessandrotestori7250
1
gravatar for theobroma22
6 months ago by
theobroma221.0k
theobroma221.0k 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 6 months ago by theobroma221.0k
1

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 months ago by alessandrotestori7250

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 months ago by theobroma221.0k

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 months ago by theobroma221.0k

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 6 months ago by WouterDeCoster24k

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 6 months ago by theobroma221.0k
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: 640 users visited in the last hour