M-estimateur

classe d'estimateurs statistitiques

En statistique, les M-estimateurs constituent une large classe de statistiques obtenues par la minimisation d'une fonction dépendant des données et des paramètres du modèle. Le processus du calcul d'un M-estimateur est appelé M-estimation. De nombreuses méthodes d'estimation statistiques peuvent être considérées comme des M-estimateurs. Dépendant de la fonction à minimiser lors de la M-estimation, les M-estimateurs peuvent permettre d'obtenir des estimateurs plus robustes que les méthodes plus classiques, comme la méthode des moindres carrés.

M-estimateur

Définition

modifier

Les M-estimateurs ont été introduits en 1964 par Peter Huber sous la forme d'une généralisation de l'estimation par maximum de vraisemblance à la minimisation d'une fonction ρ sur l'ensemble des données. Ainsi, le (ou les) M-estimateur associé aux données et à la fonction ρ est estimé par

θ̂=argminθ(∑i=1nρ(xi,θ))

Le M de M-estimateur provient donc du maximum de vraisemblance (maximum likelihood-type en anglais) et les estimateurs par maximum de vraisemblance sont un cas particulier des M-estimateurs.

La résolution du problème de minimisation passe couramment par une différentiation de la fonction cible. En effet, pour chercher θ̂, une méthode simple consiste à chercher les valeurs telles que

∂∂θ(∑i=1nρ(xi,θ))=0.

Dans le cas où cette différentiation est possible, le M-estimateur est dit de type ψ ; sinon, il est dit de type ρ.

Type ρ

modifier

Pour un entier positif r, soit (𝒳,Σ) et (Θ⊂ℝr,S) des espaces mesurables, et θ∈Θ est un vecteur de paramètres. Un M-estimateur de type ρ T est défini par une fonction mesurable ρ:𝒳×Θ→ℝ. Il envoie la fonction de répartition d'une loi de probabilités F sur 𝒳 vers la valeur T(F)∈Θ (si elle existe) qui minimise ∫𝒳ρ(x,θ)dF(x), soit :

T(F):=arg⁡minθ∈Θ∫𝒳ρ(x,θ)dF(x)

Par exemple, la fonction qui correspond à l'estimateur du maximum de vraisemblance et ρ(x,θ)=−log⁡(f(x,θ)), avec f(x,θ)=∂F(x,θ)∂x.

Type ψ

modifier

Si ρ est dérivable par rapport à θ, le calcul de θ^ est souvent plus simple. Un M-estimateur de type ψ T est défini pour une fonction mesurable ψ:𝒳×Θ→ℝr. Il envoie la fonction de répartition d'une loi de probabilités F sur 𝒳 vers la valeur T(F)∈Θ (si elle existe) qui vérifie l'équation vectorielle:

∫𝒳ψ(x,θ)dF(x)=0
∫𝒳ψ(x,T(F))dF(x)=0

Par exemple, pour l'estimateur du maximum de vraisemblance, on a ψ(x,θ)=(∂log⁡(f(x,θ))∂θ1,…,∂log⁡(f(x,θ))∂θp)T, avec uT désignant le vecteur transposé de u et f(x,θ)=∂F(x,θ)∂x.

Un tel estimateur n'est pas nécessairement un M-estimateur de type ρ, mais si ρ est dérivable et de dérivée continue par rapport à θ, alors une condition nécessaire pour qu'un M-estimateur de type ψ soit un M-estimateur de type ρ est ψ(x,θ)=∇θρ(x,θ). De telles définitions peuvent être adaptés aux échantillons finis.

Si la fonction ψ décroit vers zéro quand x→±∞, l'estimateur est dit décroissant. De tels estimateurs ont d'autres propriétés utiles, comme le rejet des valeurs aberrantes évidentes.

Exemples de M-estimateurs

modifier

Parmi les exemples connus de M-estimateurs, on peut citer :

Articles connexes

modifier

Références

modifier