Nothing
library(EMCluster, quietly = TRUE)
library(MASS, quietly = TRUE)
set.seed(1234)
### Crabs.
x <- as.matrix(crabs[, 4:8])
ret <- init.EM(x, nclass = 4, min.n = 40, min.n.iter = 100)
ret.proj <- project.on.2d(x, ret)
### Plot.
pdf("crabs_ppcontour.pdf", height = 5, width = 5)
par(mar = rep(0.1, 4))
plotppcontour(ret.proj$da, ret.proj$Pi, ret.proj$Mu, ret.proj$S,
ret.proj$class, angle = pi/6, main = "Crabs K = 4")
dev.off()
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.