#  y:  1=consumidor, 0=não consumidor
#  x1: rendimento anual (em milhares de reais)
#  x2: tamanho do lote (em hectares)
#  x3: 1=se há criação de gado, 0=caso contrário 

dados = read.table("Trator.txt",header=TRUE)

attach(dados)

library(MASS)

plot(x1,x2,col=y+1,pch=16,xlab="X1:rendimento anual (em milhares de reais)",
ylab="X2:tamanho do lote (em hectares)")

mylda   = lda(y~x1+x2)
mypred = predict(mylda)$class
table(mypred,y)

mylda1   = lda(y~x1+x2,CV=TRUE)
table(mylda1$class,y)