search

lmrodriguezr / nonpareil

Estimate metagenomic coverage and sequence diversity
http://enve-omics.ce.gatech.edu/nonpareil/
Other
42 stars 11 forks source link

Re-add factor to curve? #52

Closed sjaenick closed 2 years ago

sjaenick commented 2 years ago

Hi,

earlier releases of the package (I checked v2.3) supported a 'factor' argument for the curve that allowed to easily switch from bp to Mbp/Gbp/.. - but also mentioned that this could possibly influence the fit of the model. Would it be possible to re-introduce something similar again, maybe not for the curve but only as an additional argument to the plot function (so model fit wouldn't change)?

lmrodriguezr commented 2 years ago

Hello,

That's a very good suggestion! You can currently do it with:

Nonpareil.curve(file, xaxt = "n")
axis(1, at = 10^c(3, 6, 9 , 12), labels = paste(1, c("Kbp", "Mbp", "Gbp", "Tbp")))

I'll add that to the documentation momentarily so it's easily accessible.

Thanks! Miguel.