Error when trying to plot sample distance from and assay object generated with DESeq2
1
0
Entering edit mode
colaneri ▴ 30
@colaneri-7770
Last seen 5.1 years ago
United States

Hi  Im following RNAseq workflow for differential gene expression

Im using R command line. 

library(gplots)

library(RColorBrewer)

sampleDists<- dist(t(assay(rld)))

sampleDistsMatrix<- as.matrix(sampleDists)

xnames<- paste(colData(ddsFull)$SampleName,colData(ddsFull)$Sex,colData(ddsFull)$treatment)

dimnames(sampleDistsMatrix)<- list(xnames,xnames)

colors<- colorRampPalette(rev(brewer.pal(9,"Blues")))(255)

heatmap.2(sampleDistsMatrix,trace="none",col=colors,margins=c(10,10))

But i got this error: Error in plot.new() : internal read error in PDF_endpage Error in par(op) : parameter "i" in "mfg" is out of range what is what Im doing wrong?
deseq2 heatmap.2 • 1.5k views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 11 hours ago
United States

That seems like some error coming from heatmap.2, so I'm not sure about debugging. 

We are moving toward recommending the pheatmap package instead. Can you try this code from the workflow (search the page for "Sample distances")

http://www.bioconductor.org/help/workflows/rnaseqGene/

ADD COMMENT

Login before adding your answer.

Traffic: 436 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