Loi de von Mises

Loi de von Mises
Image illustrative de l’article Loi de von Mises
Densité de probabilité
Densité de la loi de von Mises sur [–π, π]

Image illustrative de l’article Loi de von Mises
Fonction de répartition
Répartition de la loi de von Mises sur [–π, π]

Paramètres μ {\displaystyle \mu } réel,
κ > 0 {\displaystyle \kappa >0}
Support x {\displaystyle x\in } tout intervalle de longueur
Densité de probabilité e κ cos ( x μ ) 2 π I 0 ( κ ) {\displaystyle {\frac {{\rm {e}}^{\kappa \cos(x-\mu )}}{2\pi I_{0}(\kappa )}}}
Fonction de répartition ne peut être exprimée
Espérance μ {\displaystyle \mu }
Médiane μ {\displaystyle \mu }
Mode μ {\displaystyle \mu }
Variance var ( x ) = 1 I 1 ( κ ) / I 0 ( κ ) {\displaystyle {\textrm {var}}(x)=1-I_{1}(\kappa )/I_{0}(\kappa )} (circulaire)
Entropie κ I 1 ( κ ) I 0 ( κ ) + ln [ 2 π I 0 ( κ ) ] {\displaystyle -\kappa {\frac {I_{1}(\kappa )}{I_{0}(\kappa )}}+\ln[2\pi I_{0}(\kappa )]} (différentielle)
Fonction caractéristique I | n | ( κ ) I 0 ( κ ) e i n μ {\displaystyle {\frac {I_{|n|}(\kappa )}{I_{0}(\kappa )}}{\rm {e}}^{{\rm {i}}n\mu }}
modifier Consultez la documentation du modèle

Dans la théorie des probabilités et en statistiques, la loi (distribution) de von Mises (appelée également distribution normale circulaire ou distribution de Tikhonov) est une densité de probabilité continue, nommée d'après Richard von Mises. Elle donne une bonne approximation de la loi normale périodique (en), qui est l'analogue circulaire de la loi normale. Un angle de diffusion θ {\displaystyle \theta } parcourant un cercle est une variable aléatoire suivant la loi normale périodique avec une variance non périodique qui croît linéairement en temps. D'un autre côté, la loi de von Mises est la distribution stationnaire d'un processus de diffusion et déviation sur le cercle dans un potentiel harmonique, i.e. avec une orientation guidée[1].

La loi de von Mises est la loi de probabilités à entropie maximale pour une valeur donnée de z = e i θ {\displaystyle z={\rm {e}}^{{\rm {i}}\theta }} . La loi de von Mises est un cas particulier de la loi de von Mises-Fisher sur la N-sphère.

Définition

La densité de probabilités de la loi de von Mises pour un angle x est donnée par[2] :

f ( x μ , κ ) = e κ cos ( x μ ) 2 π I 0 ( κ ) {\displaystyle f(x\mid \mu ,\kappa )={\frac {{\rm {e}}^{\kappa \cos(x-\mu )}}{2\pi I_{0}(\kappa )}}}

I0(x) est la fonction de Bessel modifiée d'ordre 0.

Les paramètres μ et 1/κ sont analogues aux μ et σ2 (moyenne et variance) de la loi normale :

  • μ est une mesure de localisation (la distribution est centrée autour de μ),
  • κ est une mesure de la concentration (une mesure réciproque de la dispersion statistique).
    • Si κ est nul, la distribution est uniforme.
    • Si κ est très grand, la distribution devient très concentrée autour de μ.

La densité de probabilités peut être exprimée comme une série de fonctions de Bessel[3]

f ( x μ , κ ) = 1 2 π ( 1 + 2 I 0 ( κ ) j = 1 I j ( κ ) cos [ j ( x μ ) ] ) {\displaystyle f(x\mid \mu ,\kappa )={\frac {1}{2\pi }}\left(1+{\frac {2}{I_{0}(\kappa )}}\sum _{j=1}^{\infty }I_{j}(\kappa )\cos[j(x-\mu )]\right)}

