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.
