library(scatterplot3d)
library(mva)
library(MASS)
library(fpc)

N<-20
c1<-cbind(	
			rep(1,N),
			rnorm(N, sd=5, mean=0),	
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=0.1, mean=1) 
		)

c2<-cbind(	rep(2,N),
			rnorm(N, sd=5, mean=0),		
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=5, mean=0),
			rnorm(N, sd=0.1, mean=-1)
		)

cx<-data.frame(rbind(c1,c2))

names(cx)<-c(	"y",
				"x1",
				"x2",
				"x3",
				"x4",
				"x5",
				"x6",
				"x7",
				"x8",
				"x9",
				"x10"
			)

par(ask=T)

plot(cx[,-1],col=cx[,1])

pca<-prcomp(cx[,-1])

v<-data.frame(pca$x)

names(v)<-c(	"PC1",
				"PC2",
				"PC3",
				"PC4",
				"PC5",
				"PC6",
				"PC7",
				"PC8",
				"PC9",
				"PC10"
			)

plot(v,col=cx[,1])

X<-discrcoord(cx[,-1],cx[,1])
d<-data.frame(X$proj)

names(d)<-c(	"LD1",
				"LD2",
				"LD3",
				"LD4",
				"LD5",
				"LD6",
				"LD7",
				"LD8",
				"LD9",
				"LD10"
			)

plot(d, col=cx[,1])

par(ask=F)