Model-based clustering
----------------------


.. parsed-literal::

        iris = read.table("http://archive.ics.uci.edu/ml/machine-learning-databases/iris/iris.data", 
            sep = ",", header = FALSE)
        names(iris) = c("sepal.length", "sepal.width", "petal.length", "petal.width", 
            "iris.type")
        
        library(mclust)
        iris_mclust = Mclust(iris[, -5], G = c(1:10))
        print(iris_mclust)
     

.. parsed-literal::

        ## 'Mclust' model object:
        ##  best model: ellipsoidal, equal shape (VEV) with 2 components
     

.. parsed-literal::

        plot(iris_mclust, "BIC")
     

.. image:: model_cluster_files/model_cluster_fig_00.png

.. parsed-literal::

        plot(iris_mclust, "classification")
     

.. image:: model_cluster_files/model_cluster_fig_01.png