#####################################################################################
#
# Analfabetos - pessoas 15 anos e mais (2000)	
# Razao entre a renda dos 10% mais ricos e 40% mais pobres (2000)	
#
#####################################################################################
rm(list=ls())

nomes = c("Analfabetismo - pessoas de 15 ou mais anos",
          "Razao renda - 10 % mais ricos e 40% mais pobres (LOG)")

data  = read.table("analfabetismo-rendapercapita.txt",header=TRUE)

y     = log(data[,3])

x     = data[,2]

reg      = lm(y~x)
ehat     = reg$res
ehat2    = ehat^2
logehat2 = log(ehat2)

pdf(file="analfabetismo-rendapercapita.pdf",width=10,height=8)

par(mfrow=c(1,1))
plot(x,y,xlab=nomes[1],ylab=nomes[2])
abline(reg$coef,col=2,lwd=2)

U = max(abs(ehat))
L = -U
plot(x,ehat,xlab=nomes[1],ylab="Residuals",ylim=c(L,U))
abline(h=0,col=2,lwd=2)

plot(x,ehat2,xlab=nomes[1],ylab="Square residuals")
abline(lm(ehat2~x),col=2,lwd=3)

plot(x,logehat2,xlab=nomes[1],ylab="Log square residuals")
abline(lm(logehat2~x),col=2,lwd=3)

dev.off()