This is a test version of Biostars. For the public version, visit https://www.biostars.org.
What are ntop on plotPCA vsd in DESeq2?

I am trying to understand the following DESeq PCA code:

plotPCA(vsd, intgroup = "condition", ntop = 500)

VERY different to the differentially expressed subset of genes from here:

de <- rownames(res[res$padj<0.05 & !is.na(res$padj), ]) 

plotPCA(vsd[de,], intgroup = "condition", ntop = 500)

What are thos ntop=500 selected at the beginning then? based on what are they TOP?

deseq2 rna-seq r pca plotpca

1 answer

Please read the help ?DESeq2::plotPCA

ntop: number of top genes to use for principal components, selected by highest row variance

Log in to answer this question.