r3fang / SnapATAC

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

ARPACK solver failed to converge in runDiffusionMaps #145

Closed wshands closed 4 years ago

wshands commented 4 years ago

I have created a pipeline based on the example https://github.com/r3fang/SnapATAC/blob/master/examples/10X_brain_5k/README.md to process my sample data.

I am not familiar with dimensionality reduction using Eigen vectors. I have encountered the following error at this point in the code, which corresponds to the step 5 in the example:

Code: dim(x.sp@bmat);

x.sp = runDiffusionMaps( obj=x.sp, input.mat="bmat", num.eigs=50 );

Output: [1] 758 199883

Epoch: computing eigen decomposition ... Error in igraph::arpack(f, extra = M, sym = sym, options = list(which = wh, : At arpack.c:1174 : ARPACK error, Maximum number of iterations reached Calls: runDiffusionMaps -> runEigDecomp -> eig_decomp -> -> .Call In addition: Warning message: In igraph::arpack(f, extra = M, sym = sym, options = list(which = wh, : At arpack.c:776 :ARPACK solver failed to converge (1001 iterations, 0/51 eigenvectors converged) Execution halted

How can I fix this? What should I focus on in my dataset to help fix the issue? What other information do you need?

Thanks!