Mercurial > repos > iuc > dropletutils
comparison scripts/dropletutils.Rscript @ 8:a9caad671439 draft default tip
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/dropletutils/ commit 9d465ee72398a8a74fe499a2e4b74f507a5845cc"
| author | iuc |
|---|---|
| date | Mon, 05 Jul 2021 13:40:00 +0000 |
| parents | 2c1200fba922 |
| children |
comparison
equal
deleted
inserted
replaced
| 7:2c1200fba922 | 8:a9caad671439 |
|---|---|
| 159 plot(brout$rank, brout$total, log = "xy", | 159 plot(brout$rank, brout$total, log = "xy", |
| 160 xlab = "(log) Rank", ylab = "(log) Total UMI Counts") | 160 xlab = "(log) Rank", ylab = "(log) Total UMI Counts") |
| 161 o <- order(brout$rank) | 161 o <- order(brout$rank) |
| 162 lines(brout$rank[o], brout$fitted[o], col = "red") | 162 lines(brout$rank[o], brout$fitted[o], col = "red") |
| 163 | 163 |
| 164 abline(h = brout$knee, col = "dodgerblue", lty = 2) | 164 abline(h = metadata(brout)$knee, col = "dodgerblue", lty = 2) |
| 165 abline(h = brout$inflection, col = "forestgreen", lty = 2) | 165 abline(h = metadata(brout)$inflection, col = "forestgreen", lty = 2) |
| 166 legend("bottomleft", lty = 2, col = c("dodgerblue", "forestgreen"), | 166 legend("bottomleft", lty = 2, col = c("dodgerblue", "forestgreen"), |
| 167 legend = c("knee", "inflection")) | 167 legend = c("knee", "inflection")) |
| 168 dev.off() | 168 dev.off() |
| 169 | 169 print(paste("knee =", metadata(brout)$knee, ", inflection = ", metadata(brout)$inflection)) |
| 170 print(paste("knee =", brout$knee, ", inflection = ", brout$inflection)) | |
| 171 } | 170 } |
| 172 | 171 |
| 173 ## Main | 172 ## Main |
| 174 set.seed(seed.val) | 173 set.seed(seed.val) |
| 175 | 174 |