glm
FLXglm(formula = . ~ ., family = c("gaussian", "binomial", "poisson", "Gamma"))
flexmix
update.formula
FLXmodel