σμ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 AcademyFeuilles de routeParcours
Ctrl K
AO
+1,2 k apprenants
ConnexionRejoindre ↗
Feuilles de routeParcours
AO
+1,2 k apprenants
Connexion
  • 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 AcademyFeuilles de routeParcours
Ctrl K
AO
+1,2 k apprenants
ConnexionRejoindre ↗
Feuilles de routeParcours
AO
+1,2 k apprenants
Connexion
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
Parcours Particle filtering and non-stationary time series›Module 2 · Filtres particulaires›Sujet 2.1·~55 minAvancé

L'échantillonnage préférentiel

Position dans le parcoursModule 2 sur 3 · Sujet 1 sur 5
2.1 L'échantillonnage préférentiel
—Score du lab— / 100XP total—Série—

Estimer avec des estimations pondérées

On veut souvent une moyenne sous une loi ppp dont on ne sait pas tirer d'échantillons. En filtrage, ppp est la croyance sur l'état caché après les données du jour. L'astuce de l'échantillonnage préférentiel est de tirer plutôt selon une loi facile ggg, appelée proposition, et de réparer l'écart avec des poids. Un point xix_ixi​ que la proposition produit trop souvent reçoit un petit poids, un point produit trop rarement reçoit un grand poids : wi=p(xi)/g(xi)w_i=p(x_i)/g(x_i)wi​=p(xi​)/g(xi​). Le plus souvent ppp n'est connue qu'à une constante près, donc on divise les poids par leur somme, ce qui les fait sommer à 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)

Dans un filtre, la proposition est la prévision (où l'état se trouve probablement avant les données) et la cible est la croyance après les données. Leur rapport est simplement la vraisemblance de l'observation : le poids d'une particule dit à quel point elle explique yyy. La taille d'échantillon effective, ESS, dit combien de tirages équipondérés vaut l'échantillon pondéré. Elle vaut NNN si tous les poids sont égaux et 111 si un seul point porte tout. Quand l'ESS est petite, l'estimation repose sur quelques points et devient bruitée. Cela arrive quand la proposition rate les zones où la cible a l'essentiel de sa masse, par exemple quand les données sont une grande surprise.

Quatre particules et une observationQuatre tirages de la prévision valent x=−1,0,1,2x=-1,0,1,2x=−1,0,1,2 et l'on observe y=1y=1y=1 avec une variance de bruit r=1r=1r=1. Les poids non normalisés sont 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, de somme 2,3482{,}3482,348. Les poids normalisés sont 0,058, 0,258, 0,426, 0,2580{,}058,\ 0{,}258,\ 0{,}426,\ 0{,}2580,058, 0,258, 0,426, 0,258. La moyenne pondérée est 0,8850{,}8850,885 et le moment d'ordre deux pondéré 1,5161{,}5161,516, donc la variance vaut 1,516−0,8852=0,7341{,}516-0{,}885^{2}=0{,}7341,516−0,8852=0,734. La taille d'échantillon effective est 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 sur 444.
  • Les poids réparent l'écart : tirer selon une ggg facile, pondérer par p/gp/gp/g, diviser par la somme.
  • Dans un filtre le poids est la vraisemblance : la particule qui explique le mieux l'observation reçoit le plus de poids.
  • Surveillez l'ESS : 1/∑w~21/\sum\tilde w^{2}1/∑w~2 chute quand quelques points dominent, signe d'une mauvaise proposition.

Lab : ce qui est noté

  • Un poids normalisé15 pts
  • La moyenne pondérée15 pts
  • La taille d'échantillon effective15 pts
  • La variance pondérée10 pts
  • Une proposition dangereuse10 pts
  • Les étapes de l'échantillonnage préférentiel10 pts
  • Associer le vocabulaire10 pts
  • Pourquoi les poids sont nécessaires15 pts
0/ 8 traités
Exercices
1Réponse numérique15 ptsSans réponse

Les poids non normalisés des quatre particules sont 0,135, 0,607, 1, 0,6070{,}135,\ 0{,}607,\ 1,\ 0{,}6070,135, 0,607, 1, 0,607, de somme 2,3482{,}3482,348. Quel est le poids normalisé de la première particule ? Donnez trois décimales.

w=w =w=accepté : ± 0.001
2Réponse numérique15 ptsSans réponse

Les particules x=−1,0,1,2x=-1,0,1,2x=−1,0,1,2 ont pour poids normalisés 0,058, 0,258, 0,426, 0,2580{,}058,\ 0{,}258,\ 0{,}426,\ 0{,}2580,058, 0,258, 0,426, 0,258. Que vaut la moyenne pondérée ∑w~ixi\sum\tilde w_i x_i∑w~i​xi​ ? Donnez trois décimales.

m=m =m=accepté : ± 0.003
3Réponse numérique15 ptsSans réponse

Avec les poids normalisés 0,058, 0,258, 0,426, 0,2580{,}058,\ 0{,}258,\ 0{,}426,\ 0{,}2580,058, 0,258, 0,426, 0,258, la somme des carrés des poids est 0,3180{,}3180,318. Que vaut la taille d'échantillon effective 1/∑w~i21/\sum\tilde w_i^{2}1/∑w~i2​ ? Donnez deux décimales.

ess=ess =ess=accepté : ± 0.02
4Réponse numérique10 ptsSans réponse

Le moment d'ordre deux pondéré ∑w~ixi2\sum\tilde w_i x_i^{2}∑w~i​xi2​ vaut 1,5161{,}5161,516 et la moyenne pondérée 0,8850{,}8850,885. Que vaut la variance pondérée ? Donnez trois décimales.

v=v =v=accepté : ± 0.004
5Choix10 ptsSans réponse

Quelle proposition risque le plus de donner une très petite taille d'échantillon effective ?

6Remettre les étapes dans l'ordre10 ptsSans réponse

Mettez dans l'ordre les étapes de l'échantillonnage préférentiel auto-normalisé.

Faites glisser pour réordonner, ou placez le focus sur un élément, Espace pour le soulever, les flèches pour le déplacer et Espace pour le déposer.

  1. 1Diviser chaque poids par la somme des poids
  2. 2Tirer des points selon la proposition facile
  3. 3Moyenner la quantité avec les poids normalisés
  4. 4Calculer un poids pour chaque point : cible sur proposition
7Associer10 ptsSans réponse

Associez chaque terme à son sens.

Choisissez le partenaire de chaque élément de gauche.

Proposition
Cible
Poids d'importance
Taille d'échantillon effective
8Expliquer avec des mots15 ptsSans réponse

En une ou deux phrases, expliquez pourquoi il faut pondérer les tirages par p/gp/gp/g quand on échantillonne la proposition ggg au lieu de la cible ppp.

Mode mathsSaisie au clavier ou en LaTeXExercice 1 sur 8 · Ctrl ⏎ pour vérifier

Connectez-vous pour exécuter votre code.

Tous les cours sont gratuits à lire. Créez un compte ou connectez-vous pour exécuter du code, lancer les tests et soumettre ce lab.

  • IDE dans le navigateur avec terminal bash — rien à installer
  • Correction instantanée, XP et votre place dans le parcours
Créer un compteSe connecter
OU
Continuer avec GitHub