EM.clusters: Classify points using EM mixture model
In ClairBee/AS.circular: ArchStats: circular functions

Description Usage Arguments Value Examples

Winner-takes-all classification of points based on a von Mises mixture model obtained using EM.vonmises.

1	EM.clusters(data, model)

`data`	Set of angular data
`model`	Fitted mixture model from `EM.vonmises`.

Vector of cluster numbers

ex1 <- c(rvonmises(120 * 0.3, mu = circular(pi/2), kappa = 10),
         rvonmises(120 * 0.7, mu = circular(pi), kappa = 3))
em1 <- EM.vonmises(ex1, k = 2)
cl <- EM.clusters(ex1, em1)
plot(ex1[cl == 1], stack = T)
points(ex1[cl == 2], col = "blue", stack = T)