Question: What to do when a bug is found in a BioConductor package (Glimma)
gravatar for Samuel
4 weeks ago by
Samuel0 wrote:

As far as I can tell, I have found a bug in Glimma.

The following code


DE <-, -1.410116), byrow=T)

colnames(DE) <- c("baseMean", "log2FoldChange")
DE$sig <- 1

glMDPlot(x=DE, xval="baseMean", yval="log2FoldChange", counts=NULL,
         anno=NULL,samples=NULL, group = NULL, status=DE$sig, 
         side.main = "baseMean",html="One_Gene_from_vector_final",
         display.columns = c("baseMean","log2FoldChange"))

produces the following graph :

glimma plot

As you can see, the baseMean is rounded (incorrectly) to 202200.

What is the etiquette in this case? Should I email the developers directly?

software error glimma • 93 views
ADD COMMENTlink modified 4 weeks ago by Gordon Smyth38k • written 4 weeks ago by Samuel0

Thanks for the bug report Samuel. As Gordon mentioned, I would prefer bugs be submitted through Github issues. The reason for the bug is that I've rounded all values to 4 significant places to avoid values with tens of digits, it's clear that this causes confusion for large numeric values so I'll update the rounding strategy soon.

ADD REPLYlink written 29 days ago by Shian Su20
Answer: What to do when a bug is found in a BioConductor package (Glimma)
gravatar for Gordon Smyth
4 weeks ago by
Gordon Smyth38k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth38k wrote:

If you type packageDescription("Glimma"), or visit the Glimma landing page at the, you will see that Glimma has an issues site for reporting bugs or other issues:

I'm not sure that rounding would be considered a bug though.

ADD COMMENTlink modified 29 days ago • written 4 weeks ago by Gordon Smyth38k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 233 users visited in the last hour