Lidando com tendência e heterocedasticidade
Aqui, vamos forçar dados não estacionários a se tornarem estacionários calculando o retorno ou a taxa de crescimento da seguinte forma.
Muitas séries temporais são geradas como $$X_t = (1 + p_t) X_{t-1}$$ o que significa que o valor da série observado no tempo \(t\) é igual ao valor observado no tempo \(t-1\) mais uma pequena variação percentual \(p_t\) no tempo \(t\).
Um exemplo determinístico simples é depositar dinheiro em um banco com juros fixos \(p\). Nesse caso, \(X_t\) é o valor da conta no período \(t\), com um depósito inicial de \(X_0\).
Normalmente, \(p_t\) é chamado de retorno ou taxa de crescimento de uma série temporal, e esse processo costuma ser estável.
Por razões que estão fora do escopo deste curso, pode-se mostrar que a taxa de crescimento \(p_t\) pode ser aproximada por $$Y_t = \log X_t - \log X_{t-1} \approx p_t.$$
No R, \(p_t\) é frequentemente calculado como diff(log(x)), e o gráfico pode ser feito em uma linha com plot(diff(log(x))).
Este exercício faz parte do curso
Modelos ARIMA em R
Instruções do exercício
- Como antes, os pacotes astsa e xts já estão pré-carregados.
- Gere um gráfico multifigura para (1) plotar os dados trimestrais do PIB dos EUA (
gnp) e observar que não é estacionário e (2) plotar a taxa de crescimento aproximada do PIB dos EUA usandodiff()elog(). - Use um gráfico multifigura para (1) plotar os fechamentos diários do DJIA (
djia$Close) e observar que não é estacionário. Os dados são um objetoxts. Em seguida, (2) plote os retornos aproximados do DJIA usandodiff()elog(). Como isso se compara à taxa de crescimento do PIB?
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# astsa and xts are preloaded
# Plot GNP series (gnp) and its growth rate
par(mfrow = c(2,1))
plot(gnp)
# Plot DJIA closings (djia$Close) and its returns
par(mfrow = c(2,1))
plot(djia$Close)