ComeçarComece de graça

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

Ver curso

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 usando diff() e log().
  • 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 objeto xts. Em seguida, (2) plote os retornos aproximados do DJIA usando diff() e log(). 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)


 
Editar e executar o código