- Information de Fisher
-
L'information de Fisher est une notion de statistique introduite par R.A. Fisher qui quantifie l'information relative à un paramètre contenue dans une distribution.
Soit f(x;θ) la distribution de vraisemblance d'une grandeur x (qui peut être multidimensionelle), paramétrée par θ. La technique d'estimation de θ par le maximum de vraisemblance, introduite par Fisher consiste à choisir la valeur maximisant la vraisemblance des observations X :
. L'information de Fisher est quant à elle définie comme la variance associée à ce maximum :
. Sommaire
Formulation discrète
Les différentes observations xi nous permettent d'échantillonner la fonction de densité de probabilité f(x;θ). Selon le théorème de Bayes, en l'absence d'a priori sur θ on a Si les observations sont décorrélées, la valeur la plus probable nous est donnée par le maximum de
∏ P(xi / θ), i qui est aussi le maximum de
λ(θ) = ∑ log P(xi / θ). i Le passage en logarithme permet de transformer le produit en somme, ce qui nous autorise à trouver le maximum par dérivation :
Cette somme correspond pour un nombre d'observations suffisamment élevé à l'espérance mathématique. La résolution de cette équation permet de trouver un estimateur de θ à partir du jeu de paramètre au sens du maximum de vraisemblance. Maintenant, la question est de quantifier la précision de notre estimation. On cherche donc à estimer la forme de la distribution de probabilité de θ autour de la valeur donnée par l'estimateur . À partir d'un développement limité à l'ordre 2, comme le terme linéaire est nul au maximum, on obtient :
où est l'information de Fisher relative à θ au point de maximum de vraisemblance. Ceci signifie que la distribution est en première approximation une gaussienne de variance :
Cette variance est appelé la borne de Cramér-Rao et constitue la meilleure précision d'estimation atteignable en absence d'a priori.
Additivité
Une des propriétés fondamentales de l'information de Fisher est son additivité. L'information résultant de deux variables aléatoires indépendantes est la somme des informations :
IX,Y(θ) = IX(θ) + IY(θ). Si on a N réalisations indépendantes obéissant une même densité de probabilité, l'information résultante est une simple mise à l'échelle de l'information individuelle.
Lorsque une statistique S(X) sur une variable aléatoire X est exhaustive, l'information relative à la statistique est égale à celle de la variable aléatoire. Autrement dit
avec égalité pour une statistique exhaustive.
Formulation multi-paramétrique
Dans le cas où la distribution de probabilité f(X) dépend de plusieurs paramètres, θ n'est plus un scalaire mais un vecteur . La recherche du maximum de vraisemblance ne se résume donc non pas à une seule équation mais à un système :
on dérive vis-à-vis des différentes composantes de . Enfin, l'information de Fisher n'est plus définie comme une variance scalaire mais comme une matrice de covariance :
Estimation et borne de Cramér-Rao
L'inverse de cette matrice permet quant à elle de déterminer les bornes de Cramér-Rao, i.e. les covariances relatives aux estimations conjointes des différents paramètres à partir des observations : en effet, le fait que tous les paramètres soient à estimer simultanément rend l'estimation plus difficile. Ce phénomène est une manifestation de ce qui est parfois appelé le « fléau de la dimension ». C'est pour cette raison que l'on utilise quand on le peut des a priori sur les paramètres (méthode d'estimation du maximum a posteriori). Ainsi, on restreint l'incertitude sur chacun des paramètres, ce qui limite l'impact sur l'estimation conjointe.
Métrique de Fisher
Cette matrice est couramment appelée la métrique d'information de Fisher. En effet, le passage de l'espace des observations à l'espace des paramètres est un changement de système de coordonnées. Dans la base des paramètres, avec comme produit scalaire la covariance, cette matrice est la métrique. Ce point de vue géométrique, introduit par C. Rao, a été ensuite largement développé par S. Amari sous la dénomination de géométrie de l'information. L'information de Fisher joue un rôle particulier en tant que métrique de part ses propriétés d'additivité et d'invariance par rapport à l'échantillonnage statistique (théorème de Chentsov ou Čencov). C'est une métrique qui est donc naturelle lorsque sont considérées des distributions de probabilité. De plus, l'approche du concept d'information sous l'angle de la géométrie différentielle permet de proposer un cadre cohérent liant différents concepts :
- divergence de Kullback-Leibler,
- entropie et principe d'entropie maximale,
- famille exponentielle des distribution,
- algorithme espérance-maximisation,
- estimation par maximum de vraisemblance.
Formulations alternatives
Il existe un certain nombre de formulations alternatives de l'information de Fisher révélant certaines propriétés intéressantes.
Information apportée par une statistique
De la même façon que l'on a défini l'information de Fisher pour le vecteur des observations X on peut définir l'information de Fisher contenue dans une statistique S(X):
Cette définition est exactement la même que celle de l'information de Fisher pour X pour un modèle multiparamétrique on remplace juste la densité de X par celle de S(X) la statistique S. Deux théorèmes illustrent l'intérêt de cette notion:
- Pour une statistique exhaustive on a IS(θ) = I(θ) ce qui permet de voir une statistique exhaustive comme une statistique comprenant toute l'information du modèle. L'on a aussi la réciproque à savoir que si IS(θ) = I(θ) alors S est exhaustif bien que cette caractérisation est rarement utilisée dans ce sens la définition grâce au critère de factorisation des statistiques exhaustives étant souvent plus maniable.
- Quelle que soit la statistique S, avec un cas d'égalité uniquement pour des statistiques exhaustives. On ne peut donc récupérer plus d'information que celle contenue dans une statistique exhaustive. Ceci explique en grande partie l'intérêt des statistiques exhaustives pour l'estimation. La relation d'ordre est ici la relation d'ordre partielle sur les matrices symétriques à savoir qu'une matrice si B-A est une matrice symétrique positive.
Références
- A. Monfort, Cours de statistique mathématique, 1982, Economica. Paris.
- C.R. Rao, Information and accuracy attainable in the estimation of statistical parameters, Bulletin of the Calcutta Mathematical Society, 37:81-91, 1945.
- S. Amari, H. Nagaoka, Methods of information geometry, Translations of mathematical monographs; v. 191, American Mathematical Society, 2000 (ISBN 978-0821805312).
- N. N. Chentsov (Čencov), Statistical Decision Rules and Optimal Inference, Translations of Mathematical Monographs; v. 53, American Mathematical Society, 1982.
Voir aussi
Liens externes
- P. Druilhet [1] Cours de statistique inférentielle.
- Portail des probabilités et des statistiques
Catégorie :- Estimation (statistique)
Wikimedia Foundation. 2010.