Maximum-Likelihood-Schätzer für minimale Exponentialverteilungen

10

Ich bin nicht sicher, wie ich dieses Problem lösen soll.

Wir haben also zwei Folgen von Zufallsvariablen, Xi und Yi für i=1,...,n . Nun sind X und Y unabhängige Exponentialverteilungen mit den Parametern λ und μ . Anstatt jedoch die Beobachtung X und Y beobachten wir stattdessen Z und W .

Z=min(Xi,Yi)W=1Zi=XiZi=YiλμZW

Jetzt weiß ich, dass das Minimum von zwei unabhängigen Exponentialen selbst exponentiell ist, wobei die Rate gleich der Summe der Raten ist, also wissen wir, dass mit dem Parameter exponentiell ist . Daher ist unser Maximum-Likelihood-Schätzer: .Zλ+μλ^+μ^=Z¯

Aber ich weiß nicht, wohin ich von hier aus gehen soll. Ich weiß, dass eine Bernoulli-Verteilung mit dem Parameter , aber ich weiß nicht, wie ich dies in eine Aussage über einen der Parameter umwandeln soll. Was würde beispielsweise der MLE in Bezug auf und / oder schätzen ? Ich verstehe, dass wenn , dann , aber es fällt mir schwer, hier herauszufinden, wie ich zu einer algebraischen Aussage kommen kann.Wp=P(Zi=Xi)W¯λμZi=Xiμ=0

UPDATE 1: Daher wurde mir in den Kommentaren gesagt, dass ich die Wahrscheinlichkeit für die gemeinsame Verteilung von und ableiten soll .ZW

Also ist wobei . Richtig? Ich weiß nicht, wie ich in diesem Fall eine gemeinsame Verteilung ableiten soll, da und nicht unabhängig sind.f(Z,W)=f(Z|W=1)p+f(Z|W=0)(1p)p=P(Zi=Xi)ZW

Dies ergibt also nach der obigen Definition vonAber was jetzt? Das bringt mich nicht weiter. Wenn ich die Schritte zur Berechnung der Wahrscheinlichkeit durchführe, erhalte ich: (unter Verwendung von und als Stichprobengröße für jeden Teil der Mischung ...)f(Zi,Wi)=pλeλzi+(1p)μeμziWmn

L(λ,μ)=pmλmeλzi+(1p)nμneμzi

logL=mlogp+mlogλλzi+nlog(1p)+nlogμμzi

Wenn ich die partiellen Ableitungen nehme, sagt mir dies, dass meine MLE-Schätzungen für und nur der Durchschnitt der von abhängigen . Das ist,λμZW

λ^=Zim

μ^=Zin

und

p^=mn+m

Ryan Simmons
quelle
1
Darf ich Sie , nachdem Sie heute eine ähnliche MLE-Frage beantwortet haben , für einige Ideen auf diese Lösung hinweisen? Die Beziehung zwischen den Fragen besteht darin, dass Ihre Daten natürlich auch in zwei disjunkte Gruppen aufgeteilt werden: diejenigen mit und diejenigen mit . Es kommt darauf an, die Wahrscheinlichkeit für eine Beobachtung der Form aufzuschreiben ; Die Symmetrie zwischen und , und , erzeugt sofort die Wahrscheinlichkeit für Daten der Form und dann können Sie loslegen. W=0W=1(Z,W)=(z,0)XYμλ(z,1)
whuber
Beeilen Sie sich nicht, die maximale Wahrscheinlichkeit zu schreiben! zuerst die gemeinsame Verteilung von und leiten Sie dann die Wahrscheinlichkeit ab, die mit der Stichprobe von , die dank der exponentiellen Annahme zufällig geschlossen ist. Dann und nur dann können Sie versuchen, die Funktion zu maximieren und damit die maximale Wahrscheinlichkeit abzuleiten. (Z,W)(Zi,W)=i)
Xi'an
@whuber: (+1) es ist in der Tat ziemlich einfach und beinhaltet die Trennung zwischen und aber beide Gruppen beinhalten sowohl als auch , da sie Informationen über beide liefern und , da . (zi,1)(zi,0) μλ XiYiWi=I(Xi<Yi)
Xi'an
2
@ Xi'an Das stimmt - und die Parallelen zum Normal-Theorie-Beispiel, das ich verlinke, bleiben bestehen, da dort beide Gruppen Informationen über den gemeinsamen Parameter (die Skala) liefern , dessen Schätzung dabei das "Pooling" von Daten beinhaltet aus den Gruppen. Hier ist zu sehen, dass uns sagt, wie die Schätzung von (die Rate oder inverse Skala für ) in separate Schätzungen von und . σW¯λ+μZλμ
whuber
Ich habe den anderen Thread gelesen, whuber, aber ich verstehe ehrlich gesagt nicht, wie ich das auf dieses Beispiel anwenden soll. Z und W sind nicht unabhängig. Wie leite ich die gemeinsame Verteilung ab?
Ryan Simmons

Antworten:

1

Ich habe nicht genug Punkte, um einen Kommentar abzugeben, deshalb werde ich hier schreiben. Ich denke, das Problem, das Sie posten, kann aus der Perspektive einer Überlebensanalyse betrachtet werden, wenn Sie Folgendes berücksichtigen:

Xi : Wahre Überlebenszeit,

Yi : ,

Beide haben eine Exponentialverteilung, wobei und unabhängig sind. Dann ist die beobachtete Überlebenszeit und der Zensurindikator.XYZiWi

Wenn Sie mit der Überlebensanalyse vertraut sind, können Sie meines Erachtens von diesem Punkt aus beginnen.

Anmerkungen: Eine gute Quelle: Analyse von Überlebensdaten durch DRCox und D.Oakes

Unten ist ein Beispiel: Angenommen, das PDF der Überlebenszeitverteilung ist . Dann ist die Überlebensfunktion: . Und die Log-Wahrscheinlichkeit ist:f(t)=ρeρtS(t)=eρt

l=ulogf(zi)+clogS(zi)

mit Summe über unzensierte Personen ( ) bzw. zensierte Personen ( ).uc

Aufgrund der Tatsache, dass wobei h (t) die Gefahrenfunktion ist, kann dies geschrieben werden:f(t)=h(t)S(t)

l=ulogh(zi)+logS(zi)

l=ulogρρzi

Und der Maximum-Likelihood-Schätzer von ist:ρ^ρ

ρ^=d/zi wobei die Gesamtzahl der Fälle vondWi=1

Jujae
quelle