Create report rule fails for test data
While testing the pipeline with Singularity support (branch) I utilized this test routine and run into an issue while generating the final report in R
:
Quitting from lines 539-549 (ncov_minipipe.Rmd)
Error in seq.default(min, max, by = by) : 'from' must be a finite number
Calls: <Anonymous> ... <Anonymous> -> f -> <Anonymous> -> seq -> seq.default
The error belongs to
ggplot(dt.coverage, aes(position, depth)) +
geom_point(size = 0.1, colour = "#4444DD") +
labs(title = "coverage distribution") +
theme(plot.title = element_text(hjust = 0.5),
axis.text.x = element_text(angle = 45, hjust = 0.5),
legend.position = "none") +
scale_y_continuous(trans = "log10") +
facet_wrap(~sample, ncol = 3, scales = "free_y")
in particular
scale_y_continuous(trans = "log10") +
that I can comment and then the script runs through.
I think this is due to the small example input resulting in too few values for log transformation.
I will try to directly fix this in the Singularity branch.