comparison probecoverage.r @ 4:daec4df60281 draft

planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/probecoverage commit 2b4a222804a40730abdaec6fb6aefc2295e89e11
author artbio
date Sun, 24 Sep 2017 18:10:12 -0400
parents 4f744d3aaf0b
children bea8435e1e79
comparison
equal deleted inserted replaced
3:4f744d3aaf0b 4:daec4df60281
38 } 38 }
39 scaleFUN <- function(x) sprintf("%.3f", x) 39 scaleFUN <- function(x) sprintf("%.3f", x)
40 40
41 ## end of function 41 ## end of function
42 ## let's do a dataframe before plotting 42 ## let's do a dataframe before plotting
43 maxdepth <- trunc(max(Table[,5:length(Table)]/(Table$end-Table$start))) + 20 43 if (args$method == 'bedtools') {
44 maxdepth <- trunc(max(Table[,5:length(Table)]/(Table$end-Table$start))) + 20
45 } else {
46 maxdepth <- trunc(max(Table[,5:length(Table)])) + 20
47 }
48
44 graphpoints <- data.frame(1:maxdepth) 49 graphpoints <- data.frame(1:maxdepth)
45 i <- 5 50 i <- 5
46 for (colonne in colnames(Table)[5:length(colnames(Table))]) { 51 for (colonne in colnames(Table)[5:length(colnames(Table))]) {
47 graphpoints <- cbind(graphpoints, mapply(cumul, 1:maxdepth, rep(i, maxdepth))) 52 graphpoints <- cbind(graphpoints, mapply(cumul, 1:maxdepth, rep(i, maxdepth)))
48 i <- i + 1 53 i <- i + 1