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