#==========================================================================================
# Aula 6 da disciplina ce223 (22/03/2011)
# Estatísticas descritivas aplicadas por colunas e linhas
#                                                               Professor Walmes M. Zeviani
#                                                                     www.leg.ufpr.br/ce223
#==========================================================================================

#------------------------------------------------------------------------------------------
# ipeadata > social > temas > eduação > analfabetos - pessoas 15 anos e mais > município

analf <- read.table("analf.txt", header=TRUE, sep="\t")
str(analf)

#------------------------------------------------------------------------------------------
# calcular o número de dados para cada ano

names(analf) <- c("sig","cod","mun","ano1970","ano1980","ano1991")
str(analf)

#------------------------------------------------------------------------------------------
# calcular a média para cada ano

apply(analf[,4:6], 2, mean, na.rm=TRUE)
apply(analf[,4:6], 2, var, na.rm=TRUE)

mmun <- apply(analf[,4:6], 1, mean, na.rm=TRUE)

#------------------------------------------------------------------------------------------
# qual o melhor índice educacional do 1 ano?

apply(analf[,4:6], 2, max, na.rm=TRUE)
z <- apply(analf[,4:6], 2, which.max)
z

analf[z,]

w <- apply(analf[,4:6], 2, which.min)
analf[w,]

#------------------------------------------------------------------------------------------
# qual o município de menor média?

str(mmun)
analf[which.min(mmun),]

#------------------------------------------------------------------------------------------
# e se eu aplicar uma função que não retorne apenas um número?

z <- apply(analf[,4:6], 2, range, na.rm=TRUE)
str(z)

#------------------------------------------------------------------------------------------

z <- apply(analf[,4:6], 2, summary, na.rm=TRUE)
z
str(z)

#------------------------------------------------------------------------------------------
# como fazer um gráfico para cada um dos anos, um histograma?

par(mfrow=c(1,3))
apply(analf[,4:6], 2, hist)

#------------------------------------------------------------------------------------------