Question: Gwas: Adjust For Population Stratification In Impute2 Output
3
gravatar for w0nderingatwork
6.7 years ago by
Germany
w0nderingatwork40 wrote:

Dear all,

I'm working with a GWAS dataset and I am not so experienced yet with imputation of missing genotypes. But I managed to phase my dataset with SHAPEIT, impute with IMPUTE2 and did an association analysis with SNPTEST, so far so good. But: how do I now adjust for population stratification? I don't find any option in SNPTEST to adjust for Genomic Control, so I guess it can be done by adjusting for principal components as covariates. But how do I calculate principal components from IMPUTE2 output (file types: .gen; .gen_info)??

Thanks a lot! Frank

gwas • 3.4k views
ADD COMMENTlink modified 6.7 years ago by kumar.vinod81280 • written 6.7 years ago by w0nderingatwork40

You can transform IMPUTE2 output with gtool. http://www.well.ox.ac.uk/~cfreeman/software/gwas/gtool.html

"GTOOL is a program for transforming sets of genotype data for use with the programs SNPTEST and IMPUTE. "

ADD REPLYlink written 6.7 years ago by Maxime Lamontagne2.2k

Hi Maxime, thanks a lot, that is a good option. The only drawback is that I loose the probability values for the imputed genotypes, because GTOOL just applies a propability threshold and will then convert the imputed genotypes to plink format such as regularly typed genotypes. That means that in subsequent association analysis in plink the uncertainty of the imputed genotypes is not considered any more. Best regards, Frank

ADD REPLYlink written 6.7 years ago by w0nderingatwork40

Calculate principle components on genotyped set, then use it for imputed set, see if it helps to reduce the inflation.

ADD REPLYlink written 6.7 years ago by zx87549.7k
1

thank you! I will try that. Best regards

ADD REPLYlink written 6.7 years ago by w0nderingatwork40
1
gravatar for kumar.vinod81
6.7 years ago by
kumar.vinod81280
New Delhi
kumar.vinod81280 wrote:

Hi, I dont know much about IMPUTE2 for population stratification but after converting your data in HapMap, Plink or flapjack format, you can import data in either TASSEL or GAPIT or in STRUCTURE, where you can run them for PCA analysis or structure analysis. Thanks,

Vinod

ADD COMMENTlink written 6.7 years ago by kumar.vinod81280
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: 2030 users visited in the last hour