Я пытаюсь соответствовать 2 компонентной модели гаммы смеси к моим данным (остаточные значения, полученные после выполнения обобщенной линейной модели), используя следующую команду (часть коды):Ошибка в gammamixEM ... Попробуйте другое количество компонентов? в R
expr_mix_gamma <- gammamixEM(expr_glm_residuals, lambda = c(0.75,0.25), k = 2, epsilon = 1e-08, maxit = 1000, maxrestarts=20, verb = TRUE)
код работает для несколько файлов генов (в цикле). он отлично работает для некоторых файлов, тогда как для других он вызывает следующую ошибку:
Error in gammamixEM(expr_glm_residuals, lambda = c(0.75, 0.25), k = 2, : Try different number of components?
Я не могу понять, что происходит. может ли кто-нибудь пролить свет на то же самое? Благодаря