σμimplied vol · KKmax(S−K,0)udβρᵢⱼfat tailsPnL · drawdownyield · TdSₜ = μSₜ dt + σSₜ dWₜC = S·N(d₁) − K·e⁻ʳᵀ·N(d₂)σₚ² = wᵀΣwd₁ = [ln(S/K) + (r + σ²⁄2)T] ⁄ σ√T∂V/∂t + ½σ²S²·∂²V/∂S² + rS·∂V/∂S − rV = 0Sharpe = (Rₚ − R_f) ⁄ σₚVaR₉₉ = μ − 2.33σρ = Cov(X,Y) ⁄ σₓσᵧE[Rₚ] = Σ wᵢ·E[Rᵢ]φ(x) = e^(−x²/2) ⁄ √2πE[Rᵢ] = R_f + βᵢ(E[Rₘ] − R_f)dXₜ = θ(μ − Xₜ)dt + σ dWₜσ²ₜ = ω + α·ε²ₜ₋₁ + β·σ²ₜ₋₁C − P = S − K·e⁻ʳᵀΔ = ∂V/∂S Γ = ∂²V/∂S²p = (e^(rΔt) − d) ⁄ (u − d)rₜ = ln(Pₜ ⁄ Pₜ₋₁)E[Xₜ₊₁ | ℱₜ] = Xₜf* = (bp − q) ⁄ bz = (x − μ) ⁄ σP(N=k) = λᵏe^(−λ) ⁄ k!Θ = −∂V/∂Tν = S·φ(d₁)·√Tmin wᵀΣw s.t. wᵀμ = μ*∫Σ∂∞λ
uantCraft AcademyRoadmapsRoteiros
Ctrl K
AO
+1,2 mil alunos
EntrarParticipar ↗
RoadmapsRoteiros
AO
+1,2 mil alunos
Entrar
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
σμimplied vol · KKmax(S−K,0)udβρᵢⱼfat tailsPnL · drawdownyield · TdSₜ = μSₜ dt + σSₜ dWₜC = S·N(d₁) − K·e⁻ʳᵀ·N(d₂)σₚ² = wᵀΣwd₁ = [ln(S/K) + (r + σ²⁄2)T] ⁄ σ√T∂V/∂t + ½σ²S²·∂²V/∂S² + rS·∂V/∂S − rV = 0Sharpe = (Rₚ − R_f) ⁄ σₚVaR₉₉ = μ − 2.33σρ = Cov(X,Y) ⁄ σₓσᵧE[Rₚ] = Σ wᵢ·E[Rᵢ]φ(x) = e^(−x²/2) ⁄ √2πE[Rᵢ] = R_f + βᵢ(E[Rₘ] − R_f)dXₜ = θ(μ − Xₜ)dt + σ dWₜσ²ₜ = ω + α·ε²ₜ₋₁ + β·σ²ₜ₋₁C − P = S − K·e⁻ʳᵀΔ = ∂V/∂S Γ = ∂²V/∂S²p = (e^(rΔt) − d) ⁄ (u − d)rₜ = ln(Pₜ ⁄ Pₜ₋₁)E[Xₜ₊₁ | ℱₜ] = Xₜf* = (bp − q) ⁄ bz = (x − μ) ⁄ σP(N=k) = λᵏe^(−λ) ⁄ k!Θ = −∂V/∂Tν = S·φ(d₁)·√Tmin wᵀΣw s.t. wᵀμ = μ*∫Σ∂∞λ
uantCraft AcademyRoadmapsRoteiros
Ctrl K
AO
+1,2 mil alunos
EntrarParticipar ↗
RoadmapsRoteiros
AO
+1,2 mil alunos
Entrar
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
Roteiro de Statistics›Módulo 3 · Modelos e regressão›Tópico 3.4·~50 minIntermediário

Multicolinearidade e regularização

Posição no roteiroMódulo 3 de 4 · Tópico 4 de 4
3.4 Multicolinearidade e regularização
—Nota do lab— / 100XP total—Sequência—

Preditores que dizem a mesma coisa

Quando dois regressores se movem quase juntos, como o momentum de uma ação em 6 e em 12 meses, a regressão não consegue dizer qual merece o crédito. É a multicolinearidade. O ajuste pode continuar bom, mas os coeficientes individuais ficam instáveis: uma pequena mudança nos dados pode subir um e descer o outro, e seus erros-padrão explodem. O fator de inflação da variância do regressor jjj é VIFj=11−Rj2\mathrm{VIF}_j=\dfrac{1}{1-R_j^2}VIFj​=1−Rj2​1​, em que Rj2R_j^2Rj2​ vem da regressão desse regressor sobre todos os outros. Um VIF acima de 5 ou 10 é um sinal de alerta.

β^ridge=(XTX+λI)−1XTy,β^lasso=arg⁡min⁡ ∥y−Xβ∥2+λ∥β∥1\hat{\boldsymbol\beta}_{\text{ridge}}=(X^{T}X+\lambda I)^{-1}X^{T}\mathbf{y},\qquad \hat{\boldsymbol\beta}_{\text{lasso}}=\arg\min\ \|\mathbf{y}-X\boldsymbol\beta\|^{2}+\lambda\|\boldsymbol\beta\|_{1}β^​ridge​=(XTX+λI)−1XTy,β^​lasso​=argmin ∥y−Xβ∥2+λ∥β∥1​(1)

