Metoda Największej Wiarygodności i rozkład Poissona

Procesy stochastyczne. Sposoby racjonalizowania wielkich ilości informacji. Matematyka w naukach społecznych.
mac18
Użytkownik
Użytkownik
Posty: 316
Rejestracja: 3 wrz 2013, o 19:09
Płeć: Mężczyzna
Lokalizacja: Warszawa
Podziękował: 88 razy
Pomógł: 23 razy

Metoda Największej Wiarygodności i rozkład Poissona

Post autor: mac18 » 31 paź 2017, o 17:15

Hej, mam problem ze zwinięciem funkcji wiarygodności próby z rozkładu Poissona.
Mam:
\(\displaystyle{ L\left( x _{1},..., x_{n}; \lambda \right) = \prod_{i=1}^{n} \frac{\lmbda ^{x _{i} } }{ x_{i}!} e ^{-\lambda}}\)
dostaję:
\(\displaystyle{ L\left( x _{1},..., x_{n}; \lambda \right) = e^{-n \lambda} \lambda^{ \sum_{i=1}^{n} x _{i}} ....}\)
i nie wiem tylko jak odpowiednio zwinąć sumę silni. Dalej sobie poradzę.

Edit:
a chyba wiem, trzeba w iloczyn:
\(\displaystyle{ L\left( x _{1},..., x_{n}; \lambda \right) = e^{-n \lambda} \cdot \lambda^{ \sum_{i=1}^{n} x _{i} } \cdot \frac{1}{ \prod_{i=1}^{n} x_{i}! }}\)

Awatar użytkownika
Premislav
Użytkownik
Użytkownik
Posty: 15211
Rejestracja: 17 sie 2012, o 13:12
Płeć: Mężczyzna
Lokalizacja: Warszawa
Podziękował: 161 razy
Pomógł: 5046 razy

Re: Metoda Największej Wiarygodności i rozkład Poissona

Post autor: Premislav » 31 paź 2017, o 17:41

Jest OK. Niemniej jednak pozwolę sobie polecić zlogarytmowanie tego, wówczas obliczenia będą przyjemniejsze (zaś logarytm naturalny jest ściśle rosnący).
Potem oczywiście rozważasz to jako funkcję \(\displaystyle{ \lambda}\), różniczkujesz etc.

mac18
Użytkownik
Użytkownik
Posty: 316
Rejestracja: 3 wrz 2013, o 19:09
Płeć: Mężczyzna
Lokalizacja: Warszawa
Podziękował: 88 razy
Pomógł: 23 razy

Re: Metoda Największej Wiarygodności i rozkład Poissona

Post autor: mac18 » 31 paź 2017, o 18:03

Ok, policzyłem to, sprawdziłem czy jest nieobciążony i zgodny (oba tak) i muszę jeszcze zobaczyć czy jest najefektywniejszy. Chcę go porównać z estymatorem, który wyjdzie z Rao-Cramera.

Liczę logarytm z funkcji gęstości a potem pochodną po lambda i dochodzę do:
\(\displaystyle{ \frac{ \partial \ln f(k,\lambda)}{ \partial \lambda} = -1 + \frac{k}{\lambda}}\)

Muszę policzyć teraz z tego wartość oczekiwaną do kwadratu i tu zaczynają się schody, bo nie miałem jeszcze tego na rachunku.
\(\displaystyle{ E[-1+ \frac{k}{\lambda}]^{2}}\). Jedyne co przychodzi mi do głowy to wyciągnąć jeden przez lamdba, ale potem nie wiem.

Edit, już wiem

Awatar użytkownika
Premislav
Użytkownik
Użytkownik
Posty: 15211
Rejestracja: 17 sie 2012, o 13:12
Płeć: Mężczyzna
Lokalizacja: Warszawa
Podziękował: 161 razy
Pomógł: 5046 razy

Re: Metoda Największej Wiarygodności i rozkład Poissona

Post autor: Premislav » 31 paź 2017, o 20:38

Podsumowanie: estymator, który wychodzi (czyli średnia), jest efektywny (tj. osiąga dolne ograniczenie wariancji z nierówności Cramera-Rao), obliczeń nie będę zapisywać, bo za długo by to trwało, a nic ciekawego tam nie ma.

ODPOWIEDZ