σμ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 Bayesian Inference›Módulo 1 · Crenças e evidências›Tópico 1.4·~1 hIntermediário

Regressão linear bayesiana

Posição no roteiroMódulo 1 de 3 · Tópico 4 de 4
1.4 Regressão linear bayesiana
—Nota do lab— / 100XP total—Sequência—

Regressão com incerteza honesta

A regressão comum dá uma única inclinação. A regressão bayesiana dá toda uma distribuição para a inclinação. Tome o retorno yyy de uma ação e o do mercado xxx com y=βx+ey=\beta x+ey=βx+e, com ruído eee de variância σ2\sigma^{2}σ2. Parta de um prior para o beta, por exemplo uma normal centrada em m0=1m_0=1m0​=1 com variância v0v_0v0​. Os dados empurram a inclinação para o valor de mínimos quadrados e o prior a puxa de volta. É a mesma média ponderada pela precisão dos tópicos anteriores. A precisão do prior é 1/v01/v_01/v0​, a dos dados é ∑x2/σ2\sum x^{2}/\sigma^{2}∑x2/σ2, e esta última é grande quando o mercado se mexeu muito.

Pn=1v0+∑x2σ2,mn=m0/v0+∑xy/σ2Pn,sd(β)=1PnP_n=\frac{1}{v_0}+\frac{\sum x^{2}}{\sigma^{2}},\qquad m_n=\frac{m_0/v_0+\sum xy/\sigma^{2}}{P_n},\qquad \mathrm{sd}(\beta)=\frac{1}{\sqrt{P_n}}Pn​=v0​1​+σ2∑x2​,mn​=Pn​m0​/v0​+∑xy/σ2​,sd(β)=Pn​​1​(1)

Há duas coisas a notar. Primeiro, um prior centrado em 000 dá exatamente a regressão ridge: a inclinação é ∑xy / (∑x2+λ)\sum xy\,/\,(\sum x^{2}+\lambda)∑xy/(∑x2+λ) com λ=σ2/v0\lambda=\sigma^{2}/v_0λ=σ2/v0​. Então a regressão ridge é uma regressão bayesiana cujo prior desconfia de inclinações grandes. Segundo, uma previsão para um novo x∗x^{*}x∗ tem duas fontes de dispersão: o ruído σ2\sigma^{2}σ2 e a incerteza da inclinação, que soma x∗2/Pnx^{*2}/P_nx∗2/Pn​. A segunda cresce com o quadrado de x∗x^{*}x∗, então previsões longe dos dados são menos seguras.

Um beta em númerosSeja ∑x2=40\sum x^{2}=40∑x2=40, ∑xy=48\sum xy=48∑xy=48 e σ2=4\sigma^{2}=4σ2=4. A inclinação de mínimos quadrados é 48/40=1,248/40=1{,}248/40=1,2. Com o prior N(1, 0,52)N(1,\,0{,}5^{2})N(1,0,52), ou seja v0=0,25v_0=0{,}25v0​=0,25: a precisão é 4+40/4=144+40/4=144+40/4=14, a média (4+12)/14=1,143(4+12)/14=1{,}143(4+12)/14=1,143 e o desvio padrão 1/14=0,2671/\sqrt{14}=0{,}2671/14​=0,267. Com um prior ridge centrado em 000 a inclinação seria 48/(40+16)=0,85748/(40+16)=0{,}85748/(40+16)=0,857. Uma previsão em x∗=2x^{*}=2x∗=2 tem média 2,292{,}292,29 e desvio padrão 4+4/14=2,07\sqrt{4+4/14}=2{,}074+4/14​=2,07.
  • Inclinação bayesiana: uma distribuição cuja precisão é a do prior mais a dos dados.
  • A regressão ridge é o caso de um prior centrado em 000, com λ=σ2/v0\lambda=\sigma^{2}/v_0λ=σ2/v0​.
  • As previsões carregam ruído mais incerteza da inclinação, e a segunda cresce com x∗2x^{*2}x∗2.

Lab: o que é avaliado

  • A precisão posterior15 pts
  • A inclinação posterior15 pts
  • A incerteza da inclinação10 pts
  • A inclinação ridge15 pts
  • A dispersão de uma previsão10 pts
  • Um prior muito confiante10 pts
  • Associar o vocabulário10 pts
  • Longe dos dados15 pts
0/ 8 respondidas
Tarefas
1Resposta numérica15 ptsSem resposta

Com ∑x2=40\sum x^{2}=40∑x2=40, σ2=4\sigma^{2}=4σ2=4 e variância prévia v0=0,25v_0=0{,}25v0​=0,25, qual é a precisão posterior 1/v0+∑x2/σ21/v_0+\sum x^{2}/\sigma^{2}1/v0​+∑x2/σ2?

p=p =p=aceito: ± 0.01
2Resposta numérica15 ptsSem resposta

Agora some a média prévia m0=1m_0=1m0​=1 e ∑xy=48\sum xy=48∑xy=48. Qual é a média posterior (m0/v0+∑xy/σ2)/14\big(m_0/v_0+\sum xy/\sigma^{2}\big)/14(m0​/v0​+∑xy/σ2)/14? Dê quatro casas decimais.

m=m =m=aceito: ± 0.001
3Resposta numérica10 ptsSem resposta

Qual é o desvio padrão posterior da inclinação, 1/141/\sqrt{14}1/14​? Dê quatro casas decimais.

s=s =s=aceito: ± 0.001
4Resposta numérica15 ptsSem resposta

Use um prior centrado em 000 com v0=0,25v_0=0{,}25v0​=0,25 e σ2=4\sigma^{2}=4σ2=4, ou seja λ=σ2/v0=16\lambda=\sigma^{2}/v_0=16λ=σ2/v0​=16. Qual é a inclinação ridge ∑xy / (∑x2+λ)\sum xy\,/\,(\sum x^{2}+\lambda)∑xy/(∑x2+λ) com ∑xy=48\sum xy=48∑xy=48 e ∑x2=40\sum x^{2}=40∑x2=40? Dê quatro casas decimais.

r=r =r=aceito: ± 0.001
5Resposta numérica10 ptsSem resposta

Para um novo x∗=2x^{*}=2x∗=2 o desvio padrão preditivo é σ2+x∗2/14\sqrt{\sigma^{2}+x^{*2}/14}σ2+x∗2/14​ com σ2=4\sigma^{2}=4σ2=4. Quanto vale? Dê três casas decimais.

p=p =p=aceito: ± 0.01
6Escolha10 ptsSem resposta

O que acontece com a inclinação posterior se a variância prévia v0v_0v0​ ficar muito pequena?

7Associar10 ptsSem resposta

Associe cada termo ao seu significado.

Escolha o par de cada item da esquerda.

Inclinação de mínimos quadrados
Regressão ridge
Precisão posterior
Variância preditiva
8Explicar com palavras15 ptsSem resposta

Em uma ou duas frases, explique por que uma regressão bayesiana dá faixas de previsão mais largas para valores de xxx longe dos dados.

Modo matemáticaEntrada por teclado ou LaTeXTarefa 1 de 8 · 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