Hi all! I have a genotype dataset in plink format. Now I want to correct for population structure with PCA in association analysis. I split my dataset to training and testing datasets. I want to do the PCA only in the training dataset and use the training dataset as a reference panel to calculate the PCs in testing datasets (I don't want to combine them before PCA). Dose anyone know how to do that? Thank you!
•
written
by
Yong •
1 answer
Yes, you can multiply the genotypes of the test data (as alt-allele counts) with those from eigen-vector from the training-set-only PCA. That's very fast to do. I use it to derive per-person placements in a PCA plot of 1kgenomes individuals, like this for inspiration, and it takes literally just a split-second to get.
I've used the prcomp package in R for that calculation, you could follow the calculation in above code-link if you want (I think the original pca-call is here), but the principle is the same for any PCA-package.
Good luck!
Log in to answer this question.
-
Should genotype PCA in eQTL analysis be performed using all SNPs? or can be with a subset of SNPs?
written by maximal_life •
-
Can PCs from PCA of "Chip Genotypes" be used for post-imputation GWAS? [to adjust for population …
written by Pine •
-
Genetic Model when running PCA
written by Mengna Zhang •
-
Whether to include reference genome in PCA for association analysis?
written by Phoenix Mu •
-
PCA: summarize results from principal components analysis and look for population heterogeneity
written by F_cm_C •
-
GWAS; high genomic inflation factor
written by Rich •
-
PLINK Principal Components not adequately controlling for population stratification in linear regre…
written by dam4l
-
SmartPCA gives very different PCs for different chromosomes
written by janeshen91 •
-
Principle Component Analysis (PCA) for GWAS (SNP genotype data)
written by samorjoy •
-
Large Genomic Inflation Factor In Gwas
written by Liyf •
