Closed RHReynolds closed 2 years ago
Hey! Thanks for your note on this, I believe your are correct and I'll just confirm before pushing a fix on the dev branch of bioconductor (master branch on github)
@bschilder do you agree this is just an oversight? There is no reason not to pass genelistSpecies to prepare_genesize_control_network from bootstrap_enrichment_test? - I know you worked on this last so just want to confirm
Hi @RHReynolds I believe you are indeed correct on this, thanks for reporting it!
@Al-Murphy that would be great if you could go ahead and push the fix.
Updated in v 1.5.3, @RHReynolds you can get this now from the github master branch or alternatively, the bioconductor devel version in the next few days (usually 3 days)
1. Bug description
When using
bootstrap_enrichment_test()
with the argumentgeneSizeControl = TRUE
a warning is given thatgenelistSpecies
has not been provided and the default human will be used. I believe this is because thegenelistSpecies
argument has not been specified/inherited inprepare_genesize_control_network()
(lines 152-7 in https://github.com/NathanSkene/EWCE/blob/master/R/bootstrap_enrichment_test.R).Console output
Expected behaviour
I would have expected that the control network uses the same species as specified for my hit list.
2. Reproducible example
Code
I used the example provided in the manual page for
bootstrap_enrichment_test()
, but added in the argument to control for GC/length.3. Session info
(Add output of the R function
utils::sessionInfo()
below. This helps us assess version/OS conflicts which could be causing bugs.)