Dane z obserwacji- szacowanie parametrów rozkładu gamma

Procesy stochastyczne. Sposoby racjonalizowania wielkich ilości informacji. Matematyka w naukach społecznych.
laser15
Użytkownik
Użytkownik
Posty: 721
Rejestracja: 13 lis 2011, o 14:58
Płeć: Mężczyzna
Lokalizacja: Kielce
Podziękował: 8 razy

Dane z obserwacji- szacowanie parametrów rozkładu gamma

Post autor: laser15 »

Witam, mam na podstawie próbki oszacować parametry rozkładu gamma metodą momentów, dane:
\(\displaystyle{ X_1=1; X_2=1,03; X_4= 1,04; X_5=1,02; X_6=1,11; X_7=0,92;
\newline X_8=1,01; X_9=0,98; X_{10}=0,95}\)


Proszę o pomoc, jak się za to zabrać?
squared
Użytkownik
Użytkownik
Posty: 1017
Rejestracja: 21 mar 2009, o 11:11
Płeć: Mężczyzna
Podziękował: 167 razy
Pomógł: 152 razy

Dane z obserwacji- szacowanie parametrów rozkładu gamma

Post autor: squared »

Rozkład Gamma ma następującą funkcję gęstości:
\(\displaystyle{ f(x)=\theta ^{k}x^{{k-1}}{\frac {\exp {\left(-x\theta \right)}}{\Gamma (k)}}\!}\),
czyli dwa parametry.

Metoda momentów poleca na rozwiązania równania:
\(\displaystyle{ M_{n} (\vec{x})=\mathbb{E}g_n(\xi)}\)

Oznaczenia:
\(\displaystyle{ M_n}\) - \(\displaystyle{ n-ty}\) moment centralny z próby.
\(\displaystyle{ g_n (\xi) = \xi^n}\)

No u Nas \(\displaystyle{ \xi}\) ma rozkład Gamma. Metoda momentów to cała grupa metod i technik, gdyż można brać dowolne \(\displaystyle{ n}\). Weźmy przykładowo \(\displaystyle{ n=1}\) i oszacujmy \(\displaystyle{ \theta}\).

\(\displaystyle{ M_{1} (\vec{x})=\mathbb{E}g_1(\xi)=\mathbb{E}(\xi)= \frac {k}{\theta } \\
x(n) = \frac{k}{\theta} \\
\theta = \frac{k}{x(n)}}\)
.

Średnią z próby, to łatwo wyznaczyć sobie, wiadomo jak to robimy. Przyjmując inne \(\displaystyle{ n}\) otrzymasz inny estymator.

Natomiast metoda w analogiczny sposób się rozszerza na dwa parametry - po prostu trzeba napisać dwa równania np.
\(\displaystyle{ M_{1} (\vec{x})=\mathbb{E}g_1(\xi)=\mathbb{E}\xi \\
M_{2} (\vec{x})=\mathbb{E}g_2(\xi)=\mathbb{E}\xi ^2}\)
.

Cała ewentualna trudność teraz sprowadza się do \(\displaystyle{ \mathbb{E}\xi ^2}\). Potem otrzymasz układ równań z dwoma niewiadomymi (nasze parametry). Rozwiązujesz i już masz estymatory.
ODPOWIEDZ