Ij(x) est la fonction de Bessel modifiée d'ordre j.

La fonction de répartition n'est pas analytique et est généralement calculée comme intégrale de la série donnée précédemment. L'intégrale indéfinie de la densité de probabilités est :

Φ ( x μ , κ ) = f ( t μ , κ ) d t = 1 2 π ( x + 2 I 0 ( κ ) j = 1 I j ( κ ) sin [ j ( x μ ) ] j ) . {\displaystyle \Phi (x\mid \mu ,\kappa )=\int f(t\mid \mu ,\kappa )\,{\rm {d}}t={\frac {1}{2\pi }}\left(x+{\frac {2}{I_{0}(\kappa )}}\sum _{j=1}^{\infty }I_{j}(\kappa ){\frac {\sin[j(x-\mu )]}{j}}\right).}

Moments

Les moments de la loi de von Mises sont habituellement calculés comme les moments de z = eix plutôt que de l'angle x. On appelle ces moments « moments circulaires ». La variance calculée à partir de ces moments est également appelée « variance circulaire ». Toutefois, on désigne par la « moyenne » l'argument de la moyenne circulaire, et non la moyenne circulaire elle-même.

Le moment d'ordre n de z est :

m n = z n = Γ z n f ( x | μ , κ ) d x = I | n | ( κ ) I 0 ( κ ) e i n μ {\displaystyle m_{n}=\langle z^{n}\rangle =\int _{\Gamma }z^{n}\,f(x|\mu ,\kappa )\,{\rm {d}}x={\frac {I_{|n|}(\kappa )}{I_{0}(\kappa )}}{\rm {e}}^{{\rm {i}}n\mu }}

où l'intégrale se fait sur tout intervalle Γ de longueur . En calculant cette intégrale, on utilise le fait que zn = cos(nx) + i sin(nx) et l'identité de Bessel[4] :

I n ( κ ) = 1 π 0 π e κ cos ( x ) cos ( n x ) d x . {\displaystyle I_{n}(\kappa )={\frac {1}{\pi }}\int _{0}^{\pi }{\rm {e}}^{\kappa \cos(x)}\cos(nx)\,{\rm {d}}x.}

La moyenne de z  est alors donnée par

m 1 = I 1 ( κ ) I 0 ( κ ) e i μ {\displaystyle m_{1}={\frac {I_{1}(\kappa )}{I_{0}(\kappa )}}{\rm {e}}^{{\rm {i}}\mu }}

est la valeur « moyenne » de x est ainsi l'argument μ. Il s'agit donc de la direction moyenne de variables aléatoires angulaires. La variance de z, ou variance circulaire de x est :

Var ( x ) = 1 E [ cos ( x μ ) ] = 1 I 1 ( κ ) I 0 ( κ ) . {\displaystyle {\textrm {Var}}(x)=1-E[\cos(x-\mu )]=1-{\frac {I_{1}(\kappa )}{I_{0}(\kappa )}}.}

Entropie

L'entropie de la loi de von Mises est donnée par[2] :

H = Γ f ( θ ; μ , κ ) ln ( f ( θ ; μ , κ ) ) d θ {\displaystyle H=-\int _{\Gamma }f(\theta ;\mu ,\kappa )\,\ln(f(\theta ;\mu ,\kappa ))\,{\rm {d}}\theta \,}

Γ est un intervalle de longueur . Le logarithme de la densité de la loi de von Mises est directement donné par :

ln ( f ( θ ; μ , κ ) ) = ln ( 2 π I 0 ( κ ) ) + κ cos ( θ ) . {\displaystyle \ln(f(\theta ;\mu ,\kappa ))=-\ln(2\pi I_{0}(\kappa ))+\kappa \cos(\theta ).}

La fonction caractéristique de la loi est donc :

