#  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)

plot(x1,y,xlab="rendimento anual (em milhares de reais)")

mylogit = glm(y~x1, data = dados, family = "binomial")
summary(mylogit)
par(mfrow=c(1,1))
plot(x1,y,xlab="rendimento anual (em milhares de reais)",ylab="Consumidor (dummy)",pch=16)
points(x1,mylogit$fitted.values,col=2,type="l",lwd=3)
points(x1,mylogit$linear.predictors,col=4,type="l",lwd=3)


mylogit = glm(y~x2, data = dados, family = "binomial")
summary(mylogit)
par(mfrow=c(1,1))
plot(x2,y,xlab="tamanho do lote (em hectares)",ylab="Consumidor (dummy)",pch=16)
points(x2,mylogit$fitted.values,col=2,pch=16)
points(x2,mylogit$linear.predictors,col=4,type="l",lwd=3)


mylogit = glm(y~x1+x2+x3, data = dados, family = "binomial")
summary(mylogit)


mylogit = glm(y~x1+x2, data = dados, family = "binomial")
summary(mylogit)