#!/usr/bin/Rscript
library(MASS)

train <- read.table("09_LDA_dataset.txt", header=T)
test <- read.table("09_LDA_newdata.txt", header=T)
colnames(test) <- paste("dim", 1:50, sep="")
plot(rep(2:50, 10), rep(1:10, each=49), xlab="Dimensions", ylab="Prediction")
for(dims in 2:50)
{
	z <- lda(labels ~ ., data = train[,1:(dims+1)])
	print(z)
	pc <- predict(z, test[,1:dims])$class
	points(rep(dims, 10), 1:10, col=as.numeric(pc), pch=19)
}