Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Issue with knownDoublets #45

Closed
bjstewart1 opened this issue Jul 19, 2021 · 3 comments
Closed

Issue with knownDoublets #45

bjstewart1 opened this issue Jul 19, 2021 · 3 comments

Comments

@bjstewart1
Copy link

bjstewart1 commented Jul 19, 2021

the knownDoublets option in scDblFinder is is throwing an error when presented with knownDoublets and also samples

Reproducible example:

library(scDblFinder)
sce <- mockDoubletSCE()
sce$type <- sce$type %in% "doublet"
sce$channel <- c(rep("sample1", floor(ncol(sce)/2)), rep("sample2", ceiling(ncol(sce)/2)))[1:ncol(sce)]
scldbl <- scDblFinder(sce = sce, 
                      samples = "channel",
                      knownDoublets = "type")


yields
Error in value[[3L]](cond) : An error occured while processing sample 'cluster1': Error in .checkColArg(sce, knownDoublets):knownDoubletsshould have a length equal to the number of columns insce.

but

scldbl <- scDblFinder(sce = sce, 
                      #samples = "channel",
                      knownDoublets = "type")


succeeds
Clustering cells... 4 clusters Creating ~5000 artifical doublets... Dimensional reduction Finding KNN... Evaluating cell neighborhoods... Training model... Finding threshold... Threshold found:0.425 19 (3.7%) doublets called

@plger
Copy link
Owner

plger commented Jul 19, 2021

thanks, will look into it in the next couple of days. I'm assuming you're using the latest github version?
Note that the current usage of knownDoublets actually decreases performance if (as is often the case) a sizeable fraction of the known doublets are homotypic. More flexibility on their usage will be added soon.

@bjstewart1
Copy link
Author

yes this is after
BiocManager::install("plger/scDblFinder")
on 19-7-21
thanks

@plger plger closed this as completed in 6b80db4 Jul 20, 2021
@plger
Copy link
Owner

plger commented Jul 20, 2021

should be fixed in the latest push

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants