In the following article. There are two ways of getting highly variable genes:
resSig <- subset(res, padj < 0.1) head(resSig[ order(resSig$log2FoldChange, decreasing=TRUE), ])
Getting genes that are statistically significant, then ordering them by the
log2FoldChange which is
The column log2FoldChange is the effect size estimate. It tells us how much the gene’s expression seems to have changed due to treatment with dexamethasone in comparison to untreated samples.
And the second one:
topVarGenes <- head(order(rowVars(assay(rld)),decreasing=TRUE),20)
I am struggling to understand what are the differences between the two ways of getting the most variable genes. What is the difference?