Plotting the most expressed genes after conducting DESeq2
1
0
Entering edit mode
amoaristotle ▴ 10
@user-24704
Last seen 3.4 years ago

So I had succeeded in conducting my DESeq2 without any problem> However, I want to go further to show 30 most expressed genes so I used the following code:

Note: The count.data.set.object is the DESeqDataset

vsd <- vst(count.data.set.object)

vsd <- varianceStabilizingTransformation(count.data.set.object)
norm.data = assay(vsd)

#for the heat map for the most expressed genes
library("pheatmap")

select <- order(rowMeans(counts(count.data.set.object,normalized=TRUE)),
                decreasing=TRUE)[1:100]

df <- as.data.frame(colData(count.data.set.object)[,c("condition")])

pheatmap(assay(vsd)[select,], cluster_rows=FALSE, show_rownames=FALSE,
         cluster_cols=FALSE, annotation_col=df)

After running the last code it the said:

Error in check.length("fill") : 
      'gpar' element 'fill' must not be length 0. 

How can I resolve this?

DESeq2 pheatmap apeglm • 1.1k views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 5 hours ago
United States

This is a pheatmap issue, not a DESeq2 issue. I'm not sure exactly but check the inputs to the function. Maybe df is not proper input here.

ADD COMMENT

Login before adding your answer.

Traffic: 833 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6