我正在尝试使用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步骤驱动程序并不相同。但这是否是我需要做的事情，只需按照https://github.com/cran/flexmix/blob/master/R/flexmix.R中所述调整一些代码并编写一些底层函数的变体来编写一个新的flexmix函数？
我是不熟悉S4对象系统的新手，因此，我将不胜感激。