σμ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 AcademyHojas de rutaRutas
Ctrl K
AO
+1,2 k estudiantes
Iniciar sesiónUnirse ↗
Hojas de rutaRutas
AO
+1,2 k estudiantes
Iniciar sesión
  • 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 AcademyHojas de rutaRutas
Ctrl K
AO
+1,2 k estudiantes
Iniciar sesiónUnirse ↗
Hojas de rutaRutas
AO
+1,2 k estudiantes
Iniciar sesión
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
Ruta de Particle filtering and non-stationary time series›Módulo 2 · Filtros de partículas›Tema 2.1·~55 minAvanzado

Muestreo por importancia

Posición en la rutaMódulo 2 de 3 · Tema 1 de 5
2.1 Muestreo por importancia
—Puntuación del lab— / 100XP total—Racha—

Estimar con conjeturas ponderadas

A menudo se quiere un promedio bajo una distribución ppp de la que no se puede muestrear. En filtrado, ppp es la creencia sobre el estado oculto tras ver los datos de hoy. El truco del muestreo por importancia es muestrear de una distribución fácil ggg, llamada propuesta, y reparar la discrepancia con pesos. Un punto xix_ixi​ que la propuesta produce demasiado a menudo recibe un peso pequeño, y uno que produce muy rara vez recibe un peso grande: wi=p(xi)/g(xi)w_i=p(x_i)/g(x_i)wi​=p(xi​)/g(xi​). Normalmente ppp solo se conoce salvo una constante, así que los pesos se dividen por su suma, con lo que suman 111.

wi=p(xi)g(xi),w~i=wi∑jwj,Ep[f]≈∑iw~i f(xi),ESS=1∑iw~i2w_i=\frac{p(x_i)}{g(x_i)},\quad \tilde w_i=\frac{w_i}{\sum_j w_j},\quad E_p[f]\approx\sum_i\tilde w_i\,f(x_i),\quad \mathrm{ESS}=\frac{1}{\sum_i\tilde w_i^{2}}wi​=g(xi​)p(xi​)​,w~i​=∑j​wj​wi​​,Ep​[f]≈i∑​w~i​f(xi​),ESS=∑i​w~i2​1​(1)

En un filtro, la propuesta es la predicción (dónde es probable que esté el estado antes de los datos) y el objetivo es la creencia tras los datos. Su cociente es simplemente la verosimilitud de la observación: el peso de una partícula dice cuán bien explica yyy. El tamaño muestral efectivo, ESS, dice a cuántas muestras equiponderadas equivale la muestra ponderada. Vale NNN si todos los pesos son iguales y 111 si un solo punto lo carga todo. Cuando el ESS es pequeño, la estimación descansa en pocos puntos y es ruidosa. Pasa cuando la propuesta no cubre donde el objetivo tiene su masa, por ejemplo cuando los datos son una gran sorpresa.

Cuatro partículas y una observaciónCuatro muestras de la predicción son x=−1,0,1,2x=-1,0,1,2x=−1,0,1,2 y se observa y=1y=1y=1 con varianza de ruido r=1r=1r=1. Los pesos sin normalizar son e−(1−x)2/2=0,135, 0,607, 1, 0,607e^{-(1-x)^2/2}=0{,}135,\ 0{,}607,\ 1,\ 0{,}607e−(1−x)2/2=0,135, 0,607, 1, 0,607, que suman 2,3482{,}3482,348. Los pesos normalizados son 0,058, 0,258, 0,426, 0,2580{,}058,\ 0{,}258,\ 0{,}426,\ 0{,}2580,058, 0,258, 0,426, 0,258. La media ponderada es 0,8850{,}8850,885 y el segundo momento ponderado 1,5161{,}5161,516, así que la varianza es 1,516−0,8852=0,7341{,}516-0{,}885^{2}=0{,}7341,516−0,8852=0,734. El tamaño muestral efectivo es 1/(0,0582+0,2582+0,4262+0,2582)=3,141/(0{,}058^2+0{,}258^2+0{,}426^2+0{,}258^2)=3{,}141/(0,0582+0,2582+0,4262+0,2582)=3,14 de 444.
  • Los pesos reparan la discrepancia: muestrear de una ggg fácil, ponderar por p/gp/gp/g, dividir por la suma.
  • En un filtro el peso es la verosimilitud: la partícula que mejor explica la observación recibe más peso.
  • Vigile el ESS: 1/∑w~21/\sum\tilde w^{2}1/∑w~2 cae cuando unos pocos puntos dominan, señal de una mala propuesta.

