Stima della massima verosimiglianza per il minimo delle distribuzioni esponenziali


10

Sono bloccato su come risolvere questo problema.

Quindi, abbiamo due sequenze di variabili casuali, e per . Ora, e sono distribuzioni esponenziali indipendenti con parametri e . Tuttavia, invece di osservare e , osserviamo invece e .Y i i = 1 , . . . , n X Y λ μ X Y Z WXiYii=1,...,nXYλμXYZW

Z=min(Xi,Yi) e W=1 se Zi=Xi e 0 se Zi=Yi . Devo trovare chiuso le forme per gli stimatori di massima verosimiglianza dei λ e μ sulla base di Z e W . Inoltre, dobbiamo dimostrare che si tratta di massimi globali.

Ora so che il minimo di due esponenziali indipendenti è esso stesso esponenziale, con il tasso uguale alla somma dei tassi, quindi sappiamo che Z è esponenziale con il parametro λ+μ . Quindi il nostro stimatore della massima verosimiglianza è: λ^+μ^=Z¯ .

Ma sono bloccato con dove andare da qui. So che W è una distribuzione di Bernoulli con il parametro p=P(Zi=Xi) , ma non so come fare per convertirlo in un'istruzione su uno dei parametri. Ad esempio, cosa stima l'MLE W¯ in termini di λ e / o μ ? Capisco che se Zi=Xi , allora μ=0 , ma sto facendo fatica a capire come elaborare qualsiasi affermazione algebrica, qui.

UPDATE 1: Così mi è stato detto nei commenti per ricavare la probabilità per la distribuzione congiunta di Z e W .

Quindi dove . Corretta? Non so in che altro modo derivare una distribuzione congiunta in questo caso, poiché e non sono indipendenti.f(Z,W)=f(Z|W=1)p+f(Z|W=0)(1p)p=P(Zi=Xi)ZW

Quindi questo ci dà, , secondo la definizione di sopra. Ma adesso cosa? Questo non mi porta da nessuna parte. Se passo attraverso i passaggi del calcolo della probabilità, ottengo: (usando e come dimensioni del campione per ogni parte della miscela ...)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

Se prendo le derivate parziali, questo mi dice che il mio MLE stime per e sono solo la media del 's subordinata . Questo è,λμZW

λ^=Zim

μ^=Zin

e

p^=mn+m


1
Avendo appena risposto a una domanda MLE simile oggi, posso indirizzarti verso quella soluzione per alcune idee? La relazione tra le domande è che i tuoi dati si suddividono naturalmente anche in due gruppi disgiunti: quelli in cui e quelli in cui . Tutto si riduce a scrivere la probabilità di un'osservazione della forma ; la simmetria tra e , e , produce immediatamente la probabilità per i dati del modulo e poi sei fuori e corri. W=0W=1(Z,W)=(z,0)XYμλ(z,1)
whuber

Non abbiate fretta di scrivere la massima probabilità! Innanzitutto, esprimere la distribuzione congiunta di , quindi dedurre la probabilità associata al campione di , che risulta essere in forma chiusa grazie all'assunto esponenziale. Quindi e solo allora puoi provare a massimizzare la funzione e quindi derivare la massima probabilità. (Z,W)(Zi,W)=i)
Xi'an,

@whuber: (+1) è piuttosto semplice e implica la separazione tra e ma entrambi i gruppi coinvolgono sia che , poiché forniscono informazioni su entrambi e , poiché . (zi,1)(zi,0) μλ XiYiWi=I(Xi<Yi)
Xi'an,

2
@ Xi'an Esatto - e i parallelismi con l'esempio di teoria normale che collego per continuare a tenere, perché entrambi i gruppi forniscono informazioni sul parametro comune (la scala), la cui stima implicherà quindi il "raggruppamento" dei dati dai gruppi. Qui si vedrà che ci dice come la stima di (il tasso, o scala inversa, per ) dovrebbe essere suddivisa in stime separate di e . σW¯λ+μZλμ
whuber

Ho letto l'altro thread, whuber, ma sinceramente non capisco come applicarlo a questo esempio. Z e W non sono indipendenti, quindi come posso derivare la distribuzione congiunta?
Ryan Simmons,

Risposte:


1

Non ho abbastanza punti per commentare, quindi scriverò qui. Penso che il problema che pubblichi possa essere visualizzato dal punto di vista dell'analisi di sopravvivenza, se si considera quanto segue:

Xi : vero tempo di sopravvivenza,

Yi : tempo di censura,

Entrambi hanno una distribuzione esponenziale con e indipendenti. Quindi è il tempo di sopravvivenza osservato e l'indicatore di censura.XYZiWi

Se hai familiarità con l'analisi della sopravvivenza, credo che tu possa iniziare da questo punto.

Note: Una buona fonte: analisi dei dati di sopravvivenza di DRCox e D.Oakes

Di seguito è riportato un esempio: supponendo che il pdf della distribuzione del tempo di sopravvivenza sia . Quindi la funzione di sopravvivenza è: . E la probabilità logaritmica è:f(t)=ρeρtS(t)=eρt

l=ulogf(zi)+clogS(zi)

con la somma delle persone senza censura ( ) e delle persone censurate ( ) rispettivamente.uc

A causa del fatto che dove h (t) è la funzione di rischio, questo può essere scritto:f(t)=h(t)S(t)

l=ulogh(zi)+logS(zi)

l=ulogρρzi

E lo stimatore di massima verosimiglianza di è:ρ^ρ

ρ^=d/zi dove è il numero totale di casi didWi=1

Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.