d    <- read.table("inp_iris.txt",header=TRUE)
ii   <- seq(1,150,2)
Ldat <- d[ii,1:4]
Lgrp <- d[ii,5]
Edat <- d[-ii,1:4]
Egrp <- d[-ii,5]

library(MASS)
res1 <-lda(Ldat,Lgrp)
res2 <-predict(res1)
plot(
  res2$x,
  xlim=c(-10,10),
  ylim=c(-3,3),
  pch=c(0,1,2)[unclass(Lgrp)],
  col=c("maroon","blue","red")[unclass(Lgrp)]
  )
y <- predict(res1,Edat)
points(
  y$x,
  pch=c(22,21,24)[unclass(Lgrp)],
  bg=c("maroon","blue","red")[unclass(Egrp)]
  )