Question

What to do when a bug is found in a BioConductor package (Glimma)

0

Entering edit mode

Sam ▴ 10

@sam-21502

Last seen 3 days ago

Jerusalem

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

The following code

library(Glimma)

DE <- as.data.frame(cbind(202197.7, -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?

glimma software error • 993 views

ADD COMMENT • link updated 4.7 years ago by Gordon Smyth 50k • written 4.7 years ago by Sam ▴ 10

1

Entering edit mode

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 REPLY • link 4.7 years ago Shian Su ▴ 40

score 1 · Answer 1 · 2019-08-17

1

Entering edit mode

Gordon Smyth 50k

@gordon-smyth

Last seen 36 minutes ago

WEHI, Melbourne, Australia

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

https://github.com/Shians/Glimma/issues

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

ADD COMMENT • link 4.7 years ago Gordon Smyth 50k