Mercurial > repos > iuc > lineagespot
diff lineagespot_verbose.R @ 1:99494998688a draft default tip
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/lineagespot commit 6a6a37f2574954dae65f9ec407fe38594ed37659
author | iuc |
---|---|
date | Sun, 25 Feb 2024 09:49:20 +0000 |
parents | 6ddf5a9ce4a5 |
children |
line wrap: on
line diff
--- a/lineagespot_verbose.R Tue Aug 08 15:12:08 2023 +0000 +++ b/lineagespot_verbose.R Sun Feb 25 09:49:20 2024 +0000 @@ -13,8 +13,8 @@ # Rscript ${__tool_directory__}/lineagespot_verbose.R --in_vcf ${__tool_directory__}/test-data/extdata/vcf-files --in_gff3 ${__tool_directory__}/test-data/extdata/NC_045512.2_annot.gff3 --in_ref ${__tool_directory__}/test-data/extdata/ref --in_voc "B.1.617.2, B.1.1.7, B.1.351, P.1" --in_threshold 0.8 # Set up R error handling to go to stderr options(show.error.messages = FALSE, error = function() { - cat(geterrmessage(), file = stderr()) - q("no", 1, FALSE) + cat(geterrmessage(), file = stderr()) + q("no", 1, FALSE) }) # Import required libraries @@ -22,9 +22,9 @@ library_path <- .libPaths() suppressPackageStartupMessages({ - library("getopt", lib.loc = library_path) - library("data.table", lib.loc = library_path) - library("lineagespot", lib.loc = library_path) + library("getopt", lib.loc = library_path) + library("data.table", lib.loc = library_path) + library("lineagespot", lib.loc = library_path) }) @@ -36,24 +36,26 @@ # Get options using the spec as defined by the enclosed list # Read the options from the default: commandArgs(TRUE) option_specification <- matrix(c( - "in_vcf", "vcf", 1, "character", - "in_gff3", "gff3", 1, "character", - "in_ref", "ref", 1, "character", - "in_voc", "voc", 2, "character", - "in_threshold", "thr", 2, "double" + "in_vcf", "vcf", 1, "character", + "in_gff3", "gff3", 1, "character", + "in_ref", "ref", 1, "character", + "in_voc", "voc", 2, "character", + "in_threshold", "thr", 2, "double" ), byrow = TRUE, ncol = 4) options <- getopt(option_specification) if (!is.null(options$in_voc) && is.character(options$in_voc)) { - options$in_voc <- unlist(strsplit(options$in_voc, split = ",")) + options$in_voc <- unlist(strsplit(options$in_voc, split = ",")) } -result <- lineagespot(vcf_folder = options$in_vcf, - ref_folder = options$in_ref, - gff3_path = options$in_gff3, - voc = options$in_voc, - AF_threshold = options$in_threshold) +result <- lineagespot( + vcf_folder = options$in_vcf, + ref_folder = options$in_ref, + gff3_path = options$in_gff3, + voc = options$in_voc, + AF_threshold = options$in_threshold +) # Write output to new file which will be recognized by Galaxy