Hi,
When I performed LinkPeaks with a single gene name like:
LinkPeaks(
object = seurat.obj,
peak.assay = "ATAC",
expression.assay = "RNA",
genes.use = "INS")
It returned:
No peaks fall within distance threshold
Have you set the proper genome and seqlevelsStyle for ATAC assay?"
I found the problem is caused by line 51 and 57. When genes.use is a single character, the expression.data will become a vector without a row.name.
expression.data <- expression.data[genes.keep, ]
genes <- rownames(x = expression.data)
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