Maximum Likelihood Estimator - Konfidenzintervall

11

Wie kann ich ein asymptotisches Konfidenzintervall für einen realen Parameter erstellen, beginnend mit dem MLE für diesen Parameter?

Fabrizio
quelle
Eine Möglichkeit, dieses Problem
Mir ist aufgefallen, dass es Stimmen gibt, um diese Frage als zu weit gefasst zu schließen, aber es gibt einen allgemeinen Satz über das asymptotische Verhalten von MLEs, der kurz und bündig formuliert werden kann. Ich habe eine knappe Antwort gegeben, die ich etwas später erweitern werde.
Scortchi - Monica wieder einsetzen

Antworten:

13

Verwenden Sie die Tatsache, dass für eine iid-Stichprobe der Größe unter bestimmten Regelmäßigkeitsbedingungen der MLE ein konsistenter Schätzer des wahren Parameters und seiner asymptotisch normalen Verteilung ist, wobei die Varianz durch den Kehrwert von bestimmt wird Fischerinformation:nθ^θ0

n(θ^θ0)N(0,1I1(θ0))
wobei die Fisher-Information aus einer einzelnen Stichprobe ist. Die beobachteten Informationen im MLE tendieren asymptotisch zu den erwarteten Informationen, sodass Sie Konfidenzintervalle (z. B. 95%) mit berechnen könnenI1(θ0)I(θ^)

θ^±1.96nI1(θ^)

Wenn beispielsweise eine Poisson-Variable mit Null-Kürzung ist, können Sie eine Formel für die beobachteten Informationen in Bezug auf die MLE erhalten (die Sie numerisch berechnen müssen): X

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

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

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

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

Bemerkenswerte Fälle, die von den Regelmäßigkeitsbedingungen ausgeschlossen sind, umfassen solche, in denen

  • Der Parameter bestimmt die Unterstützung der Daten, z. B. Stichproben aus einer gleichmäßigen Verteilung zwischen nichts undθθ
  • Die Anzahl der Störparameter nimmt mit der Stichprobengröße zu
Scortchi - Monica wieder einsetzen
quelle
Gilt diese Methode unverändert, wenn es Einschränkungen für , z. B. ? Was ist mit einem MLE für Parameter , so dass und ? θθ[0,1]Nθii=0,...,N1i=0N1θi=1θi[0,1]
quant_dev
1
Wenn , dh der wahre Wert ist nicht gleich einer der Grenzen. θ(0,1)
Scortchi - Monica wieder einsetzen
Wenn undWürde das nicht bedeuten, dass die normale Annäherung nicht anwendbar ist und ich mehr Proben benötige? θ(0,1)σ(θ^)>|θ^|
quant_dev
Ja, es ist nur ein asymptotisches Konfidenzintervall.
Scortchi - Monica wieder einsetzen
1
@quant_dev: Nein: Sie würden nach einer Transformation der Parameter suchen, die die normale Näherung anständig gemacht haben - oder eine andere Methode verwenden.
Scortchi - Monica wieder einsetzen