This is a test version of Biostars. For the public version, visit https://www.biostars.org.
PCA plot for TOP 500 genes + prcomp

Hi all,

Do you know how do I plot PCA for top500 genes by prcomp() the same as PCAplot in Desq2?

I would like to plot the PCA for top 500 genes for PCA3 and PCA4.

Thanks

rna-seq
• 6,362 views • link
written by mkh • mkh

IIRC prcomp gives you a rotation matrix, from which you can get genes that contribute the most to PC3 and PC4. You could then use that subset to plot.

• 1 views • link
written by Ram Ram

pcaExplorer is a cool package which can do it.

You can find more information here: pcaExplorer

• 1 views • link
updated by GenoMax GenoMax • written by mkh • mkh

As per Ram, you can easily do it with the rotation matrix from prcomp. For example:

pca <- prcomp(x)

# indices of top 500 genes on PC1 by absolute value
values1 <- order(abs(pca$rotation[,1]), decreasing=TRUE)[1:500]

# indices of top 500 genes on PC2 by absolute value
values2 <- order(abs(pca$rotation[,2]), decreasing=TRUE)[1:500]

et cetera

Note that DESeq2's plotPCA() function removes genes based on low variance prior to performing PCA. So, which do you want:

  1. remove genes based on low variance prior to performing PCA (via prcomp())?
  2. identify the top 500 genes contributing to variance in your dataset along PC3 and PC4 after you have performed PCA on your entire dataset (as per my code, above)?
• 1 views • link
written by Kevin Blighe Kevin Blighe

This is very helpful, thanks! I'm really new to prcomp and am having trouble figuring out what the x represents? Do you perhaps have a longer stretch of code I could look at? Thanks!

• 0 views • link
written by elisa.t.zhang • elisa.t.zhang

In this case, x represents your input data-matrix. By the way, my PCAtools package was accepted to Bioconductor since this post, and will be officially released at the end of this month. Details here: https://github.com/kevinblighe/PCAtools

• 0 views • link
written by Kevin Blighe Kevin Blighe

It is awesome!!! Thank you!

• 0 views • link
written by mkh • mkh

0 answers

No answers yet.

Log in to answer this question.

More posts like this