# NOT RUN {
data(latrendData)
mKML <- lcMethodKML(response = "Y", id = "Id", time = "Time")
kml1 <- latrend(mKML, nClusters = 1, latrendData)
kml2 <- latrend(mKML, nClusters = 2, latrendData)
kml3 <- latrend(mKML, nClusters = 3, latrendData)
gmm <- latrend(lcMethodLcmmGMM(fixed = Y ~ Time, mixture = ~ Time,
id = "Id", time = "Time"), latrendData)
models <- lcModels(kml1, kml2, kml3, gmm)
subset(models, nClusters > 1 & .method == 'kml')
# }
Run the code above in your browser using DataLab