# HG changeset patch # User artbio # Date 1615158788 0 # Node ID 67213c4eefa6b041df2f95ff244c5f83290944c1 # Parent e66b91475343bb9e861337e6dcad2cc0fd5f91d3 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit 10ad3a0ca7cd23ad1e0940844147e1d1b3d069f0" diff -r e66b91475343 -r 67213c4eefa6 sanitize_snv.py --- a/sanitize_snv.py Sun Mar 07 14:21:06 2021 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,24 +0,0 @@ -import sys - -handle = open(sys.argv[1], 'r') -out = open(sys.argv[2], 'w') -for line in handle: - if line[0] == '#': - out.write(line) - continue - linelist = line.split('\t') - refcol = linelist[0].split('chr') - infocol = linelist[7].split('INDEL') - if len(infocol) > 1: - continue - if len(refcol) > 1: - refcol = refcol[1] - else: - refcol = refcol[0] - if refcol not in ['1', '2', '3', '4', '5', '6', '7', '8', '9', '10', - '11', '12', '13', '14', '15', '16', '17', '18', '19', - '20', '21', '22', 'X', 'Y']: - continue - else: - linelist[0] = refcol - out.write('\t'.join(linelist)) diff -r e66b91475343 -r 67213c4eefa6 segmentation_sequenza.R --- a/segmentation_sequenza.R Sun Mar 07 14:21:06 2021 +0000 +++ b/segmentation_sequenza.R Sun Mar 07 23:13:08 2021 +0000 @@ -1,10 +1,13 @@ -options(warn = -1) +options(show.error.messages = F, error = function() { + cat(geterrmessage(), file = stderr()); q("no", 1, F) }) # load packages that are provided in the conda env library(optparse) library(sequenza) library(BiocParallel) library(tidyverse) +library(readr) + option_list <- list( make_option( c("-i", "--input"), @@ -41,8 +44,14 @@ "chr10", "chr11", "chr12", "chr13", "chr14", "chr15", "chr16", "chr17", "chr18", "chr19", "chr20", "chr21", "chr22") +Sys.setenv(TZDIR = "/usr/share/zoneinfo/") +Sys.setenv(TZ = "US/Eastern") +options(tz = "US/Eastern") +Sys.timezone(location = TRUE) + + segfile <- sequenza.extract(data_file, - verbose = TRUE, + verbose = FALSE, chromosome.list = chrom_list) ## Estimation of cellularity and ploidy diff -r e66b91475343 -r 67213c4eefa6 sequenza_index.xml --- a/sequenza_index.xml Sun Mar 07 14:21:06 2021 +0000 +++ b/sequenza_index.xml Sun Mar 07 23:13:08 2021 +0000 @@ -1,4 +1,4 @@ - + diff -r e66b91475343 -r 67213c4eefa6 sequenza_to_hrdtools_input.R --- a/sequenza_to_hrdtools_input.R Sun Mar 07 14:21:06 2021 +0000 +++ b/sequenza_to_hrdtools_input.R Sun Mar 07 23:13:08 2021 +0000 @@ -1,4 +1,5 @@ -options(warn = -1) +options(show.error.messages = F, error = function() { + cat(geterrmessage(), file = stderr()); q("no", 1, F) }) # load packages that are provided in the conda env @@ -61,6 +62,6 @@ print(reformatted) reformatted %>% - write.table(opt$output, quote=F, row.names=F, sep="\t") + write.table(opt$output, quote = F, row.names = F, sep = "\t") message(sprintf("Output written to %s", opt$output)) diff -r e66b91475343 -r 67213c4eefa6 snvtocnv.xml --- a/snvtocnv.xml Sun Mar 07 14:21:06 2021 +0000 +++ b/snvtocnv.xml Sun Mar 07 23:13:08 2021 +0000 @@ -1,16 +1,16 @@ - + macro.xml - samtools sequenza-utils r-sequenza r-optparse bioconductor-biocparallel r-tidyverse + tzdata