n <- 250 ## create artificial data
set.seed(0)
x <- rbind(matrix(rnorm(n*2), n, 2) %*% diag(c(1,9)),
matrix(rnorm(n*2), n, 2) %*% diag(c(1,9))[,2:1])
xclass <- c(rep(1,n),rep(2,n))
par(pty = "s")
mclust2Dplot(x, classification = xclass, type="classification", ask=FALSE)
odd <- seq(from = 1, to = 2*n, by = 2)
even <- odd + 1
testMclustDA <- mclustDA(trainingData = x[odd, ], labels = xclass[odd],
testData = x[even,])
clEven <- testMclustDA$testClassification ## classify training set
compareClass(clEven,xclass[even])
plot(testMclustDA, trainingData = x[odd, ], labels = xclass[odd],
testData = x[even,])
Run the code above in your browser using DataLab