Estimateur du maximum de vraisemblance - intervalle de confiance

12

Comment puis-je construire un intervalle de confiance asymptotique pour un paramètre réel, à partir du MLE pour ce paramètre?

Fabrizio
la source
Une façon d'aborder ce problème est d'utiliser la méthode delta: en.wikipedia.org/wiki/Delta_method
J'ai remarqué qu'il y a des votes pour clore cette question comme trop large, mais il y a un théorème général sur le comportement asymptotique des MLE qui peut être résumé succinctement. J'ai mis une réponse concise que je développerai un peu plus tard.
Scortchi - Réintégrer Monica

Réponses:

13

Utilisez le fait que pour un échantillon iid de taille , étant donné certaines conditions de régularité, le MLE est un estimateur cohérent du vrai paramètre , et de sa distribution asymptotiquement normale, avec une variance déterminée par l'inverse de la Informations sur le pêcheur:& thetav & thetav 0nθ^θ0

I 1 ( θ 0 )I( θ )

n(θ^θ0)N(0,1I1(θ0))
où est l'information Fisher d'un seul échantillon. Les informations observées au MLE tendent asymptotiquement aux informations attendues, vous pouvez donc calculer (disons 95%) des intervalles de confiance avecI1(θ0)I(θ^)

θ^±1.96nI1(θ^)

Par exemple, si est une variable de Poisson tronquée, vous pouvez obtenir une formule pour les informations observées en termes de MLE (que vous devez calculer numériquement): X

f(x)=eθθxx!(1eθ)

(θ)=θ+xlogθlog(1eθ)

d(θ)dθ=1+xθeθ1eθ

I1(θ^)=d2(θ^)(dθ^)2=xθ^eθ^(1eθ^)2

Les cas notables exclus par les conditions de régularité comprennent ceux où

  • le paramètre détermine le support des données, par exemple l'échantillonnage à partir d'une distribution uniforme entre rien etθθ
  • le nombre de paramètres de nuisance augmente avec la taille de l'échantillon
Scortchi - Réintégrer Monica
la source
Cette méthode s'applique-t-elle sans modification lorsqu'il existe des contraintes sur , par exemple ? Qu'en est-il d'un MLE pour paramètres , tels que et ? θθ[0,1]Nθii=0,...,N1i=0N1θi=1θi[0,1]
quant_dev
1
Si , c'est-à-dire que la vraie valeur n'est pas égale à l'une des bornes. θ(0,1)
Scortchi - Réintégrer Monica
Si et, cela ne signifierait-il pas que l'approximation normale n'est pas applicable et j'ai besoin de plus d'échantillons? θ(0,1)σ(θ^)>|θ^|
quant_dev
Oui, ce n'est qu'un intervalle de confiance asymptotique.
Scortchi - Réintégrer Monica
1
@quant_dev: Non: vous chercheriez une transformation du ou des paramètres qui ont rendu l'approximation normale décente - ou utiliser une autre méthode.
Scortchi - Réintégrer Monica