Question: ntop parameter in plotPCA(). DESeq2
0
gravatar for Mthabisi Moyo
14 months ago by
Mthabisi Moyo10 wrote:

Given that the PCA plot is likely to change somewhat depending on the number of genes you decide to specify with the ntop parameter, are there any recommendations on how to best set this value besides arbitrarily setting it at the default of 500/1000? Could including all genes have a negative effect if a lot of the genes have low variance?

Thanks!

 

deseq2 • 494 views
ADD COMMENTlink modified 14 months ago by Michael Love24k • written 14 months ago by Mthabisi Moyo10
Answer: ntop parameter in plotPCA(). DESeq2
1
gravatar for Michael Love
14 months ago by
Michael Love24k
United States
Michael Love24k wrote:

"Could including all genes have a negative effect if a lot of the genes have low variance?"

Not really. The low count noise is dealt with by vst() or rlog(). It's up to you what ntop to choose, choosing a small-ish number, e.g. 500, focuses the PCA plot on the most variable genes across samples, which is often the DE genes across condition (although this information is not used to pick the genes in the PCA plot).

ADD COMMENTlink written 14 months ago by Michael Love24k
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 16.09
Traffic: 141 users visited in the last hour