f ( θ ; μ , κ ) = 1 2 π ( 1 + 2 n = 1 ϕ n cos ( n θ ) ) {\displaystyle f(\theta ;\mu ,\kappa )={\frac {1}{2\pi }}\left(1+2\sum _{n=1}^{\infty }\phi _{n}\cos(n\theta )\right)}

avec ϕ n = I | n | ( κ ) / I 0 ( κ ) {\displaystyle \phi _{n}=I_{|n|}(\kappa )/I_{0}(\kappa )} . Les calculs donnent alors :

H = ln ( 2 π I 0 ( κ ) ) κ ϕ 1 = ln ( 2 π I 0 ( κ ) ) κ I 1 ( κ ) I 0 ( κ ) {\displaystyle H=\ln(2\pi I_{0}(\kappa ))-\kappa \phi _{1}=\ln(2\pi I_{0}(\kappa ))-\kappa {\frac {I_{1}(\kappa )}{I_{0}(\kappa )}}}

On retrouve, pour κ = 0, que la loi de von Mises est la loi uniforme circulaire (en) et l'entropie devient maximale et vaut ln(2π).

Voir aussi

Références

  • (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Von Mises distribution » (voir la liste des auteurs).
  1. (en) H. Risken, The Fokker–Planck Equation, Springer, , 472 p. (ISBN 978-3-540-61530-9, lire en ligne).
  2. a et b (en) Kantilal Mardia (en) et Peter E. Jupp, Directional Statistics, Wiley, (ISBN 978-0-471-95333-3, lire en ligne).
  3. Voir (en) Milton Abramowitz et Irene Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables [détail de l’édition] (lire en ligne) §9.6.34.
  4. Voir Abramowitz Stegun, §9.6.19.
  • (en) Mardia, « Algorithm AS 86: The von Mises Distribution Function », Applied Statistics, vol. 24, 1975, p. 268-272
  • (en) Hill, « Algorithm 518, Incomplete Bessel Function I0: The von Mises Distribution », ACM Transactions on Mathematical Software, vol. 3, n° 3, september 1977, p. 279-284
  • (en) D. Best et N. Fisher, « Efficient simulation of the von Mises distribution », Applied Statistics, vol. 28, 1979, p. 152-157
  • (en) M. Evans, N. Hastings et B. Peacock, « von Mises Distribution », ch. 41 in Statistical Distributions, 3e éd., New York, Wiley, 2000
  • (en) Nicholas I. Fisher, Statistical Analysis of Circular Data, New York, Cambridge, 1993
  • (en) M. Evans, N. Hastings et B. Peacock, Statistical Distributions, 2e éd., John Wiley and Sons, 1993 (ISBN 0-471-55951-2), chap. 39
  • (en) Graham Borradaile, Statistics of Earth Science Data : their distribution in time, space, and orientation, Berlin, Springer, , 351 p. (ISBN 978-3-540-43603-4, lire en ligne)
v · m
Lois discrètes
à support fini
0 paramètre de forme
1 paramètre de forme
2 paramètres de forme
3 paramètres de forme
N paramètres de forme
à support infini
0 paramètre de forme
1 paramètre de forme
2 paramètres de forme
3 paramètres de forme
Lois absolument continues
à support compact
0 paramètre de forme
1 paramètre de forme
2 paramètres de forme
3 paramètres de forme
à support semi-infini
0 paramètre de forme
1 paramètre de forme
2 paramètres de forme
3 paramètres de forme
N paramètres de forme
à support infini
0 paramètre de forme
1 paramètre de forme
2 paramètres de forme
Autres types de lois
Lois à support mixte continu-discret
Lois à support variable
Lois multidimensionnelles
Discrètes
Continues
Matricielles
Lois directionnelles
Univariantes
  • von Mises
Sphériques bidimensionnelles
Toroïdales bidimensionnelles
Multidimensionnelles
Lois singulières
Familles de lois
  • icône décorative Portail des probabilités et de la statistique