# przyklad 2.1

library(e1071)

attach(irys)
tp<-factor(irys$klasa)
xp<-seq(0.0,2.5,length=100)
yp<-seq(1.0,7.0,length=100)
irysT<-expand.grid(sp=xp,dp=yp)
np<-length(xp)

irys2<-irys[,3:4]

# LDA 
irys.lda<-lda(irys2,tp)
Z<-predict(irys.lda,irysT)
plot(sp,dp,type="n",ylab='Długość płatka',xlab='Szerokość płatka',xlim=c(0, 2.5), ylim=c(1.0,7.0),main="LDA")
text(sp,dp, as.character(tp))
zp<-unclass(Z$class)
zp<-Z$post[,1]-pmax(Z$post[,2],Z$post[,3])
contour(xp,yp,matrix(zp,np),add=T,levels=0.5,labex=0,drawlabels=F)

# QDA 
irys.qda<-qda(irys2,tp)
Z<-predict(irys.qda,irysT)
plot(sp,dp,type="n",ylab='Długość płatka',xlab='Szerokość płatka',xlim=c(0, 2.5), ylim=c(1.0,7.0), main="QDA")
text(sp,dp, as.character(tp))
zp<-unclass(Z$class)
zp<-Z$post[,2]-pmax(Z$post[,3],Z$post[,1])
contour(xp,yp,matrix(zp,np),add=T,levels=0.5,labex=0,drawlabels=F)