Ampliando a série de desemprego da PNAD Contínua com o R - Análise Macro
Ampliando a série de desemprego da PNAD Contínua com o R
By Vitor Wilher 14 de junho de 2017 No Comments
Desde o ano passado, os economistas brasileiros estão orfãos da Pesquisa Mensal de Emprego do IBGE. Em seu lugar, o instituto tem divulgado mensalmente a PNAD Contínua com base em trimestres móveis. O problema nessa mudança é que a série da PNAD é curta, o que impede esses mesmos economistas de fazer algum exercício mais sofisticado com os dados. Esse, a propósito, é um dos problemas que temos ajudado a resolver com os cursos aplicados da Análise Macro. Vamos mostrar aqui como podemos usar o R para estender a amostra da taxa de desemprego da PNAD com base na taxa de desemprego da PME, tendo assim uma série com mais observações.
Abaixo levamos os dados para o R:
pme = ts(read.csv2('pme.csv', header=T, sep=';', dec=',')[,-1],
start=c(2002,03), freq=12)
pnad = BETS.get(24369)
data = cbind(pme,pnad)
Uma vez importados os dados, observamos como esperado que as mesmas possuem tamanhos diferentes. Um primeiro ajuste, portanto, é igualar a amostra das séries. Isso é feito abaixo.
desemprego = ts(data[complete.cases(data),], start=c(2012,03), freq=12)
Abaixo um gráfico com as séries.
O gráfico nos mostra que a série da PNAD Contínua tem uma média maior do que a da PME. Ambas as séries caminham, entretanto, de forma parecida, a despeito da diferença da amostra. Abaixo as estatísticas descritivas das duas séries.
Taxa de Desemprego: PNAD Contínua vs. PME.
pme 48 5,73 1,01 4,30 8,20
pnad 48 7,51 0,89 6,20 10,20
Uma vez definido o conjunto de dados, passemos ao exercício em si. Como o nosso objetivo é basicamente expandir a taxa de desemprego da PNAD Contínua, vamos então regredir a mesma contra a taxa da PME.
reg = lm(pnad~pme, data=desemprego)
A tabela 2 resume a regressão.
PNAD Contínua vs. PME
pme 0.842***
Constant 2.682***
Adjusted R2 0.919
Residual Std. Error 0.253 (df = 46)
F Statistic 531.093*** (df = 1; 46)
Em outras palavras, a taxa de desemprego da PNAD estimada e a da PME se relacionam pela equação
Abaixo, plotamos a série original da taxa de desemprego da PNAD Contínua contra a taxa gerada pelo modelo.
Uma vez estimado o modelo, podemos agora utilizá-lo para construir nossa série de desemprego da PNAD Contínua para trás, com base na taxa de desemprego da PME. Isso é feito com o código abaixo.
### Criar Série PNAD Contínua Ampliada
pnad.pme = ts(data[!complete.cases(data[,2]),], start=c(2002,03), freq=12)
pnad.pme[,2] = coef(reg)[1] +coef(reg)[2]*pnad.pme[,1]
data[1:nrow(pnad.pme),2] = pnad.pme[1:nrow(pnad.pme),2]
pnad = data[,2]
E agora a taxa de desemprego que construímos.
Podemos, agora, dessazonalizar a nossa taxa de desemprego, utilizando o pacote seasonal e colocar abaixo as duas séries.
Temos agora duas séries, uma normal e outra dessazonalizada, que podemos utilizar em exercícios futuros. Observe, entretanto, que a mesma é apenas uma aproximação estatística da PNAD Contínua através da PME, ok? Até o próximo exercício e não esqueça de dar uma olhada nos cursos aplicados de R da Análise Macro! 🙂