What are ntop on plotPCA vsd in DESeq2?
1
1
Entering edit mode
3.4 years ago
ecg1g15 ▴ 30

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 • 2.9k views
ADD COMMENT
2
Entering edit mode
3.4 years ago
ATpoint 81k

Please read the help ?DESeq2::plotPCA

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

ADD COMMENT

Login before adding your answer.

Traffic: 1604 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6