Modelo Café Arábica - Mínimos Quadrados Ordinários

UNIVERSIDADE FEDERAL DA PARAÍBA
PROFESSOR: Sinézio Fernandes Maia

Autor

Josué de Meneses Lopes

Data de Publicação

26 de abril de 2024

Atividade em Mínimos Quadrados Ordinários (MQO)

A atividade abaixo se refere a busca de um modelo robusto e não tendencioso em MQO, para encontrar a quantidade ótima de contratos de hedge para o mercado futuro de café arábica.

Configurações iniciais

Para evitar notações científicas:

options(scipen = 9999)
options(max.print = 100000)

Pacotes útilizados:

library(readxl)
library(tidyverse)
library(deflateBR)
library(dynlm)
library(ggplot2)
library(fBasics)
library(lmtest)
library(whitestrap)
library(FinTS)

Base de dados

dados <- read_excel("/Users/josue/Documents/R/Derivativos/modelo mqo/SJCFUT (Mensal) (1).xlsx")

Plotagem dos dados

ggplot(dados, aes(x = Data)) +
  geom_line(aes(y = Abertura, color = "Spot")) +
  geom_line(aes(y = Fechamento, color = "Futuro")) +
  labs(title = "CAFÉ ARÁBICA",
       x = "Data",
       y = "Preço",
       color = "Variáveis") +
  scale_color_manual(values = c("Spot" = "blue", "Futuro" = "red")) +
  theme_minimal()

Transformando em Séries Temporais e Deflacionando pelo pacote ‘deflateBR’

# Transformando em séries temporais
spot <- ts(dados$Abertura, start= c(2003,01), end=c(2024,03), frequency = 12)
futuro <- ts(dados$Fechamento, start= c(2003,01), end=c(2024,03), frequency = 12)

# Definindo o tempo
times <- seq(as.Date("2003/1/2"), by = "month", length.out = 255)

# Deflacionando
spotR <- deflate(spot, nominal_dates = times, real_date = "02/2024", index = "ipca")
futuroR <- deflate(futuro, nominal_dates = times, real_date = "02/2024", index = "ipca")

Plotando os Gráficos de Preço Spot e Futuro

par(mfrow = c(1,2))

plot(spotR, main = "CAFÉ ARÁBICA - Spot"); plot(futuroR, main = "CAFÉ ARÁBICA - Futuro")

Tranformando dem logaritmo

O objetivo é corrigir a heterocedásticidade

lspot <- log(spotR)
lfuturo <- log(futuroR)

plot(lspot, main = "CAFÉ ARÁBICA - log SPOT"); plot(lfuturo, main = "CAFÉ ARÁBICA - log FUTURO")

Fazendo a regressão em diferença

Para corrigir os problemas causados pela autocorrelação nós fazemos a regressão em diferença

lag <- stats::lag
Regressao <- dynlm(diff(lspot)~diff(lag(lfuturo,0)))
Regressao


Time series regression with "ts" data:
Start = 2003(2), End = 2024(3)

Call:
dynlm(formula = diff(lspot) ~ diff(lag(lfuturo, 0)))

Coefficients:
          (Intercept)  diff(lag(lfuturo, 0))  
            -0.003896               0.063751

Testes para resíduos

par(mfrow = c(1,1))
residuos <- residuals(Regressao); head(residuos)

   fev 2003    mar 2003    abr 2003    mai 2003    jun 2003    jul 2003 
-0.08058037  0.03593724 -0.02248757 -0.03612945  0.02368021  0.08939638

plot(residuos, type="l", col="red")
abline(h=0, col="blue", lw=3)

par(mfrow=c(1,2))
hist(residuos, main="", col="cadetblue", prob=T, xlab = names(residuos)[1], breaks = 30)
curve(expr=dnorm(x,mean=mean(residuos),sd=sd(residuos)),col="red",add= TRUE, lwd=2)
qqnorm(residuos, col="blue")
qqline(residuos, col="red")

Teste para a presença de normalidade

Teste de Jarque Bera:

jarqueberaTest(residuos)


Title:
 Jarque - Bera Normalality Test

Test Results:
  STATISTIC:
    X-squared: 2.4802
  P VALUE:
    Asymptotic p Value: 0.2894

Teste de Shapiro Wilk:

shapiro.test(residuos)


    Shapiro-Wilk normality test

data:  residuos
W = 0.99384, p-value = 0.3852

Testes para presença de heteroscedasticidade

Teste de Goldfeld-Quandt:

length(residuos); length(residuos)*0.15

[1] 254

[1] 38.1

gqtest(Regressao, fraction = 38.1, alternative = "greater")


    Goldfeld-Quandt test

data:  Regressao
GQ = 0.94757, df1 = 106, df2 = 105, p-value = 0.6088
alternative hypothesis: variance increases from segment 1 to 2

Teste de Breusch-Pagan-Godfrey:

  bptest(Regressao)


    studentized Breusch-Pagan test

data:  Regressao
BP = 0.55691, df = 1, p-value = 0.4555

Teste de White:

library(whitestrap)
white_test(Regressao)

White's test results

Null hypothesis: Homoskedasticity of the residuals
Alternative hypothesis: Heteroskedasticity of the residuals
Test Statistic: 1.73
P-value: 0.420027

Testes de autocorrelação

Teste de Durbin Watson:

dwtest(Regressao)


    Durbin-Watson test

data:  Regressao
DW = 2.1779, p-value = 0.9226
alternative hypothesis: true autocorrelation is greater than 0

Teste de Breusch-Godfrey (1978):

bgtest(Regressao, order = 4)


    Breusch-Godfrey test for serial correlation of order up to 4

data:  Regressao
LM test = 4.8329, df = 4, p-value = 0.3049

Teste de Arch:

ArchTest(residuos, lags = 4)


    ARCH LM-test; Null hypothesis: no ARCH effects

data:  residuos
Chi-squared = 4.6801, df = 4, p-value = 0.3217

Conclusão

Corrigindo todos os problemas de heterocedásticidade, autocorrelação e mantendo uma distribuição normal, podemos concluir que para cada unidade de futuro, há uma perda de 13,50% no preço spot e aumenta o risco em 0,98%.