MLE von

7

Sei eine Zufallsstichprobe aus einer Verteilung mit pdf X1,X2,X3,...,Xn

f(x;α,θ)=ex/θθαΓ(α)xα1I(0,)(x),α,θ>0

Finden Sie den Maximum-Likelihood-Schätzer für und . SeiαθΨ(α)=dΓ(α)dα

Mein Versuch,

L(α,θ)=i=1nf(xi)=i=1nexi/θθαΓ(α)xiα1=1Γn(α)θnα(i=1nxi)α1exp(i=1nxiθ)
(α,θ)=nlog(Γ(α))nαlog(θ)+(α1)i=1nlog(xi)1θi=1nxiδ(α,θ)δθ=nαθ+1θ2i=1nxi=01θ2i=1nxi=nαθθ^=i=1nxinα=1αx¯
d(α,θ^)dα=nΓ(α)Γ(α)nlog(1αx¯)+i=1nlog(xi)=0=nΓ(α)Γ(α)+nlog(α)nlog(x¯)+i=1nlog(xi)=0log(α)Γ(α)Γ(α)=log(x¯)i=1nlog(xi)n

Ich konnte nicht mehr fortfahren, um das α . Zweitens weiß ich nicht, wie man Ψ(α)=dΓ(α)dα wie in der Frage angegeben. Hoffe jemand kann es mir erklären.

Danke im Voraus.

Mathxx
quelle
3
Man müsste auf numerische Methoden zurückgreifen, um die MLE von . α
Hartnäckig
Sie schreiben alsΨ(α)Γ(α)
Taylor
1
Könnte Ihnen gesagt worden sein, dass Sie die Digammafunktion anstelle Ihres ? Das würde für mich sinnvoller seinψ(x)=Γ(x)Γ(x)Ψ
am
@Chaconne Ich denke, die Frage hat Tippfehler. Es sollte psi anstelle von Psi sein.
Mathxx
@ Taylor Ich dachte ? ψ(α)=Γ(α)Γ(α)
Mathxx

Antworten:

6

Lassen Sie also ist die Digammafunktion (ich verwende anstelle Ihres ).ψ(α)=Γ(α)Γ(α)ψψΨ

Durch die AM-GM-Ungleichung also (wobei und sind fast sicher definiert). Darüber hinaus gilt die Gleichheit nur für was ein Ereignis mit einer Wahrscheinlichkeit von , sodass fast sicher ist.

x¯(ixi)1/n
logx¯logx¯0
logx¯logxix1==xn0logx¯logx¯>0

Der Einfachheit halber nehme ich .y=logx¯logx¯

Betrachten Sie on . Dies ist stetig und so dass durch den Zwischenwertsatz jede reelle Zahl in trifft . Dies bedeutet insbesondere, dass dh es gibt mindestens einen Punkt in , der , da .f(α)=log(α)ψ(α)(0,)

limα0f(α)=
limαf(α)=0
f(0,)
f1({y})
(0,)yy>0

Darüber hinaus stellt sich heraus, dass auf als injektiv ist so dass es tatsächlich ein eindeutiges mit .f(0,)f<0α^f(α^)=y

Um dieses finden, sind jedoch numerische Methoden erforderlich, wie @StubbornAtom sagt.α^

jld
quelle
4
+1 Schöne Analyse. Angesichts der Asymptotik von und würde ich wärmstens empfehlen, die Lösung von indem die Wurzel von numerisch gefunden wird da dies extrem linear sein wird (mit Steigung bei großen Werten). Γψ,log(x)ψ(x)=C
1log(x)ψ(x)1C,
2
whuber