Hello list I'm having some problems with the baySeq package for RNA-Seq analysis. Basically if I filter my data (non-normalised counts at the gene level) to get rid of genes with low expression (<1 count / million as per edgeR manual - this is just a filter choice) I cannot get the plotMA.CD command in baySeq to work. It fails with: Error in plot.window(...) : need finite 'ylim' values In addition: Warning message: In max(abs((log2(Adata) - log2(Bdata))[-union(Azeros, Bzeros)])) : no non-missing arguments to max; returning -Inf