Lab: qué se califica

  • Un peso normalizado15 pts
  • La media ponderada15 pts
  • El tamaño muestral efectivo15 pts
  • La varianza ponderada10 pts
  • Una propuesta peligrosa10 pts
  • Los pasos del muestreo por importancia10 pts
  • Emparejar el vocabulario10 pts
  • Por qué hacen falta los pesos15 pts
0/ 8 respondidos
Ejercicios
1Respuesta numérica15 ptsSin responder

Los pesos sin normalizar de las cuatro partículas son 0,135, 0,607, 1, 0,6070{,}135,\ 0{,}607,\ 1,\ 0{,}6070,135, 0,607, 1, 0,607 y suman 2,3482{,}3482,348. ¿Cuál es el peso normalizado de la primera partícula? Dé tres decimales.

w=w =w=aceptado: ± 0.001
2Respuesta numérica15 ptsSin responder

Las partículas x=−1,0,1,2x=-1,0,1,2x=−1,0,1,2 tienen pesos normalizados 0,058, 0,258, 0,426, 0,2580{,}058,\ 0{,}258,\ 0{,}426,\ 0{,}2580,058, 0,258, 0,426, 0,258. ¿Cuánto vale la media ponderada ∑w~ixi\sum\tilde w_i x_i∑w~i​xi​? Dé tres decimales.

m=m =m=aceptado: ± 0.003
3Respuesta numérica15 ptsSin responder

Con los pesos normalizados 0,058, 0,258, 0,426, 0,2580{,}058,\ 0{,}258,\ 0{,}426,\ 0{,}2580,058, 0,258, 0,426, 0,258, la suma de los cuadrados de los pesos es 0,3180{,}3180,318. ¿Cuánto vale el tamaño muestral efectivo 1/∑w~i21/\sum\tilde w_i^{2}1/∑w~i2​? Dé dos decimales.

ess=ess =ess=aceptado: ± 0.02
4Respuesta numérica10 ptsSin responder

El segundo momento ponderado ∑w~ixi2\sum\tilde w_i x_i^{2}∑w~i​xi2​ es 1,5161{,}5161,516 y la media ponderada 0,8850{,}8850,885. ¿Cuánto vale la varianza ponderada? Dé tres decimales.

v=v =v=aceptado: ± 0.004
5Elección10 ptsSin responder

¿Qué propuesta es más probable que dé un tamaño muestral efectivo diminuto?

6Ordenar los pasos10 ptsSin responder

Ponga en orden los pasos del muestreo por importancia autonormalizado.

Arrastra para reordenar, o enfoca un elemento, pulsa Espacio para levantarlo, las flechas para moverlo y Espacio para soltarlo.

  1. 1Dividir cada peso por la suma de los pesos
  2. 2Muestrear puntos de la propuesta fácil
  3. 3Promediar la cantidad con los pesos normalizados
  4. 4Calcular un peso para cada punto: objetivo entre propuesta
7Emparejar10 ptsSin responder

Empareje cada término con su significado.

Elige la pareja de cada elemento de la izquierda.

Propuesta
Objetivo
Peso de importancia
Tamaño muestral efectivo
8Explicar con palabras15 ptsSin responder

En una o dos frases, explique por qué hay que ponderar las muestras por p/gp/gp/g cuando se muestrea la propuesta ggg en lugar del objetivo ppp.

Modo matemáticasEntrada por teclado o LaTeXEjercicio 1 de 8 · Ctrl ⏎ para comprobar

Inicia sesión para ejecutar tu código.

Todas las clases se pueden leer gratis. Crea una cuenta o inicia sesión para ejecutar código, lanzar las pruebas y enviar este lab.

  • IDE en el navegador con terminal bash — sin instalar nada
  • Corrección instantánea, XP y tu lugar en la ruta
Crear una cuentaIniciar sesión
O
Continuar con GitHub