σμ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.3·~55 minAvanzado

Esquemas de remuestreo y degeneración

Posición en la rutaMódulo 2 de 3 · Tema 3 de 5
2.3 Esquemas de remuestreo y degeneración
—Puntuación del lab— / 100XP total—Racha—

Cuando los pesos se colapsan

Los pesos se multiplican de un día a otro. Una partícula que explica mal unas pocas observaciones recibe un peso diminuto, y los pesos diminutos nunca se recuperan. Tras suficientes días una partícula lleva casi todo el peso y las demás son trabajo perdido. Es la degeneración de pesos, que mide el tamaño muestral efectivo 1/∑w~21/\sum\tilde w^{2}1/∑w~2. En SPY con 100010001000 partículas y sin remuestreo, vale 960960960 tras el primer día, 368368368 tras diez, 525252 tras cincuenta y unos 1,31{,}31,3 tras 250250250 días. Un filtro reducido a una partícula es una sola conjetura, no una nube.

ESS=1∑iw~i2,ESS<N2,u+iN,  i=0,…,N−1\mathrm{ESS}=\frac{1}{\sum_i\tilde w_i^{2}},\qquad \mathrm{ESS}<\frac N2,\qquad \frac{u+i}{N},\ \ i=0,\dots,N-1ESS=∑i​w~i2​1​,ESS<2N​,Nu+i​,  i=0,…,N−1(1)

El remuestreo cura la degeneración: se saca una nube nueva de NNN partículas en la que una partícula de peso w~i\tilde w_iw~i​ aparece en promedio Nw~iN\tilde w_iNw~i​ veces, y luego se da a cada copia el mismo peso 1/N1/N1/N. Hay varias formas de sacarla. El multinomial saca NNN uniformes independientes. El sistemático saca un solo uniforme uuu y usa NNN punteros equiespaciados (u+i)/N(u+i)/N(u+i)/N sobre los pesos acumulados. El estratificado saca un uniforme en cada una de las NNN rebanadas iguales. El residual da primero a cada partícula sus copias enteras, ⌊Nw~i⌋\lfloor N\tilde w_i\rfloor⌊Nw~i​⌋, y solo sortea el resto. El sistemático es el más barato y el más usado.

Remuestreo sistemático a manoTomemos cuatro partículas en −1,0,1,2-1,0,1,2−1,0,1,2 con pesos 0,058, 0,258, 0,426, 0,2580{,}058,\ 0{,}258,\ 0{,}426,\ 0{,}2580,058, 0,258, 0,426, 0,258. Los pesos acumulados son 0,058, 0,316, 0,742, 10{,}058,\ 0{,}316,\ 0{,}742,\ 10,058, 0,316, 0,742, 1. Saquemos u=0,3u=0{,}3u=0,3. Los punteros son (0,3+i)/4=0,075, 0,325, 0,575, 0,825(0{,}3+i)/4=0{,}075,\ 0{,}325,\ 0{,}575,\ 0{,}825(0,3+i)/4=0,075, 0,325, 0,575, 0,825. El puntero 0,0750{,}0750,075 cae en la segunda casilla, 0,3250{,}3250,325 y 0,5750{,}5750,575 en la tercera, 0,8250{,}8250,825 en la cuarta. La nube nueva es 0,1,1,20,1,1,20,1,1,2: la partícula en −1-1−1 desaparece y la pesada en 111 se copia dos veces, cerca de sus 4×0,426=1,74\times0{,}426=1{,}74×0,426=1,7 copias esperadas. Caso peligroso: con pesos 0,97, 0,01, 0,01, 0,010{,}97,\ 0{,}01,\ 0{,}01,\ 0{,}010,97, 0,01, 0,01, 0,01 el tamaño muestral efectivo es 1,061{,}061,06 y cualquier juego de punteros elige solo la primera partícula, así que la nube se reduce a un valor distinto.
  • Remuestrear cuando el ESS es bajo: regla común ESS<N/2\mathrm{ESS}<N/2ESS<N/2. Remuestrear cada día desperdicia diversidad.
  • El sistemático usa un solo número aleatorio: los punteros (u+i)/N(u+i)/N(u+i)/N sobre los pesos acumulados dan un sorteo barato y de poco ruido.
  • El precio es el empobrecimiento: muchas copias de pocas partículas dejan la nube con poca variedad, que el siguiente paso de propagación debe reconstruir.

Lab: qué se califica

  • Copias esperadas15 pts
  • Un puntero sistemático10 pts
  • Pesos colapsados15 pts
  • Supervivientes distintas10 pts
  • Hay que remuestrear hoy10 pts
  • Remuestreo sistemático10 pts
  • Emparejar los esquemas15 pts
  • Por qué remuestrear15 pts
0/ 8 respondidos
Ejercicios
1Respuesta numérica15 ptsSin responder

Con N=4N=4N=4 partículas, una partícula con peso normalizado 0,4260{,}4260,426 aparece Nw~N\tilde wNw~ veces tras el remuestreo. ¿Cuánto es ese número? Dé dos decimales.

n=n =n=aceptado: ± 0.02
2Respuesta numérica10 ptsSin responder

El remuestreo sistemático con N=4N=4N=4 y u=0,3u=0{,}3u=0,3 usa los punteros (u+i)/N(u+i)/N(u+i)/N. ¿Cuánto vale el segundo puntero, i=1i=1i=1? Dé tres decimales.

p=p =p=aceptado: ± 0.001
3Respuesta numérica15 ptsSin responder

Los pesos normalizados son 0,97, 0,01, 0,01, 0,010{,}97,\ 0{,}01,\ 0{,}01,\ 0{,}010,97, 0,01, 0,01, 0,01. ¿Cuánto vale el tamaño muestral efectivo 1/∑w~21/\sum\tilde w^{2}1/∑w~2? Dé dos decimales.

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

Los pesos son 0,97, 0,01, 0,01, 0,010{,}97,\ 0{,}01,\ 0{,}01,\ 0{,}010,97, 0,01, 0,01, 0,01 y el remuestreo sistemático usa N=4N=4N=4 punteros con u=0,5u=0{,}5u=0,5. ¿Cuántas partículas originales distintas aparecen en la nube nueva?

k=k =k=aceptado: ± 0.01
5Elección10 ptsSin responder

Hay N=1000N=1000N=1000 partículas y el tamaño muestral efectivo de hoy es 620620620. Con la regla «remuestrear si ESS < N/2», ¿qué hace?

6Ordenar los pasos10 ptsSin responder

Ponga en orden los pasos del remuestreo sistemático.

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

  1. 1Para cada puntero tomar la primera partícula cuyo peso acumulado lo alcanza
  2. 2Sumar los pesos en pesos acumulados
  3. 3Colocar N punteros equiespaciados (u + i) / N
  4. 4Sacar un solo número uniforme u
  5. 5Dar a cada copia el peso 1 / N
7Emparejar15 ptsSin responder

Empareje cada esquema de remuestreo con su forma de sortear.

Elige la pareja de cada elemento de la izquierda.

Multinomial
Sistemático
Estratificado
Residual
8Explicar con palabras15 ptsSin responder

En una o dos frases, explique por qué un filtro de partículas necesita remuestreo y qué cuesta.

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