Ajustement de loi de probabilité
L'ajustement de la loi de probabilité ou simplement l'ajustement de la loi est l'ajustement d'une loi de probabilité à une série de données concernant la mesure répétée d'un phénomène aléatoire. L'ajustement de la loi a pour but de prédire la probabilité ou de prévoir la fréquence d'occurrence de l'ampleur du phénomène dans un certain intervalle.
Il existe de nombreuses lois de probabilité, dont certaines peuvent être ajustées plus étroitement à la fréquence observée des données que d'autres, selon les caractéristiques du phénomène et de la loi. La loi donnant un ajustement serré est supposée conduire à de bonnes prédictions. Dans l'ajustement, il faut donc sélectionner une famille de lois qui convient bien aux données.
Sélection de la loi
modifier![](http://upload.wikimedia.org/wikipedia/commons/thumb/7/74/Normal_Distribution_PDF.svg/220px-Normal_Distribution_PDF.svg.png)
Le choix de la loi appropriée dépend de la présence ou de l'absence de symétrie de l'ensemble de données par rapport à la tendance centrale.
Lois symétriques
modifierLorsque les données sont réparties symétriquement autour de la moyenne alors que la fréquence d'occurrence des données plus éloignées de la moyenne diminue, on peut par exemple sélectionner la loi normale, la loi logistique ou la loi de Student. Les deux premières sont très similaires, tandis que la dernière, avec un degré de liberté, a des "queues plus lourdes", ce qui signifie que les valeurs les plus éloignées de la moyenne se produisent relativement plus souvent (c'est-à-dire que le kurtosis est plus élevé). La loi de Cauchy est également symétrique.
Lois asymétriques à droite
modifier![](http://upload.wikimedia.org/wikipedia/commons/thumb/f/f8/Negative_and_positive_skew_diagrams_%28English%29.svg/220px-Negative_and_positive_skew_diagrams_%28English%29.svg.png)
Lorsque les valeurs les plus grandes ont tendance à être plus éloignées de la moyenne que les valeurs les plus petites, on a une loi asymétrique vers la droite (c'est-à-dire qu'il y a une asymétrie positive), on peut par exemple sélectionner la loi log-normale (c'est-à-dire les logarithmes des données sont normalement distribuées), la loi log-logistique (c'est-à-dire que les logarithmes des données suivent une loi logistique), la loi de Gumbel, la loi exponentielle, la loi de Pareto, la loi de Weibull, la loi de Burr ou la loi de Fréchet. Les quatre dernières lois sont bornées à gauche.
Lois asymétriques à gauche
modifierLorsque les petites valeurs ont tendance à être plus éloignées de la moyenne que les grandes valeurs, on a une loi asymétrique vers la gauche (c'est-à-dire qu'il y a une asymétrie négative), on peut par exemple sélectionner la loi du χ2
1 (c'est-à-dire le carré des données suit une loi normale)[1], la loi de Gumbel réfléchie[1], la loi de Dagum (loi de Burr en miroir) ou la loi de Gompertz, qui est limitée à gauche.
Techniques d'ajustement
modifierLes techniques d'ajustement de loi suivantes existent[2] :
- Méthodes paramétriques, par lesquelles les paramètres de la loi sont calculés à partir des séries de données[3]. Les méthodes paramétriques sont :
- Méthode des moments
- Estimation de l'espacement maximal (en)
- Méthode des L-moments (en)[4]
- Méthode du maximum de vraisemblance[5]
![](http://upload.wikimedia.org/wikipedia/commons/thumb/9/9e/FitGumbelDistr.tif/lossless-page1-220px-FitGumbelDistr.tif.png)
- Diagramme quantile-quantile puis analyse de régression, utilisant une transformation de la fonction de répartition de sorte qu'une relation linéaire soit trouvée entre la répartition et les valeurs des données, qui peuvent également devoir être transformées, selon la loi de probabilité sélectionnée. Dans cette méthode, la répartition doit être estimée par la position du tracé[6]
Généralisation des lois
modifierIl est d'usage de transformer les données de manière logarithmique pour ajuster des lois symétriques (comme la normale et la logistique) aux données obéissant à une loi positivement asymétrique (c'est-à-dire asymétrique vers la droite, avec une moyenne supérieure au mode, et avec une queue à droite plus longue que la queue de gauche), voir la loi log-normale et la loi log-logistique. Un effet similaire peut être obtenu en prenant la racine carrée des données.
Pour ajuster une loi symétrique aux données obéissant à une loi biaisée négativement (c'est-à-dire biaisée vers la gauche, avec une moyenne inférieure mode, et avec une queue droite plus courte que la queue gauche), on pourrait utiliser les valeurs au carré des données pour accomplir l'ajustement.
Plus généralement, on peut élever les données à une puissance p afin d'ajuster des lois symétriques à des données obéissant à une loi d'asymétrie quelconque, où p < 1 lorsque l'asymétrie est positive et p > 1 lorsque l'asymétrie est négative. La valeur optimale de p doit être trouvée par une méthode numérique . La méthode numérique peut consister à supposer une plage de valeurs p, puis à appliquer la procédure d'ajustement de loi à plusieurs reprises pour toutes les valeurs p supposées, et enfin à sélectionner la valeur de p pour laquelle la somme des carrés des écarts des probabilités calculées à partir des fréquences mesurées (test du χ²) est minimum, comme c'est le cas dans CumFreq.
La généralisation améliore la flexibilité des lois de probabilité et augmente leur applicabilité dans l'ajustement de loi[6].
La polyvalence de la généralisation permet, par exemple, d'adapter des ensembles de données distribués approximativement normalement à un grand nombre de lois de probabilité différentes, tandis que des lois asymétriques négatives peuvent être ajustées à des lois de Gumbel invesée et du χ2
1.
Inversion de l'asymétrie
modifier![](http://upload.wikimedia.org/wikipedia/commons/thumb/a/a9/Gumbel_distribution_and_Gumbel_mirrored.png/220px-Gumbel_distribution_and_Gumbel_mirrored.png)
Les lois asymétriques peuvent être réfléchies en remplaçant dans l'expression mathématique de la fonction de répartition (F) par son complément : F' = 1–F, obtenant la fonction de répartition complémentaire (également appelée fonction de survie) qui donne une image miroir. De cette manière, une loi asymétrique vers la droite est transformée en une loi asymétrique vers la gauche et vice versa.
La technique d'inversion d'asymétrie augmente le nombre de lois de probabilité disponibles pour l'ajustement de loi et élargit les opportunités d'ajustement de loi.
Déplacement des lois
modifierCertaines lois de probabilité, comme l'exponentielle, ne prennent pas en charge les valeurs de données (X) égales ou inférieures à zéro. Pourtant, lorsque des données négatives sont présentes, de telles lois peuvent toujours être utilisées en remplaçant X par Y = X – Xm, où Xm est la valeur minimale de X. Ce remplacement représente un déplacement de la lois de probabilité dans le sens positif, c'est-à-dire vers la droite, car Xm est négatif. Après avoir terminé l'ajustement de la loi de X, les valeurs X correspondantes sont trouvées à partir de X = Y + Xm, ce qui représente un décalage arrière de la loi dans le sens négatif, c'est-à-dire vers la gauche. La technique de déplacement de loi augmente les chances de trouver une loi de probabilité bien ajustée.
Lois composites
modifier![](http://upload.wikimedia.org/wikipedia/commons/thumb/3/3a/SanLor.jpg/220px-SanLor.jpg)
L'option existe d'utiliser deux loi de probabilité différentes, une pour la plage de données inférieure et une pour la plage supérieure, comme la loi de Laplace. Les plages sont séparées par un point d'arrêt. L'utilisation de telles loi de probabilité composites (discontinues) peut être opportune lorsque les données du phénomène étudié ont été obtenues dans deux ensembles de conditions différentes[6].
Incertitude de prédiction
modifier![](http://upload.wikimedia.org/wikipedia/commons/thumb/0/0f/BinomialConfBelts.jpg/220px-BinomialConfBelts.jpg)
Les prévisions d'occurrence basées sur des lois de probabilité ajustées sont sujettes à l'incertitude, qui découle des conditions suivantes :
- La vraie loi de probabilité des événements peut s'écarter de la loi ajustée, car la série de données observées peut ne pas être totalement représentative de la probabilité réelle d'occurrence du phénomène en raison d'une erreur aléatoire
- L'occurrence d'événements dans une autre situation ou dans le futur peut s'écarter de la loi ajustée car cette occurrence peut également être sujette à une erreur aléatoire
- Un changement des conditions environnementales peut entraîner une modification de la probabilité d'occurrence du phénomène
![](http://upload.wikimedia.org/wikipedia/commons/thumb/5/51/SampleFreqCurves.tif/lossy-page1-220px-SampleFreqCurves.tif.jpg)
Une estimation de l'incertitude dans le premier et le second cas peut être obtenue avec la loi de probabilité binomiale en utilisant par exemple la probabilité de dépassement pe (c'est-à-dire la chance que l'événement X soit supérieur à une valeur de référence Xr de X) et la probabilité de non-dépassement pn (c'est-à-dire la probabilité que l'événement X soit inférieur ou égal à la valeur de référence Xr, on l'appelle aussi répartition). Dans ce cas, il n'y a que deux possibilités : soit il y a dépassement, soit il y a non-dépassement. Cette dualité est la raison pour laquelle la loi binomiale est applicable.
Avec la loi binomiale, on peut obtenir un intervalle de prédiction. Un tel intervalle estime également le risque d'échec, c'est-à-dire la probabilité que l'événement prédit reste toujours en dehors de l'intervalle de confiance. L'analyse de confiance ou de risque peut inclure la période de retour T = 1/pe comme cela se fait en hydrologie.
![](http://upload.wikimedia.org/wikipedia/commons/thumb/f/f6/CumList.png/220px-CumList.png)
![](http://upload.wikimedia.org/wikipedia/commons/thumb/1/17/GEVdistrHistogr%2BDensity.png/220px-GEVdistrHistogr%2BDensity.png)
Qualité de l'ajustement
modifierEn classant la qualité de l'ajustement des différentes lois, on peut se faire une idée de la loi qui est acceptable et de celle qui ne l'est pas.
Histogramme et fonction de densité
modifierÀ partir de la fonction de répartition, on peut dériver un histogramme et la fonction de densité.
Notes et références
modifier- (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Probability distribution fitting » (voir la liste des auteurs).
- Les histogrammes de fréquences asymétriques à gauche (négativement) peuvent être ajustés à une loi du χ2
1 ou de Gumbel réfléchie. - (en) H.P.Ritzema, Drainage Principles and Applications, vol. 16, Wageningen, International Institute for Land Reclamation and Improvement (ILRI), , 175–224 p. (ISBN 9070754339), « Frequency and Regression Analysis »
- (en) H. Cramér, Mathematical methods of statistics, Princeton Univ. Press,
- Hosking, « L-moments: analysis and estimation of distributions using linear combinations of order statistics », Journal of the Royal Statistical Society, Series B, vol. 52, , p. 105–124 (JSTOR 2345653)
- (en) Aldrich, « R. A. Fisher and the making of maximum likelihood 1912–1922 », Statistical Science, vol. 12, no 3, , p. 162–176 (DOI 10.1214/ss/1030037906, MR 1617519)
- Software for Generalized and Composite Probability Distributions. International Journal of Mathematical and Computational Methods, 4, 1-9
- Intro to composite probability distributions
- Frequency predictions and their binomial confidence limits. In: International Commission on Irrigation and Drainage, Special Technical Session: Economic Aspects of Flood Control and non-Structural Measures, Dubrovnik, Yugoslavia, 1988. On line
- M.A. Benson, Flood frequency analysis, vol. 1543-A, T. Dalrymple, , p. 51-71, « Characteristics of frequency curves based on a theoretical 1000 year record ».