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

N<-200
c1<-cbind(rep(0,N),rnorm(N, sd=5, mean=0),rnorm(N, sd=0.1, mean=1), rep(1,N))
c2<-cbind(rep(0,N),rnorm(N, sd=5, mean=0),rnorm(N, sd=0.1, mean=-1), rep(2,N))
cx<-data.frame(rbind(c1,c2))
names(cx)<-c("x1","x2","x3","y")

par(mfcol=c(1,1),ask=T)

#scatterplot3d(cx[,-4],color=cx[,4])

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

pca<-prcomp(cx[,-4])
v<-data.frame(pca$x[,1],pca$x[,2],pca$x[,3])
names(v)<-c("PC1","PC2","PC3")
plot(v,col=cx[,4])

X<-discrcoord(cx[,-4],cx[,4])
data<-data.frame(X$proj)
names(data)<-c("LD1","LD2","LD3")
plot(data, col=cx[,4])

par(ask=F)