tanaylab / mcATAC

Metacell analysis for ATAC data
https://tanaylab.github.io/mcATAC/
Other
1 stars 0 forks source link

scc_to_peaks malfunction? #63

Closed yonatans2 closed 7 months ago

yonatans2 commented 7 months ago
library(mcATAC)
wd <- '/net/mraid14/export/tgdata/users/yonshap/proj/mmcortex'
setwd(wd)
my_genome <- "mm10"
gset_genome(my_genome)
mca <- readRDS('./output/mcatac/mmcortex_mcatac_feat_peaks.rds')
scc <- scc_read(path = './data/frag_reads_28122022/')

 ✔ Succesfully read a ScCounts object from ./data/frag_reads_28122022/

>  scpeaks <- scc_to_peaks(sc_counts = scc, peaks = dplyr::select(mca@peaks, chrom, start, end, peak_name))
Error in data.frame(i = all_0_idxs, j = 1, x = 0) :
  arguments imply differing number of rows: 0, 1