########################################################
# baixar dados de fecund. para munic. bras. em 2010
# visualizar a distrib. de taxas de fecund. para
# mulheres de 20-24 anos, para o Brasil inteiro e 
# para a mesoregião Central Potiguar
########################################################

# local dos dados na internet
url = 'http://ufrn.schmert.net/dados/fecund.2010.csv'

# baixando os dados para um data.frame
D = read.csv( url)

# baixar o pacote dplyr (é necessario fazer isso só uma vez)
install.packages('dplyr')  

# use o pacote
library('dplyr')

# estimativas 5f20 para todos os municipios
tx20 = D$Cr20 / D$Mu20
head(tx20)

# distribuicao de estimativas para todos os municipios 
plot( density( tx20), main='Todos os munic. bras.')

# filtrar e selecionar D: só municipios na mesoreg. 'Central Potiguar'
#  e só variáveis cujos nomes contêm 'Mu' ou 'Cr'

CentPotig = filter( D, Nome_Meso == 'Central Potiguar')
CentPotig = select( tmp, contains('Mu'), contains('Cr'))

# estimativas 5f20 para todos os municipios
tx20 = CentPotig$Cr20 / CentPotig$Mu20
head(tx20)

# distribuicao de estimativas para todos os municipios 
plot( density( tx20), main ='Munic. na mesoregião Central Potiguar')

# comentário do professor: também pode-se fazer tudo nas linhas 30-38
# assim no dplyr

CentPotig = filter( D, Nome_Meso == 'Central Potiguar') %>%
              select(contains('Mu'), contains('Cr')) %>%
              mutate( tx20 = Cr20/Mu20)

# distribuicao de estimativas para todos os municipios 
plot( density( CentPotig$tx20), main ='Munic. na mesoregião Central Potiguar')