Question: deseq 2x2: highly variable genes are not among differentially expressed gene
gravatar for rahel14350
19 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 19 months ago by Michael Love11k • written 19 months ago by rahel1435010
gravatar for Michael Love
19 months ago by
Michael Love11k
United States
Michael Love11k 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 19 months ago by Michael Love11k
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: 290 users visited in the last hour