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 |
