## ----echo=FALSE, message=FALSE------------------------------------------------
knitr::opts_chunk$set(error=FALSE, message=FALSE, warning=FALSE)
library(BiocStyle)

## -----------------------------------------------------------------------------
library(scRNAseq)
sce <- BachMammaryData(samples="G_2")

set.seed(1000)
sce <- sce[,sample(ncol(sce), 500)]

## -----------------------------------------------------------------------------
set.seed(123)
library(scrapper)
sce <- scrapper::analyze.se(sce)$x
library(scater)
plotTSNE(sce, colour_by="graph.cluster", text_by="graph.cluster")

## -----------------------------------------------------------------------------
library(scDblFinder)
tab <- findDoubletClusters(sce, sce$graph.cluster)
tab

## ----echo=FALSE---------------------------------------------------------------
# Sanity check that one of the clusters is a good doublet candidate.
# If this fails, we probably need to pick a more demonstrative example.
stopifnot(tab[1,"num.de"]==0)

## -----------------------------------------------------------------------------
sessionInfo()