A regularização corrige a instabilidade acrescentando uma penalização ao tamanho dos coeficientes. O ridge (L2) soma λ∥β∥2\lambda\|\boldsymbol\beta\|^2λ∥β∥2: a solução é a ordinária com λI\lambda IλI somado a XTXX^{T}XXTX, o que condiciona bem a matriz (relacione com a habilidade de álgebra linear) e encolhe cada coeficiente para zero. O lasso (L1) soma λ∥β∥1\lambda\|\boldsymbol\beta\|_1λ∥β∥1​ e pode zerar exatamente alguns coeficientes, selecionando preditores. Ambos trocam um pouco de viés por uma grande queda de variância, e λ\lambdaλ é escolhido por validação cruzada, o próximo tópico.

Números para preditores correlacionadosSe regredir o preditor jjj sobre os outros dá Rj2=0,9R_j^2=0{,}9Rj2​=0,9, então VIFj=1/(1−0,9)=10\mathrm{VIF}_j=1/(1-0{,}9)=10VIFj​=1/(1−0,9)=10: seu erro-padrão é cerca de 10≈3,2\sqrt{10}\approx3{,}210​≈3,2 vezes maior do que com preditores não correlacionados. Inversamente, VIF=5\mathrm{VIF}=5VIF=5 significa Rj2=0,8R_j^2=0{,}8Rj2​=0,8. Com XTX=(4334)X^{T}X=\begin{pmatrix}4&3\\3&4\end{pmatrix}XTX=(43​34​), o ridge com λ=1\lambda=1λ=1 usa (5335)\begin{pmatrix}5&3\\3&5\end{pmatrix}(53​35​). Em uma dimensão com XTX=4X^{T}X=4XTX=4 e XTy=8X^{T}y=8XTy=8, o MQO dá β^=2\hat\beta=2β^​=2 enquanto o ridge com λ=4\lambda=4λ=4 dá 8/(4+4)=18/(4+4)=18/(4+4)=1: a penalização reduz o coeficiente pela metade.
  • Multicolinearidade: coeficientes instáveis e erros-padrão inflados, previsões muitas vezes ainda boas.
  • Ridge: encolhe todos os coeficientes, nunca exatamente a zero. Lasso: pode zerar coeficientes exatamente.
  • Escolha a penalização λ\lambdaλ por validação cruzada, nunca olhando o ajuste de treino.
Não leia um coeficiente quando seu gêmeo está no modeloCom dois fatores quase duplicados, cada coeficiente sozinho não é interpretável: um positivo grande e um negativo grande podem se cancelar e dar uma soma razoável. Julgue o par em conjunto, elimine ou combine um deles, ou regularize. As previsões do modelo podem continuar boas, porque a combinação é bem determinada mesmo quando cada parte não é.

Lab: o que é avaliado

  • Um fator de inflação da variância15 pts
  • Do VIF ao R quadrado10 pts
  • A matriz ridge15 pts
  • Um coeficiente encolhido15 pts
  • Qual penalização seleciona variáveis15 pts
  • O que a colinearidade faz15 pts
  • Por que o ridge ajuda15 pts
0/ 7 respondidas
Tarefas
1Resposta numérica15 ptsSem resposta

Regredir um preditor sobre todos os outros dá R quadrado = 0,9. Qual é o seu VIF?

VIF=VIF =VIF=aceito: ± 0.001
2Resposta numérica10 ptsSem resposta

Um preditor tem VIF = 5. Qual é o R quadrado da sua regressão sobre os outros preditores?

R2=R^2 =R2=aceito: ± 0.001
3Matriz · vetor15 ptsSem resposta

Seja XTX=(4334)X^{T}X=\begin{pmatrix}4&3\\3&4\end{pmatrix}XTX=(43​34​). Preencha XTX+λIX^{T}X+\lambda IXTX+λI para λ=1\lambda=1λ=1.

XTX+I=X^T X + I =XTX+I=
4Resposta numérica15 ptsSem resposta

Com um único regressor, X^T X = 4 e X^T y = 8. O MQO dá 2. O que o ridge dá com lambda = 4, usando X^T y / (X^T X + lambda)?

beta=beta =beta=aceito: ± 0.001
5Escolha15 ptsSem resposta

Qual penalização pode zerar alguns coeficientes exatamente?

6Escolha15 ptsSem resposta

Selecione todos os efeitos de uma forte multicolinearidade.

7Explicar com palavras15 ptsSem resposta

Em duas frases, diga por que a regressão ridge ajuda quando os preditores são fortemente correlacionados.

Modo matemáticaEntrada por teclado ou LaTeXTarefa 1 de 7 · Ctrl ⏎ para verificar

Entre para executar o seu código.

Todas as aulas são gratuitas para ler. Crie uma conta ou entre para executar código, rodar os testes e enviar este lab.

  • IDE no navegador com terminal bash — nada para instalar
  • Correção instantânea, XP e o seu lugar no roteiro
Criar uma contaEntrar
OU
Continuar com o GitHub