R Сложности с моделированием логистической смеси, logisregmixEM

Question

Для исследования мне нужно применить моделирование смеси, более конкретно: logisregmixEM из пакета mixtools. Я застреваю в приложении на R. Я не могу получить результат логистической регрессии, логарифмической вероятности и пропорций смешивания.

Я провел исследование для него и сумел найти оригинальный код того, как он должен работать, я также нашел «пример». ниже приведены коды. Первый блок - пример того, как они создают двухкомпонентную модель двоичной логистической регрессии. второй блок - это фреймы данных, в которых my y является двоичным с 0 и 1, а x со всеми независимыми переменными. последний блок кодов - это функция logisregmixEM, которая также не работает для меня.

beta <- matrix(c(-10, .1, 20, -.1), 2, 2)
x <- runif(500, 50, 250)
x1 <- cbind(1, x)
xbeta <- x1 %*% beta
w <- rbinom(500, 1, .3)
y <- w*rbinom(500, size = 1, prob = (1/(1 + exp(-xbeta[, 1])))) +
  (1 - w) * 
  rbinom(500, size = 1, prob =(1/(1 + exp(-xbeta[, 2]))))
out.2 <- logisregmixEM(y, x, 
                       beta = beta,
                       lambda = c(.3, .7),
                       verb = TRUE, 
                       epsilon = 1e-03)
out.2



y <- as.numeric(data$Churn) - 1
X <- model.matrix(Churn ~ . , data = data)


logisregmixEM(y, X, 
              N = NULL, 
              lambda = NULL, 
              beta = NULL, 
              k = 2,
              addintercept = TRUE, 
              epsilon = 1e-08, 
              maxit = 10000, 
              verb = FALSE)

Я надеюсь, что вы можете помочь с рабочими и интерпретируемыми результатами, поэтому я нахожу свои коэффициенты логистической регрессии, вероятность логарифма ипропорции смешивания. Любая помощь очень ценится!