when plot l=plotHistogram(l,suva_expr,clusters=2,zscoreThreshold=4,mice)
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!
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=plotHistogram(l,suva_expr,clusters=2,zscoreThreshold=4,mice)
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!