Question: ntop parameter in plotPCA(). DESeq2
0
gravatar for Mthabisi Moyo
11 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 • 391 views
ADD COMMENTlink modified 11 months ago by Michael Love23k • written 11 months ago by Mthabisi Moyo10
Answer: ntop parameter in plotPCA(). DESeq2
1
gravatar for Michael Love
11 months ago by
Michael Love23k
United States
Michael Love23k 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 11 months ago by Michael Love23k
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: 301 users visited in the last hour