BIC to choose the number of components in a model based clustering using mixtures of von Mises-Fisher distributions
bic.mixvmf(x, G = 5, n.start = 20)
A plot of the BIC values and a list including:
The BIC values for all the models tested.
The run time of the algorithm. A numeric vector. The first element is the user time, the second element is the system time and the third element is the elapsed time.
A matrix containing directional data.
The maximum number of clusters to be tested. Default value is 5.
The number of random starts to try. See also R's built-in function kmeans
for more information about this.
Michail Tsagris.
R implementation and documentation: Michail Tsagris mtsagris@uoc.gr and Giorgos Athineou <gioathineou@gmail.com>.
If the data are not unit vectors, they are transformed into unit vectors.
Hornik, K. and Grun, B. (2014). movMF: An R package for fitting mixtures of von Mises-Fisher distributions. Journal of Statistical Software, 58(10):1--31.
mixvmf.mle, rmixvmf, mixvmf.contour
x <- as.matrix( iris[, 1:4] )
x <- x / sqrt( rowSums(x^2) )
bic.mixvmf(x)
Run the code above in your browser using DataLab