Mercurial > repos > iuc > dropletutils
comparison scripts/dropletutils.Rscript @ 7:2c1200fba922 draft
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/dropletutils/ commit 0a1586a00f94215b8182423e4b2f34fdf706a47d"
| author | iuc |
|---|---|
| date | Thu, 07 Jan 2021 19:36:54 +0000 |
| parents | 8855361fcfc5 |
| children | a9caad671439 |
comparison
equal
deleted
inserted
replaced
| 6:8855361fcfc5 | 7:2c1200fba922 |
|---|---|
| 155 | 155 |
| 156 brout <- do.call(barcodeRanks, c(bparams)) | 156 brout <- do.call(barcodeRanks, c(bparams)) |
| 157 | 157 |
| 158 png(files$plot) | 158 png(files$plot) |
| 159 plot(brout$rank, brout$total, log = "xy", | 159 plot(brout$rank, brout$total, log = "xy", |
| 160 xlab = "(log) Rank", ylab = "(log) Total Number of Barcodes") | 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 = brout$knee, col = "dodgerblue", lty = 2) |
| 165 abline(h = brout$inflection, col = "forestgreen", lty = 2) | 165 abline(h = brout$inflection, col = "forestgreen", lty = 2) |