r3fang / SnapATAC

Analysis Pipeline for Single Cell ATAC-seq
GNU General Public License v3.0
301 stars 125 forks source link

scaleCountMatrix function error in step 13 #47

Open znavidi opened 5 years ago

znavidi commented 5 years ago

Hi,

I am following this sample analysis steps: https://github.com/r3fang/SnapATAC/blob/master/examples/10X_P50/README.md and till the step 13 seems it works fine!

I face this error: Error in .Ops.recycle.ind(e1, len = l2) : vector too long in Matrix - vector operation

while running this command: x.sp = scaleCountMatrix( obj=x.sp, cov=SnapATAC::rowSums(x.sp, mat="bmat"), mat="gmat", method = "RPM" )

Apparently same errors in google are related to the sparse matrix. but there is no specific solution for that. Have you ever seen this error or do you have any idea how it could be solved? because I am running it on a single cell data form 10x Genomics: https://support.10xgenomics.com/single-cell-atac/datasets/1.1.0/atac_v1_pbmc_5k and it seems you have analyzed these datasets before. I would appreciate any help.

Best