4 svar
223 visningar
beep är nöjd med hjälpen
beep 21 – Fd. Medlem
Postad: 13 nov 2019 17:20

Maximum Likelihood Skattning

Hej,

En stokastisk variabel har fördelningsfunktion 

F(x)=xa, 0x1

där a>0.

Man gör 10 oberoende observationer av variabeln och erhåller

0.57   0.81   0.63   0.44   0.31   0.91   0.36   0.65   0.74   0.99 (Summa = 6.41)

ML-skatta a!

Försök till lösning

Täthetsfunktionen ges av 

f(x)=axa-1

Och ML funktionen blir då

L(a)=i=110  a xia-1.

Sedan vet jag inte riktigt hur jag ska fortsätta, har testat två olika "vägar", en av dem är att gå jobba med logaritmlikelihoodfunktionen 

l(a)=lnL(a)=ln(i=110  a xia-1) =10ln(a)+(a-1)ln(xii=110)

och sedan derivera. Men då utnyttjar jag inte observationsvärdena eller att summan av dem blir 6.41. Något förslag på hur jag ska fortsätta? Tack

parveln 703 – Fd. Medlem
Postad: 13 nov 2019 18:03

Vad menar du med att du inte utnyttjar observationsvärdena? Observationsvärdena är det du skrivit som x_i

beep 21 – Fd. Medlem
Postad: 13 nov 2019 18:33
parveln skrev:

Vad menar du med att du inte utnyttjar observationsvärdena? Observationsvärdena är det du skrivit som x_i

Juste, jag blandade ihop funktionsvärdet f(x)med observationsvärdet xdärav

min kommentar. Jag vet inte hur jag ska fortsätta eftersom derivatan

l'(a)=10/a+ln(6.41)

inte kan bli noll om a>0.

Arktos 3927
Postad: 13 nov 2019 18:37 Redigerad: 13 nov 2019 20:05

Kolla derivatan en gång till.
(EDIT  Jag menade kolla  log-likelihood-funktionen!)
Borde det inte stå summa log av observationsvärdena?

beep 21 – Fd. Medlem
Postad: 13 nov 2019 18:59
Arktos skrev:

Kolla derivatan en gång till.
Borde det inte stå summa log av observationsvärdena?

Juste, gjorde istället tvärtom. Det ska alltså bli

l(a)=10ln(a)+(a-1)i=1nln(xi)
och derivatan blir då

l'(a)=10/a+inln(xi).

Räknade även ut summan av log av observationerna (den sista termen) som summerades till -5.0849. 

Då får vi att derivatan är noll när a är lika med 1.9665, vilket överensstämmer med svaret i facit. 

Tackar för hjälpen!

Svara Avbryt
Close