Mercurial > repos > artbio > probecoverage
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 |