lmrodriguezr / nonpareil

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

Stats summary return a matrix of lists #62

Open fgvieira opened 1 year ago

fgvieira commented 1 year ago

I am trying to make a snakemake wrapper to plot nonpareil curves:

[...]
curves <- Nonpareil.curve.batch(snakemake@input[["npo"]],
                label = params[["label"]],
                labels = params[["labels"]],
                col = params[["col"]],
                enforce.consistency = params[["enforce_consistency"]],
                star = params[["star"]],
                correction.factor = params[["correction_factor"]],
                weights.exp = params[["weights_exp"]],
                skip.model = params[["skip_model"]],
                plot = FALSE
                )

# Get stats
stats <- summary(curves)
[...]

stats is a numeric matrix but, when skip.model=True, stats is a matrix of numeric lists. Is this intended?

lmrodriguezr commented 2 months ago

Hi @fgvieira That doesn't sound intentional, I'll address this next. Thank you!