Hoje de manhã publiquei um código para um gráfico de área utilizando os indexadores da dívida bruta brasileira. O gráfico, entretanto, não estava empilhado, apenas sobreposto. Como a soma da participação dos indexadores tem que dá 100%, um gráfico empilhado talvez fique mais ilustrativo, não é mesmo? Para isso, é preciso fazer algumas modificações no código. Ilustro abaixo.
## Carregar pacotes library(ggplot2) library(XLConnect) library(reshape2) library(xts) ## Importar dados temp = tempfile() download.file('http://www.bcb.gov.br/ftp/notaecon/Partggp.zip',temp) data = unzip(temp, files='Partggp.xls') data = loadWorkbook(data) dbgg = readWorksheet(data, sheet = 1, header = TRUE, colTypes = 'numeric') ## Retirar linhas e colunas que não interessam dbgg = dbgg[, -c(1,2,6,10,15,16,17)] # Retirar colunas desimportantes dbgg = dbgg[complete.cases(dbgg),-1] # Retirar linhas com NA e coluna 1 ## Nomear colunas colnames(dbgg) = c('Cambial Interna', 'Cambial Externa', 'IGP-M', 'IGP-DI', 'IPCA', 'SELIC', 'TJLP', 'TR', 'PRÉ-FIXADO') ## Criar vetor de datas para o gráfico dates = seq(as.Date('2006-12-01'), as.Date('2016-12-01'), by='1 month') ## Ordenar séries conforme vetor de datas e criar novo data frame dbgg = xts(dbgg, order.by=dates) dbgg = data.frame(time = index(dbgg), melt(as.data.frame(dbgg))) colnames(dbgg) = c('time', 'Indexador', 'value') ## Gerar gráfico theme_set(theme_bw()) ggplot(dbgg, aes(x = time, y = value)) + geom_area(aes(colour = Indexador, fill = Indexador))+ xlab('')+ylab('Participação Percentual')+ labs(title='Indexadores da Dívida Bruta brasileira', caption='analisemacro.com.br')
E o resultado vai abaixo...
No nosso Curso de Análise de Conjuntura usando o R os alunos aprendem a coletar, tratar, analisar e apresentar dados macroeconômicos usando o poder do R/RStudio e do Beamer/LaTeX. Saiba mais sobre esse curso inovador clicando no botão abaixo!
Relatórios, apresentações e exercícios macroeconométricos usando extensivamente o R são feitos no âmbito do Clube do Código, o espaço de compartilhamento de códigos da Análise Macro. Ainda não conhece o Clube?! Saiba mais abaixo.