Displaying 1 result from an estimated 1 matches for "poisonce".
Did you mean:
poison
2011 Feb 28
1
mixture models/latent class regression comparison
...trate the comparison:
library(flexmix)
library(mmlcr)
data(NPreg) # from package flexmix
m1 <- flexmix(yp ~ x, k=2, data=NPreg, model=FLXMRglm(family='poisson'))
NPreg$id <- 1:200 # mmlcr requires an id column
m2 <- mmlcr(outer=~1|id, components=list(list(formula=yp~x,
class="poisonce")), data=NPreg, n.groups=2)
# summary and coefficients for flexmix model
summary(m1)
summary(refit(m1))
# summary and coefficients for mmlcr model
summary(m2)
m2
Regards,
Carson
P.S. I have attached a copy of the mmlcr package with a modified
mmlcr.poisonce function due to errors in the v...