Comments (6)
Hi,
Is Parse the combinatorial indexing method (i.e. what used to be called split-seq before it was commercial)?
If they did some assays with different amounts of loading you can perhaps use that to estimate an expected rate, but if not it's not a big deal, in a typical context the importance of the expected doublet rate. You can leave dbr
to the default and set dbr.sd=1
(i.e. very high uncertainty of the dbr
), and basically those parameters won't have an impact, and the threshold will be set based on the misclassification rate.
You can also have a look at the histogram of the doublet score, unless you have a very complex dataset one typically observes a very strong bimodal distribution, so that it doesn't matter much where exactly one sets the threshold in between.
Hope this helps,
Pierre-Luc
from scdblfinder.
Hi,
Yes. It is four-step barcoding to label cells in different wells.
I left dbr
default and this ended up with ~2600 doublet but once I set up dbr = 0.03
, it ended up with ~900 doublet because I thought that Parse has much lower doublet rate 10x has and I should set up smaller dbr
.
I'll draw a histogram to see distribution of doublet rate, thank you for the suggestion.
Best.
~Kaan
from scdblfinder.
Right, you probably do have a lot more cells than in a single 10x capture, right?
Try with a high dbr.sd
, and let me know about your histogram of doublet scores!
from scdblfinder.
The number of cells is 12k in data for now but data with ~ 1 million cells will be produced soon, so yes, number of cell will be way higher in Parse data compare to 10x.
Once dbr.sd = 1
set in scDblFinder()
function, 2906 out of 12k cells detected as doublet which is much closer to what I expected.
Here, density of doublet score across cells:
It seems that doublet score is fairly discriminative between singlets and doublets by regarding histogram plot of doublet score. However, when I put score threshold of 0.598 manually, this ended up with almost the same number of doublets (2908) as what I found automatically with dbr.sd = 1
. Is that threshold of 0.598 feasible in doublet detection?
Best.
~Kaan
from scdblfinder.
I don't understand your question: the dbr
and dbr.sd
(along with the misclassification rate) are used to establish the threshold. Of course if you manually use the threshold one expects to find the same thing (I assume the difference is due to rounding)...
You've got very clear singlets and very clear doublets, and some ambiguous cells in the middle, which will be influenced by the thresholding. Where exactly you should place the threshold also depends on how bad missing doublets or singlets would be given your aims, but if you don't have strong opinion it's reasonable to leave it where the thresholding procedure set it.
from scdblfinder.
My expectation was higher threshold in detection of doublets, e.g., around 0.9 but I'm surprised that it was ~0.6. I agree on difference in threshold due to rounding.
I see, I do not have strong opinion so I'll use dbr.sd = 1
and whatever threshold set up will be considered in doublet detection.
Thank you!
from scdblfinder.
Related Issues (20)
- Error running scDblFinder HOT 4
- how to get layers info after running scDblFinder HOT 2
- Multiple errors: samples argument and Matrix HOT 4
- Removing genes from count matrix before doublet detection HOT 3
- SCTransform example in vignette not working HOT 1
- Wild Variability in doublet scores HOT 2
- Running scDblFinder before or after normalization? HOT 2
- Error running scDblFinder HOT 3
- non-interpretable error message when providing non-integer as nfeatures argument HOT 1
- Clarify recommended nFeatures/artificialDoublets for scATACseq doublet removal in vignette HOT 2
- Unable to run scDblFinder - "as_cholmod_sparse" does not exist HOT 4
- Unable to install scDblFinder HOT 2
- Ambient RNA Removal HOT 4
- Compatibility between Seurat, matrix and scDblFinder versions HOT 4
- error in scATAC HOT 1
- Does ‘Size factors should be positive’ error matter? HOT 4
- "did not converge" Error on cellbender3 HOT 8
- How can I add results to back my original file HOT 2
- Error running scDblFinder HOT 7
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from scdblfinder.