Я пытаюсь настроить конечную регрессию смеси, используя flexmix.Я хотел бы настроить каждый компонент для использования различных объясняющих переменных.
Например, вот пример из справки flexmix:
data("NPreg", package = "flexmix")
ex1 <- flexmix(yn ~ x + I(x^2), data = NPreg, k = 2)
matplot(NPreg$x, fitted(ex1), pch = 16, type = "p")
points(NPreg$x, NPreg$yn)
В этом случае я хотел бы настроить один компонент для оценки:
yn ~ x #k=1
, в то время как другой использовать:
yn ~ x + I(x^2) #k=2
Я считаю, что это не то же самое, что написать новый M-step драйвер.Но это то, что мне нужно сделать, написав новую функцию flexmix, настроив часть кода, как описано в https://github.com/cran/flexmix/blob/master/R/flexmix.R, и написав несколько вариантов базовых функций?
Я новичок в пониманииОбъектные системы S4, поэтому я буду признателен за любые комментарии или предложения.