area plot Archives - Análise Macro

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...