## ----echo=FALSE,eval=FALSE----------------------------------------------------
# options(width=80)

## ----eval=FALSE---------------------------------------------------------------
# library(catdata)
# data(knee)

## ----eval=FALSE---------------------------------------------------------------
# knee <- reshape(knee, direction="long", varying=list(5:8), v.names="R",
#                 timevar="Time")
# 
# knee$RD <- rep(0, length(knee$R))
# knee$RD[knee$R>2] <- 1
# 
# knee$Age <- knee$Age - 30
# knee$Age2<-knee$Age^2

## ----eval=FALSE---------------------------------------------------------------
# knee$Th <- as.factor(knee$Th)
# knee$Sex <- as.factor(knee$Sex)

## ----eval=FALSE---------------------------------------------------------------
# library(flexmix)

## ----eval=FALSE---------------------------------------------------------------
# kneeflex2 <-stepFlexmix(cbind(RD,1-RD) ~ 1 | N,	model = FLXMRglmfix(family =
#   "binomial", fixed= ~ Th + Sex + Age + Age2), k = 2, nrep = 5, data = knee)
# 
# kneeflex3 <-stepFlexmix(cbind(RD,1-RD) ~ 1 | N,	model = FLXMRglmfix(family =
#   "binomial", fixed= ~ Th + Sex + Age + Age2), k = 3, nrep = 5, data = knee)
# 
# kneeflex4 <-stepFlexmix(cbind(RD,1-RD) ~ 1 | N,	model = FLXMRglmfix(family =
#   "binomial", fixed= ~ Th + Sex + Age + Age2), k = 4, nrep = 5, data = knee)

## ----eval=FALSE---------------------------------------------------------------
# summary(kneeflex2)@BIC
# summary(kneeflex3)@BIC
# summary(kneeflex4)@BIC

## ----eval=FALSE---------------------------------------------------------------
# summary(refit(kneeflex2))

## ----echo=FALSE,eval=FALSE----------------------------------------------------
# detach(package:flexmix)