CummeRbund: I can't create a dispersion plot - help please?
0
1
Entering edit mode
6.4 years ago
BioinfGuru ★ 1.6k

Hi all,

I am now somewhat familiar with cummeRbund and have all my plots already except for one: a dispersion plot. I know dispersionPlot only runs on a cuffData object, however I don;t know how to create a cuffData object. Ive looked at the cuffData help file but the usage is not clear to me.

My code

library("cummeRbund")

library("plyr")

library("Hmisc")

setwd("[pathto]/r_working_directory")

disp<-dispersionPlot(genes(cuff))

disp

So the code is successful (creates a "disp" value in the environment window)

however when i try to veiw it I get the following error:

Error in $<-.data.frame(*tmp*, "SCALE_X", value = 1L) : replacement has 1 row, data has 0 In addition: Warning message: In max(panels$ROW) : no non-missing arguments to max; returning -Inf

I really have no idea how to fix this. Do need to create a cuffData object, if so how? Or is there a simpler fix?

RNA-Seq cummeRbund dispersionPlot cuffData class • 1.9k views
1
Entering edit mode

The gtf file is merged.gtf file from cuffmerge, right? If not, can you try with merge.gtf file?

0
Entering edit mode

Right you are,.... I had the wrong gtf.......still didnt solve the problem though

0
Entering edit mode

Can you check if cuff is created properly?

cuff

If it is made properly, you should see number of samples, genes, isoforms, TSS, CDSs etc.

0
Entering edit mode

I've actually found out that it has to do with the versions I'm using. Thank you.