Hi,
I wonder is there bug in the Signac code or in data objects, but I get following error:
peaks.to.genes <- LinkPeaks(seurat.object, peak.assay = "peaks", expression.assay = "RNA", n_sample=10, genes.use = G.a[1:10])
Testing 3 genes and 196807 peaks
Error in if (sum(x = links.keep) == 0) { :
missing value where TRUE/FALSE needed
In addition: Warning message:
In cor(x = as.matrix(x = bg.access), y = as.matrix(x = gene.expression), :
the standard deviation is zero
This is really hard to traceback why this error happens? Where exactly it finds Na where logical is needed? That seurat.object is otherwise working fine. And error is not dependent on that gene amount, I just dropped it to very small to quickly confirm repeat of this problem.
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4