Skip to content

Commit 76a88c5

Browse files
committed
Accommodate changes from DropletUtils::emptyDrops()
- Details of changes: MarioniLab/DropletUtils#118 - Since the default is now `alpha=Inf`, the claims around `alpha` values and sanity check no longer make sense, so just remove them.
1 parent c60ce4b commit 76a88c5

File tree

1 file changed

+0
-2
lines changed

1 file changed

+0
-2
lines changed

inst/book/droplet-processing.Rmd

Lines changed: 0 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -303,7 +303,6 @@ hist(log10(e.out.hto$Total[is.cell.hto]), xlab="Log[10] HTO count", main="")
303303

304304
While both approaches are valid, we tend to favor the cell calls derived from the gene matrix as this directly indicates that a cell is present in the droplet.
305305
Indeed, at least a few libraries have very high total HTO counts yet very low total gene counts (Figure \@ref(fig:hto-total-comp)), suggesting that the presence of HTOs may not always equate to successful capture of that cell's transcriptome.
306-
HTO counts also tend to exhibit stronger overdispersion (i.e., lower `alpha` in the `emptyDrops()` calculations), increasing the risk of violating `emptyDrops()`'s distributional assumptions.
307306

308307
```{r hto-total-comp, fig.cap="Total HTO counts plotted against the total gene counts for each library in the cell line mixture dataset. Each point represents a library while the dotted lines represent the thresholds below which libraries were assumed to be empty droplets."}
309308
table(HTO=is.cell.hto, Genes=is.cell, useNA="always")
@@ -316,7 +315,6 @@ abline(h=metadata(e.out.hto)$lower, col="blue", lwd=2, lty=2)
316315

317316
```{r, echo=FALSE}
318317
# Sanity checks for trash-talk above.
319-
stopifnot(metadata(e.out.hto)$alpha < 50 * metadata(e.out.gene)$alpha)
320318
stopifnot(sum(e.out.gene$Total < 100 & e.out.hto$Total > 1000) > 50)
321319
```
322320

0 commit comments

Comments
 (0)