Question: Genotype PCA from genotype matrix
3
gravatar for i.sudbery
5 weeks ago by
i.sudbery9.1k
Sheffield, UK
i.sudbery9.1k wrote:

If I have a individual x locus genotype matrix, what is the easiest way to calculate genotype principal components? Can I just use prcomp in R?

This lovely post Produce PCA bi-plot for 1000 Genomes Phase III - Version 2 goes through the process of generating a PCA from VCF files, but I wondered if I needed to use PLINK if I already had the matrix?

snp genotype pca • 161 views
ADD COMMENTlink modified 5 weeks ago • written 5 weeks ago by i.sudbery9.1k

I still don't know if PLINK, GEMMA, SNPRelate, etc. are doing anything other than generating a standard covariance matrix from the genotypes, and then just using that. Similar question here: https://bioinformatics.stackexchange.com/questions/4030/pca-on-genotype-matrix-with-multiple-alleles

ADD REPLYlink written 5 weeks ago by Kevin Blighe65k
1
gravatar for i.sudbery
5 weeks ago by
i.sudbery9.1k
Sheffield, UK
i.sudbery9.1k wrote:

I don't know if this is absolutely the most optimal way to do it, but we did straightforward prcomp on the genotype matrix and the first 2 PCs more or less recapitulated what we know about the population structure of the samples.

ADD COMMENTlink written 5 weeks ago by i.sudbery9.1k
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: 1114 users visited in the last hour