estados = c("Rondônia","Acre","Amazonas","Roraima","Pará","Amapá","Tocantins",
"Maranhão","Piauí","Ceará","Rio Grande do Norte","Paraíba","Pernambuco","Alagoas",
"Sergipe","Bahia","Minas Gerais","Espírito Santo","Rio de Janeiro","São Paulo","Paraná",
"Santa Catarina","Rio Grande do Sul","Mato Grosso do Sul","Mato Grosso","Goiás","Distrito Federal")

data = read.csv("IDHM-Brasileiro-dados.csv",header=TRUE)

n = 5565

data = data[(n+1):(2*n),]

est = unique(data[,2])

cbind(estados,est)

nmunicip = NULL
for (estado in est)
  nmunicip = c(nmunicip,sum(data[,2]==estado))
  
cbind(estados,est,nmunicip)

cbind(estados,table(data[,2]))