Mercurial > repos > proteore > proteore_prot_features
annotate protein_features.R @ 1:bfc679370c64 draft
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
author | proteore |
---|---|
date | Fri, 16 Feb 2018 04:06:16 -0500 |
parents | |
children | 867d47ff782c |
rev | line source |
---|---|
1
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
1 # Read file and return file content as data.frame |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
2 readfile = function(filename, header) { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
3 if (header == "true") { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
4 # Read only first line of the file as header: |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
5 headers <- read.table(filename, nrows = 1, header = FALSE, sep = "\t", stringsAsFactors = FALSE, fill = TRUE, na.strings=c("", "NA"), blank.lines.skip = TRUE) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
6 #Read the data of the files (skipping the first row) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
7 file <- read.table(filename, skip = 1, header = FALSE, sep = "\t", stringsAsFactors = FALSE, fill = TRUE, na.strings=c("", "NA"), blank.lines.skip = TRUE) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
8 # Remove empty rows |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
9 file <- file[!apply(is.na(file) | file == "", 1, all), , drop=FALSE] |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
10 #And assign the header to the data |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
11 names(file) <- headers |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
12 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
13 else { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
14 file <- read.table(filename, header = FALSE, sep = "\t", stringsAsFactors = FALSE, fill = TRUE, na.strings=c("", "NA"), blank.lines.skip = TRUE) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
15 # Remove empty rows |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
16 file <- file[!apply(is.na(file) | file == "", 1, all), , drop=FALSE] |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
17 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
18 return(file) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
19 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
20 |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
21 protein_features = function() { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
22 args <- commandArgs(TRUE) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
23 if(length(args)<1) { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
24 args <- c("--help") |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
25 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
26 |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
27 # Help section |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
28 if("--help" %in% args) { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
29 cat("Selection and Annotation HPA |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
30 Arguments: |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
31 --inputtype: type of input (list of id or filename) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
32 --input: input |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
33 --nextprot: path to nextprot information file |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
34 --column: the column number which you would like to apply... |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
35 --header: true/false if your file contains a header |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
36 --type: the type of input IDs (UniProt/EntrezID) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
37 --argsP1: IsoPoint,SeqLength,MW |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
38 --argsP2: Chr,SubcellLocations |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
39 --argsP3: Diseases |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
40 --output: text output filename \n") |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
41 q(save="no") |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
42 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
43 |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
44 # Parse arguments |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
45 parseArgs <- function(x) strsplit(sub("^--", "", x), "=") |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
46 argsDF <- as.data.frame(do.call("rbind", parseArgs(args))) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
47 args <- as.list(as.character(argsDF$V2)) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
48 names(args) <- argsDF$V1 |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
49 |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
50 inputtype = args$inputtype |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
51 if (inputtype == "copypaste") { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
52 input = strsplit(args$input, " ")[[1]] |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
53 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
54 else if (inputtype == "tabfile") { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
55 filename = args$input |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
56 ncol = args$column |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
57 # Check ncol |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
58 if (! as.numeric(gsub("c", "", ncol)) %% 1 == 0) { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
59 stop("Please enter an integer for level") |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
60 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
61 else { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
62 ncol = as.numeric(gsub("c", "", ncol)) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
63 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
64 header = args$header |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
65 # Get file content |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
66 file = readfile(filename, header) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
67 # Extract Protein IDs list |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
68 input = c() |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
69 for (row in as.character(file[,ncol])) { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
70 input = c(input, strsplit(row, ";")[[1]][1]) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
71 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
72 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
73 nextprot_file = args$nextprot |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
74 nextprot = human_id_map = read.table(nextprot_file, header = TRUE, sep = "\t", stringsAsFactors = FALSE, fill = TRUE, na.strings = "") |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
75 typeid = args$type |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
76 P1_args = strsplit(args$argsP1, ",")[[1]] |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
77 P2_args = strsplit(args$argsP2, ",")[[1]] |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
78 P3_args = strsplit(args$argsP3, ",")[[1]] |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
79 output = args$output |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
80 |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
81 # Change the sample ids if they are uniprot ids to be able to match them with |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
82 # Nextprot data |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
83 if (typeid=="uniprot"){ |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
84 input = gsub("^","NX_",input) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
85 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
86 |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
87 # Select user input protein ids in nextprot |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
88 if ((length(input[input %in% nextprot[,1]]))==0){ |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
89 write.table("None of the input ids are can be found in Nextprot",file=filename,sep="\t",quote=FALSE,col.names=TRUE,row.names=FALSE) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
90 } else { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
91 names = c() |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
92 res = matrix(nrow=length(input), ncol=0) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
93 |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
94 # Get information from neXtProt |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
95 if (length(P1_args)>0) { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
96 for (arg in P1_args) { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
97 names = c(names, arg) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
98 info = nextprot[match(input, nextprot["NextprotID"][,]),][arg][,] |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
99 res = cbind(res, info) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
100 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
101 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
102 if (length(P2_args)>0) { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
103 for (arg in P2_args) { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
104 names = c(names, arg) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
105 info = nextprot[match(input, nextprot["NextprotID"][,]),][arg][,] |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
106 res = cbind(res, info) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
107 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
108 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
109 if (length(P3_args)>0) { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
110 for (arg in P3_args) { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
111 names = c(names, arg) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
112 info = nextprot[match(input, nextprot["NextprotID"][,]),][arg][,] |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
113 res = cbind(res, info) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
114 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
115 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
116 |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
117 # Write output |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
118 if (inputtype == "copypaste") { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
119 res = cbind(as.matrix(input), res) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
120 names = c(typeid, names) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
121 colnames(res) = names |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
122 write.table(res, output, row.names = FALSE, sep = "\t", quote = FALSE) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
123 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
124 else if (inputtype == "tabfile") { |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
125 names = c(names(file), names) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
126 output_content = cbind(file, res) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
127 colnames(output_content) = names |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
128 write.table(output_content, output, row.names = FALSE, sep = "\t", quote = FALSE) |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
129 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
130 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
131 |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
132 } |
bfc679370c64
planemo upload commit 5774fd6a5a746f36f6bf4671a51a39ea2b978300-dirty
proteore
parents:
diff
changeset
|
133 protein_features() |