queimadas · R

Queimadas e Quebras Estruturais

Meu amigo Cristiano Oliveira, da FURG, fez um ótimo texto sobre as queimadas usando boa econometria na rede social que não mais uso. Podemos aproveitar sua análise para brincar um pouco com o R.

Primeiro, é possível ver que os dados apresentam sazonalidade marcante.

Rplot01

As médias sazonais?

sazon

Como qualquer um pode perceber, as queimadas aumentam sempre no mês de Setembro. Dito isso, algo displicentemente, busquei quebras endógenas na série. Eis o que encontrei.Rplot

O algoritmo encontrou duas quebras estruturais na série (pensando em termos de média), mas os intervalos de confiança não me inspiram confiança. As quebras seriam algo como Junho de 2002 e Agosto de 2007 seriam as supostas quebras.

Contudo, buscar quebras estruturais em séries sazonais não é algo trivial como aponta este artigo e este pequeno post.

A propósito, agradeço ao Cristiano pelo ótimo post. É inspirador ver meus colegas ajudando na melhoria da qualidade do debate público.

Anúncios