Указание смеси означает при использовании regmixEM или flemix в R
В настоящее время я пытаюсь приспособить модели регрессии смеси к своим данным. К сожалению, алгоритм не сходится при k>2. Поэтому я хотел бы дать лучшие начальные значения вместо тех, которые выбирают R. Тем не менее, синтаксис:
regmixEM(y, x, lambda = NULL, beta = NULL, sigma = NULL, k = 2,
addintercept = TRUE, arbmean = TRUE, arbvar = TRUE,
epsilon = 1e-08, maxit = 10000, verb = FALSE)
указание mu не сработало, поскольку это не тот параметр, который я могу дать. Есть ли способ указать параметры? Я знаю, что есть пакет под названием flemix. Синтаксис функции flemix:
flexmix(formula, data = list(), k = NULL, cluster = NULL,
model = NULL, concomitant = NULL, control = NULL,
weights = NULL)
## S4 method for signature 'flexmix'
summary(object, eps = 1e-4, ...)
Читая руководство, я могу найти только способ манипулирования вероятностями членства, но не мю.
Есть ли способ заставить алгоритм сходиться?
Я благодарен за любую помощь, спасибо.