diazlab / CONICS

CONICS: COpy-Number analysis In single-Cell RNA-Sequencing
73 stars 28 forks source link

error on plotall #28

Closed FuqianShi closed 3 years ago

FuqianShi commented 4 years ago

Hello, I use mm10 region, and prepared the matrix

F1_AAACCCATCTCTCTAA.1 F1_AAACGAATCAGATGCT.1 F1_AAACGCTGTACTCCGG.1 0610005C13Rik 0 0 0 0610009B22Rik 0 0 0 0610009O20Rik 0 0 0 0610010F05Rik 0 0 0 0610010K14Rik 0 0 0

normFactor success

when plot l=plotAll(suva_expr,normFactor,regions,gene_pos,"SUVA_CNVs")

I got an error: [1] "Fitting GMM for chr1 0:195471971 iteration 1" number of iterations= 33 [1] "Fitting GMM for chr1 0:195471971 iteration 2" number of iterations= 42 [1] "Fitting GMM for chr1 0:195471971 iteration 3" number of iterations= 42 [1] "Fitting GMM for chr1 0:195471971 iteration 4" number of iterations= 44 [1] 19 [1] "Fitting GMM for chr2 0:182113224 iteration 1" [1] "EM algorithm did not converge for region 2 0 182113224" [1] "Fitting GMM for chr2 0:182113224 iteration 2" [1] "EM algorithm did not converge for region 2 0 182113224" [1] "Fitting GMM for chr2 0:182113224 iteration 3" [1] "EM algorithm did not converge for region 2 0 182113224" [1] "Fitting GMM for chr2 0:182113224 iteration 4" [1] "EM algorithm did not converge for region 2 0 182113224" Error in hist.default(chr_exp, breaks = 50, main = paste("Chr: ", chr, : character(0)

Would you please help to fix the problem for me? Thank you!

MelissaSaichi commented 3 years ago

Hi @FuqianShi , I am actually getting the same error as you. Did you manage to solve this issue ?