Question: deseq 2x2: highly variable genes are not among differentially expressed gene
gravatar for rahel14350
21 months ago by
United States
rahel1435010 wrote:


I did deseq 2x2 model. When I do a pheatmap for highly variable genes (using rld) the outcome are the genes that they are not even among the differentially expressed genes.

Here is my command for pheatmap:

topVarGenes <- head(order(rowVars(assay(rld)),decreasing=TRUE),20)
mat <- assay(rld)[ topVarGenes, ]
mat <- mat - rowMeans(mat)
df <-[,c("group","condition")])
pheatmap(mat, annotation_col=df)

For example for the gene with highest rld over samples, the basemean in dds is 98000 but the pvalue and padj is NA.

Can anyone help me with this issue? 

Many thanks in advance,


ADD COMMENTlink modified 21 months ago by Michael Love12k • written 21 months ago by rahel1435010
gravatar for Michael Love
21 months ago by
Michael Love12k
United States
Michael Love12k wrote:

topVarGenes looks at the row variance of the transformed values regardless of which samples come from which group. Differential expression testing asks whether the difference across group is large relative to the within-group variance. So these are different ways of ranking genes.

ADD COMMENTlink written 21 months ago by Michael Love12k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 111 users visited in the last hour