Mercurial > repos > galaxyp > custom_pro_db_annotation_data_manager

annotate data_manager/customProDB_annotation.R @ 6:c922df2f9c2a draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
author galaxyp
date Tue, 31 May 2016 13:02:29 -0400
parents a0ca66b32792
children 7399b23db662
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
1 #!/usr/bin/env Rscript
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
2
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
3 initial.options <- commandArgs(trailingOnly = FALSE)
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
4 script_parent_dir <- dirname(sub("--file=", "", initial.options[grep("--file=", initial.options)]))
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
5
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
6 ## begin warning handler
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
7 withCallingHandlers({
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
8
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
9 library(methods) # Because Rscript does not always do this
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
10
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
11 options('useFancyQuotes' = FALSE)
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
12
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
13 suppressPackageStartupMessages(library("optparse"))
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
14 suppressPackageStartupMessages(library("RGalaxy"))
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
15
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
16
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
17 option_list <- list()
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
18 option_list$dbkey <- make_option('--dbkey', type='character')
5
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
19 option_list$dbsnp <- make_option('--dbsnp', type='character')
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
20 option_list$cosmic <- make_option('--cosmic', type='logical')
0
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
21 option_list$outputFile <- make_option('--outputFile', type='character')
3
af0a098e15bd planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 2
diff changeset
22 option_list$dbkey_description <- make_option('--dbkey_description', type='character')
0
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
23
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
24 opt <- parse_args(OptionParser(option_list=option_list))
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
25
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
26
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
27 customProDB_annotation <- function(
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
28 dbkey = GalaxyCharacterParam(required=TRUE),
5
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
29 dbsnp_str = GalaxyCharacterParam(required=FALSE),
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
30 cosmic = GalaxyLogicalParam(required=FALSE),
3
af0a098e15bd planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 2
diff changeset
31 dbkey_description = GalaxyCharacterParam(required=FALSE),
0
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
32 outputFile = GalaxyOutput("output","json"))
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
33 {
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
34 if (!file.exists(outputFile))
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
35 {
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
36 gstop("json params file does not exist")
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
37 }
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
38
3
af0a098e15bd planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 2
diff changeset
39 if (length(dbkey_description) < 1)
af0a098e15bd planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 2
diff changeset
40 {
af0a098e15bd planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 2
diff changeset
41 dbkey_description = dbkey
af0a098e15bd planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 2
diff changeset
42 }
af0a098e15bd planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 2
diff changeset
43
5
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
44 if (length(dbsnp_str) > 0)
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
45 {
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
46 dbsnp = dbsnp_str
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
47 }
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
48 else
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
49 {
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
50 dbsnp = NULL
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
51 }
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
52
6
c922df2f9c2a planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 5
diff changeset
53 if (length(cosmic) > 0 && cosmic && grepl("^hg", dbkey))
5
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
54 {
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
55 use_cosmic = TRUE
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
56 }
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
57 else
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
58 {
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
59 use_cosmic = NULL
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
60 }
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
61
0
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
62 suppressPackageStartupMessages(library(rjson))
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
63 params = fromJSON(file=outputFile)
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
64 target_directory = params$output_data[[1]]$extra_files_path
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
65 dir.create(target_directory)
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
66
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
67 tryCatch(
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
68 {
3
af0a098e15bd planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 2
diff changeset
69 file.remove(outputFile)
0
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
70 }, error=function(err)
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
71 {
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
72 gstop("failed to remove json params file after reading")
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
73 })
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
74
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
75 ucscTableCodingFastaURL = paste("http://genome.ucsc.edu/cgi-bin/hgTables?db=", dbkey, "&hgSeq.cdsExon=on&hgSeq.granularity=gene&hgSeq.casing=exon&hgSeq.repMasking=lower&hgta_doGenomicDna=get+sequence&hgta_group=genes&hgta_track=refGene&hgta_table=refGene&hgta_regionType=genome", sep="")
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
76 ucscTableProteinFastaURL = paste("http://genome.ucsc.edu/cgi-bin/hgTables?db=", dbkey, "&hgta_geneSeqType=protein&hgta_doGenePredSequence=submit&hgta_track=refGene&hgta_table=refGene", sep="")
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
77 codingFastaFilepath = paste(target_directory, "/", dbkey, ".cds.fa", sep="")
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
78 proteinFastaFilepath = paste(target_directory, "/", dbkey, ".protein.fa", sep="")
3
af0a098e15bd planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 2
diff changeset
79
0
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
80 suppressPackageStartupMessages(library(customProDB))
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
81 options(timeout=3600)
3
af0a098e15bd planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 2
diff changeset
82
0
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
83 cat(paste("Downloading coding FASTA from:", ucscTableCodingFastaURL, "\n"))
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
84 download.file(ucscTableCodingFastaURL, codingFastaFilepath, quiet=T, mode='wb')
3
af0a098e15bd planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 2
diff changeset
85
0
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
86 cat(paste("Downloading protein FASTA from:", ucscTableProteinFastaURL, "\n"))
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
87 download.file(ucscTableProteinFastaURL, proteinFastaFilepath, quiet=T, mode='wb')
3
af0a098e15bd planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 2
diff changeset
88
af0a098e15bd planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 2
diff changeset
89 cat(paste("Preparing Refseq annotation files\n"))
5
a0ca66b32792 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 98f3323d4b92a4a3668788b3e1bad1a6529fd6df-dirty
galaxyp
parents: 4
diff changeset
90 customProDB::PrepareAnnotationRefseq(genome=dbkey, CDSfasta=codingFastaFilepath, pepfasta=proteinFastaFilepath, annotation_path=target_directory, dbsnp=dbsnp, cosmic=use_cosmic)
0
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
91
3
af0a098e15bd planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 2
diff changeset
92 outputPath = paste(dbkey, "/customProDB", sep="")
0
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
93 output = list(data_tables = list())
4
b42974daa659 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents: 3
diff changeset
94 output[["data_tables"]][["customProDB"]]=c(path=outputPath, name=dbkey_description, dbkey=dbkey, value=dbkey)
0
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
95 write(toJSON(output), file=outputFile)
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
96 }
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
97
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
98
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
99 params <- list()
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
100 for(param in names(opt))
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
101 {
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
102 if (!param == "help")
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
103 params[param] <- opt[param]
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
104 }
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
105
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
106 setClass("GalaxyRemoteError", contains="character")
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
107 wrappedFunction <- function(f)
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
108 {
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
109 tryCatch(do.call(f, params),
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
110 error=function(e) new("GalaxyRemoteError", conditionMessage(e)))
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
111 }
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
112
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
113
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
114 suppressPackageStartupMessages(library(RGalaxy))
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
115 do.call(customProDB_annotation, params)
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
116
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
117 ## end warning handler
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
118 }, warning = function(w) {
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
119 cat(paste("Warning:", conditionMessage(w), "\n"))
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
120 invokeRestart("muffleWarning")
663ee21a8609 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tools/bumbershoot/data_manager_customProDB commit 141369f97aa2804d2bbfd9ed620ea2a5574994c2-dirty
galaxyp
parents:
diff changeset
121 })