library(MASS)
library(mva)

rawdata<-matrix(scan("tab1_1.dat"),ncol=3,byrow=T)
group <- rawdata[,1]
X <- 100 * rawdata[,2:3]
d <- data.frame(rawdata)
names(d)<-c("y","x1","x2")
d$x1 = d$x1 * 100
d$x2 = d$x2 * 100

v<-prcomp(X)

par(mfcol=c(2,1))
plot(data.frame(v,rep(0,length(v[1]))), col=group+1)

g<-lda( y ~ x1 + x2, data=d)
v2 <- predict(g, d)

par(mfcol=c(3,1))

plot(X, type="n")
text(X[group==0,], col=1)
text(X[group==1,], col=2)

plot(data.frame(v$x,rep(0,length(v$x[1]))), type="n")
text(matrix(c(v$x[group==0], rep(0, length(v$x[group==0]))), ncol=2), col=1)
text(matrix(c(v$x[group==1], rep(0, length(v$x[group==1]))), ncol=2), col=2)
title("Projection onto 1st PC")
mtext(paste("beta_1=",format(b[1],digits=2), "beta_2=",format(b[2],digits=2)))

plot(data.frame(v2$x,rep(0,length(v[1]))), type="n")
title("Projection onto 1st Canonical (LDA)")
mtext(paste("beta_1=",format(g$sc[1],digits=2), "beta_2=",format(g$sc[2],digits=2)))
text(matrix(c(v2$x[group==0], rep(0, length(v2$x[group==0]))), ncol=2), col=1)
text(matrix(c(v2$x[group==1], rep(0, length(v2$x[group==1]))), ncol=2), col=2)