annotate scripts/contamination_plot.R @ 1:253c9448f524 draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
author iuc
date Mon, 03 Jun 2019 14:55:24 -0400
parents 582b7bd4ae4c
children 5ade5cf200da
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
1 #!/usr/bin/env R
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
2
1
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
3 suppressPackageStartupMessages(require(ggplot2))
0
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
4
1
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
5 log10histoPlot <- function(title="", columncounts){
0
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
6 #' Log10 histogram plot
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
7 #'
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
8 #' @param columncounts colSums(input_matrix)
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
9 #' @param title Title of plot
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
10 #' @return ggplot grob
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
11 dfer <- data.frame(colcounts=log10(columncounts))
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
12 p1 <- ggplot(dfer, aes(x=colcounts)) +
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
13 geom_histogram(binwidth = 0.05, color="black",fill="white") +
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
14 theme(plot.title = element_text(hjust = 0.5)) +
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
15 labs(title=title, y="Frequency", x="Library Size (Log10)")
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
16
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
17 return(p1)
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
18 }
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
19
1
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
20 ## This is calculated by the first call to contaminationPlot
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
21 ## and then re-used by the second call.
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
22 ylim.max = NULL
0
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
23
1
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
24 contaminationPlot <- function(title, columndata, barcode.data, plate.data, RAW)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
25 {
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
26 coldata = data.frame(colsums=columndata)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
27 maxval = max(coldata)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
28 # Set once and once only
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
29 if (is.null(ylim.max)){
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
30 ylim.max <<- maxval + 500
0
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
31 }
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
32
1
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
33 drawPlates <- function(plate.data){ # 0 384 768
0
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
34
1
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
35 plate.boundaries <- plate.data["filtered.plates",]
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
36 if (RAW){
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
37 plate.boundaries <- plate.data["unfilter.plates",]
0
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
38 }
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
39
1
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
40 plate.minyval <- -200
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
41 plate.color <- "grey"
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
42 plate.text.color <- "red"
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
43 plate.text.alpha <- 0.5
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
44 plate.height <- 2* maxval / 5
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
45 plate.spacing <- 12
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
46 plate.height.text <- plate.height - 3000
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
47 plate.text.size <- 3
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
48
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
49
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
50 zzz <- lapply(names(plate.boundaries), function(plate.name){
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
51 plate.num = as.integer(sub("P","",plate.name))
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
52 plate.xval = plate.boundaries[[plate.name]]
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
53
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
54 ## If not first plate -- inner boundary, print next plate name to the right
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
55 if (plate.name != names(plate.boundaries)[length(plate.boundaries)]){
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
56 gplot <<- gplot + annotate("text", x=plate.xval + plate.spacing, label=paste("Plate", plate.num + 1),
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
57 size=plate.text.size, y=plate.height*0.7, angle=-90,
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
58 color=plate.text.color, alpha=plate.text.alpha)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
59 }
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
60 ## If not last plate -- inner boundary, print previous name to the left
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
61 if (plate.name != names(plate.boundaries)[1]){
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
62 gplot <<- gplot + annotate("text", x=plate.xval - plate.spacing, label=paste("Plate", plate.num),
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
63 size=plate.text.size, y=plate.height*0.7, angle=90,
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
64 color=plate.text.color, alpha=plate.text.alpha)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
65 }
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
66 gplot <<- gplot + geom_segment(aes(x=plate.xval, xend=plate.xval, y=plate.minyval, yend=plate.height), color=plate.color, lty=1, size=1)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
67 })
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
68
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
69 ## Connecting bar underneath
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
70 gplot <<- gplot + geom_segment(aes(x=min(unname(unlist(plate.boundaries))), xend=max(unname(unlist(plate.boundaries))),
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
71 y=plate.minyval, yend=plate.minyval), color=plate.color, lty=1, size=1)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
72 }
0
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
73
582b7bd4ae4c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit 6f73edc667e61fabdab8b24a7ff40942588fee5b
iuc
parents:
diff changeset
74
1
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
75 drawBatches <- function(barcode.data, plate.data){
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
76 batch.minyval = -200
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
77 divide.color <- "grey"
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
78 divide.style <- 2 ## dashed
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
79 divide.size <- 0.3
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
80 boundary.color <- "blue"
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
81 boundary.style <- 1
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
82 boundary.size <- 0.5
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
83 divide.text.true.color <- "blue"
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
84 divide.text.false.color <- "black"
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
85 divide.text.name.color <- "grey"
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
86 batch.text.height <- maxval * 4/5
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
87 batch.text.tf.height <- maxval * 2/5
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
88 batch.text.tf.spacing <- 12
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
89 batch.text.alpha <- 0.8
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
90 batch.text.size <- 3
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
91 batch.label.text.size <- 6
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
92 batch.height <- maxval * 4/5
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
93 divide.height <- maxval * 4/5
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
94 batch.spacing <- 24
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
95
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
96 boundary.pos <- barcode.data["filtered_positions",]
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
97 plate.pos <- unname(unlist(plate.data["filtered.plates",])) ## just want values
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
98
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
99 if (RAW){
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
100 boundary.pos <- barcode.data["unfilter_positions",]
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
101 divider.pos <- barcode.data["filter_in_unfilter",]
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
102 plate.pos <- unname(unlist(plate.data["unfilter.plates",])) ## just want values
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
103 }
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
104
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
105
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
106 zzz <- lapply(names(boundary.pos), function(batch.name){
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
107 batch.num = as.integer(sub("B","",batch.name))
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
108 batch.xval = unname(unlist(boundary.pos[[batch.name]]))
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
109
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
110 # Plot boundary line only if not intersecting with a plate line
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
111 if (!(batch.xval %in% plate.pos)){
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
112 gplot <<- gplot + geom_segment(aes(x=batch.xval, xend=batch.xval, y=batch.minyval, yend=batch.height),
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
113 color=boundary.color, lty=boundary.style, size=boundary.size)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
114 }
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
115
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
116 ## Plot labels (except P0)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
117 if (batch.name != "B0"){
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
118 if (RAW){
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
119 divide.xval = unname(unlist(divider.pos[[batch.name]]))
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
120
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
121 ## Plot the divider line
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
122 gplot <<- gplot + geom_segment(aes(x=divide.xval, xend=divide.xval, y=batch.minyval, yend=divide.height),
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
123 color=divide.color, lty=divide.style, size=divide.size)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
124
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
125 ## Plot the True/False labels
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
126 gplot <<- gplot + annotate("text", x=divide.xval - batch.text.tf.spacing, label="True positives",
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
127 size=batch.text.size, y=batch.text.tf.height, angle=+90,
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
128 color=divide.text.true.color, alpha=batch.text.alpha)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
129 gplot <<- gplot + annotate("text", x=divide.xval + batch.text.tf.spacing, label="False positives",
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
130 size=batch.text.size, y=batch.text.tf.height, angle=-90,
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
131 color=divide.text.false.color, alpha=batch.text.alpha)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
132 ## Plot the Batch names
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
133 gplot <<- gplot + annotate("text", x=divide.xval, label=batch.name,
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
134 size=batch.label.text.size, y=batch.text.height, angle=0,
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
135 color = divide.text.name.color, alpha=batch.text.alpha)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
136 } else {
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
137 ## Plot the Batch names between blue lines
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
138 gplot <<- gplot + annotate("text", x=batch.xval, label=batch.name,
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
139 size=batch.label.text.size, y=batch.text.height, angle=0,
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
140 color = divide.text.name.color, alpha=batch.text.alpha)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
141 }
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
142 }
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
143 })
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
144 }
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
145
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
146
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
147 plotCells <- function(coldata){
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
148 gplot <<- gplot + geom_point(data=coldata, aes(x=1:nrow(coldata), y=coldata$colsums), pch = 16, cex = 1)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
149 }
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
150
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
151 plotTheme <- function(){
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
152 gplot <<- gplot + theme(plot.title = element_text(hjust = 0.5),
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
153 axis.ticks.x=element_blank(), axis.ticks.y=element_blank(),
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
154 axis.text.x=element_blank()) +
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
155 labs(title=paste("Contamination Plot\n", title), y="Library Size", x="Barcode Index") +
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
156 coord_cartesian(ylim = c(-200, ylim.max)) +
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
157 scale_x_continuous(breaks=NULL)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
158 }
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
159
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
160
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
161 gplot <- ggplot()
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
162 drawPlates(plate.data)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
163 drawBatches(barcode.data, plate.data)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
164 plotCells(coldata)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
165 plotTheme()
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
166
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
167 return(gplot)
253c9448f524 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/crosscontamination_barcode_filter commit f967afe562781e5c8ed4e24e9d1e0bc3ebb29401
iuc
parents: 0
diff changeset
168 }