Skip to content
Snippets Groups Projects
Verified Commit fba1fce0 authored by Laurent Modolo's avatar Laurent Modolo
Browse files

add cpu args to kmerclust

parent 1a3ed11b
Branches
No related tags found
No related merge requests found
library(kmerclust)
library(ggplot2)
args <- commandArgs(trailingOnly = TRUE)
data <- readr::read_tsv(args[2], show_col_types = FALSE)
format(object.size(data), units = "Mb")
annotation <- parse_annotation(args[1])
count <- annotate_counts(annotation, data, args[3])
save(count, file = paste0(args[3], ".Rdata"))
res <- compare_models(count %>% dplyr::ungroup(), nboot = 100, bootsize = 0.1, core = as.numeric(args[4]))
save(res, file = paste0(args[3], "_boot.Rdata"))
res %>%
ggplot(aes(x = name, y = BIC)) +
geom_violin()
ggsave(paste0(args[3], "_BIC.pdf"))
res %>%
ggplot(aes(x = name, y = WSS_f / BSS)) +
geom_violin()
ggsave(paste0(args[3], "_WSS_BSS.pdf"))
kmerclust @ bd49510e
Subproject commit fa432617ab1e7b8eea3c4d671ce9e8948653e6f6
Subproject commit bd49510e4e4931796e2ed7811333479e74f52269
......@@ -16,7 +16,7 @@ process KMERCLUST {
script:
def args = task.ext.args ?: ''
"""
Rscript ../bin/kmerclust.R ${params.csv} ${csv} ${specie}
Rscript ../bin/kmerclust.R ${params.csv} ${csv} ${specie} ${task.cpus}
cat <<-END_VERSIONS > versions.yml
"${task.process}":
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment