Hi, 

I am trying to make a MA plot with my RNAseq data. So, I wrote the following script:

b<-readCufflinks(system.file("extdata", package="cummeRbund"))
> genes<-b@genes
> d<-MAplot(genes,'LC','LT')

However, the following error arise:

Error in dat$x + pseudocount : non-numeric argument to binary operator

Does anyone knows how to solve this problem.

Thanks in advance.

Pablo

The error arises because you aren't using your own data, you are using the example data that comes with cummeRbund. If you want to use your own data, you have to tell readCufflinks() where those data are.