# NOT RUN {
# Package loading
require(Rmixmod)
# Data loading (categorical data)
data("congress")
# Model-based clustering with 4 components
set.seed(123)
res <- mixmodCluster(congress[,-1], 4, strategy = mixmodStrategy(nbTryInInit = 500, nbTry=25))
# Inference of the parameters used for results visualization
# (specific for Rmixmod results)
# It is better because probabilities of classification are generated
# by using the model parameters
resvisu <- clusvisMixmod(res)
# Component interpretation graph
plotDensityClusVisu(resvisu)
# Scatter-plot of the observation memberships
plotDensityClusVisu(resvisu, add.obs = TRUE)
# }
Run the code above in your browser using DataLab