Please can anyone help figure out where I am wrong here. I was trying to plot heatmap from RNA-seq data that I analysed using DESeq2. I want to filter out TopVarGenes from the list of DEGs obtained but whenever I type the following code:
library("genefilter") topVarGenes <- head(order(rowVars(assay(rld)), decreasing = TRUE), 20) mat <- assay(rld)[ topVarGenes, ] mat <- mat - rowMeans(mat) anno <- as.data.frame(coldata(rld)[, c("condition")]) pheatmap(mat, annotation_col = anno)
It returns with the following errors
Error in check.length("fill") : 'gpar' element 'fill' must not be length 0
Where have I missed it then?.