install.packages("maptools")
require(maptools)
library(maptools)



br<-readShapeSpatial(file.choose())



dados<-read.csv2(file.choose())


attach(dados)
fix(dados)


br.dados<-attr(br, "data")


br.dados$indice<-1:dim(br.dados)[1]


dados$UF<-as.character(dados$UF)
br.dados$UF<-as.character(br.dados$UF)


br.dados<-merge(br.dados, dados)
fix(br.dados)

br.dados<-br.dados[order(br.dados$indice),]


br.dados$idh<-cut(br.dados$idh,8)
summary(br.dados$idh)


palette(c("gray100", "gray90", "gray70", "gray50", "gray30", "gray20", "gray10", "gray1"))


attr(br, "data") <- br.dados
plot(br, col = br$idh)


title("IDH dos Estados Brasileiros")
legend(-80, -11, bty = "n", fill = 8:1, cex = 0.8,legend = levels
(br$idh)[8:1])

points(c(-34.53, -46.37), c(-7.07, -23.32), pch = 23, cex = 0.6,col = "red", bg = "yellow")
text(c(-34.53, -46.37), c(-7.07, -23.32), pos = c(4, 2), cex = 0.6, labels = c("João Pessoa", "São Paulo"))