Poziom zaawansowany

Zmienna losowa i wartość oczekiwana

Zmienna losowa przypisuje liczbę każdemu wynikowi doświadczenia — wygraną, czas oczekiwania, liczbę usterek. Jej rozkład zapisuje się tabelą, a dwie liczby streszczają całą tabelę: wartość oczekiwana mówi, wokół czego wyniki krążą, a wariancja — jak daleko od tego uciekają.

Zanim zaczniesz

Ten temat opiera się na wcześniejszych zagadnieniach. Zanim zaczniesz, warto przerobić poniższe lekcje — dzięki nim wszystko pójdzie gładko:

Gdzie się tego używa

Konkretne sytuacje, w których liczysz dokładnie tak, jak uczy ta lekcja:

  • Los na loterii
    Los kosztuje 5 zł. Główna wygrana 1000 zł trafia się z prawdopodobieństwem 1/2000, a pociesznie 50 zł z prawdopodobieństwem 1/200. Wartość oczekiwana wygranej wynosi 1000/2000 + 50/200 = 0,50 + 0,25 = 0,75 zł, więc oczekiwany wynik gracza to 0,75 − 5 = −4,25 zł na każdym losie. Przy stu losach rocznie daje to spodziewaną stratę 425 zł, i żadna seria szczęścia tego nie zmienia — zmienia tylko to, kiedy się to ujawni.
  • Wycena składki ubezpieczeniowej
    Ubezpieczyciel szacuje, że szkoda zdarza się u 2 procent klientów rocznie, a jej średni koszt to 30 000 zł. Wartość oczekiwana wypłaty na jedną polisę wynosi 0,02 · 30 000 = 600 zł, więc składka 900 zł zostawia 300 zł na koszty i zysk. Cała aktuarialna wycena polisy jest liczeniem wartości oczekiwanej zmiennej losowej „wypłata”.
  • Ruletka europejska
    Zakład 1 zł na jeden numer wypłaca 35 zł do 1, a numerów jest 37 razem z zerem. Wartość oczekiwana wynosi 35 · 1/37 − 1 · 36/37 = −1/37, czyli około −0,027 zł na każdą postawioną złotówkę. To jest przewaga kasyna: 2,7 procenta obrotu, niezależnie od strategii obstawiania i niezależnie od tego, ile razy się gra.
  • Ryczałt za serwis sprzętu
    Serwis wie, że dana maszyna psuje się w ciągu roku z prawdopodobieństwem 0,15, a naprawa kosztuje 800 zł. Wartość oczekiwana kosztu to 0,15 · 800 = 120 zł na maszynę, więc ryczałtowa opłata 200 zł rocznie pokrywa koszt z zapasem. Klientowi opłaca się ona wtedy, gdy woli pewne 200 zł od loterii, w której z prawdopodobieństwem 15 procent zapłaci 800 zł — i to jest cała różnica między wartością oczekiwaną a ryzykiem.

Wszystkie wzory

  • Zmienna losowa

    X:ΩRX: \Omega \to \mathbb{R}

    funkcja przypisująca liczbę wynikowi doświadczenia

  • Rozkład zmiennej

    ipi=1,pi=P(X=xi)\sum_{i} p_i = 1, \qquad p_i = P(X = x_i)

    tabela wartości i ich prawdopodobieństw

  • Wartość oczekiwana

    E(X)=ixipiE(X) = \sum_{i} x_i \, p_i

    średnia ważona prawdopodobieństwami

  • Wariancja zmiennej

    σ2(X)=i(xiE(X))2pi\sigma^2(X) = \sum_{i} \big(x_i - E(X)\big)^2 p_i

    średni kwadrat odchylenia od E(X)

  • Wzór rachunkowy

    σ2(X)=E(X2)(E(X))2\sigma^2(X) = E(X^2) - \big(E(X)\big)^2

    jeden przebieg tabeli zamiast dwóch

  • Odchylenie standardowe

    σ(X)=σ2(X)\sigma(X) = \sqrt{\sigma^2(X)}

    w tej samej jednostce co X

  • Dystrybuanta

    F(x)=P(Xx)F(x) = P(X \leq x)

    niemalejąca, od 0 do 1

  • Przekształcenie liniowe

    E(aX+b)=aE(X)+b,σ2(aX+b)=a2σ2(X)E(aX + b) = a\,E(X) + b, \qquad \sigma^2(aX + b) = a^2 \sigma^2(X)

    przesunięcie nie zmienia rozrzutu

Wynikiem rzutu dwiema kostkami jest para oczek, ale w zadaniu prawie zawsze chodzi o ich sumę. Wynikiem gry jest przebieg rozgrywki, ale interesuje nas wygrana. W obu przypadkach zdarzenie tłumaczymy na liczbę — i ten przepis tłumaczenia nazywa się zmienną losową.

Zmienna losowa to funkcja

Zmienną losową nazywamy funkcję, która każdemu zdarzeniu elementarnemu przypisuje liczbę:

X:ΩRX: \Omega \to \mathbb{R}

Nazwa myli podwójnie. To nie jest zmienna w sensie algebraicznym — nie stoi w równaniu i nie ma się jej wyznaczać. To też nie jest liczba losowa — funkcja jest całkowicie zdeterminowana, losowy jest tylko jej argument.

Dla rzutu dwiema kostkami Ω\Omega ma 3636 elementów, a XX = suma oczek przypisuje parze (3,4)(3, 4) liczbę 77. Zmienna „upraszcza" doświadczenie: zamiast 3636 par mamy 1111 liczb — ale za cenę tego, że nie są one już jednakowo prawdopodobne.

Rozkład jako tabela

Rozkładem zmiennej dyskretnej nazywamy listę jej wartości wraz z prawdopodobieństwami. Zapisuje się go tabelą:

xix_i00551010
pip_i14\tfrac{1}{4}12\tfrac{1}{2}14\tfrac{1}{4}

Warunek poprawności jest jeden i zawsze warto go sprawdzić:

pi0,ipi=1p_i \geq 0, \qquad \sum_{i} p_i = 1

Tabela wyczerpuje wszystkie możliwości, więc prawdopodobieństwa muszą się sumować dokładnie do jedynki. Ten sam rozkład rysuje się słupkami:

−10123456789101100,20,40,6wygrana (zł)p1/41/21/4gra A
Rozkład wygranej w grze A. Wysokości słupków sumują się do jedynki — to ta sama tabela, tylko narysowana.

Wartość oczekiwana

Wartością oczekiwaną zmiennej losowej nazywamy średnią jej wartości ważoną prawdopodobieństwami:

E(X)=ixipiE(X) = \sum_{i} x_i \, p_i

Fizycznie jest to środek ciężkości rozkładu: punkt, w którym rysunek ze słupkami dałby się podeprzeć palcem. Dla gry A:

E(X)=014+512+1014=0+2,5+2,5=5E(X) = 0 \cdot \frac{1}{4} + 5 \cdot \frac{1}{2} + 10 \cdot \frac{1}{4} = 0 + 2{,}5 + 2{,}5 = 5

Dwie uwagi, obie ważniejsze niż sam rachunek.

Wartość oczekiwana nie musi być możliwa. Dla rzutu kostką E(X)=1+2++66=3,5E(X) = \tfrac{1+2+\dots+6}{6} = 3{,}5, choć kostka nigdy nie pokazuje 3,53{,}5. To środek ciężkości, a nie prognoza pojedynczego wyniku.

Wartość oczekiwana to nie jest średnia z próby. Średnia arytmetyczna liczy się z tego, co wypadło, i zmienia się z próby na próbę. Wartość oczekiwana liczy się z modelu, czyli z prawdopodobieństw, i jest jedna. Po tysiącu rzutów kostką średnia będzie blisko 3,53{,}5, po dziesięciu może wyjść 2,82{,}8 — a E(X)E(X) w obu przypadkach nadal wynosi 3,53{,}5.

Los na loterii kosztuje 5 zł. Wygrana 1000 zł zdarza się z prawdopodobieństwem 1/2000, a 50 zł z prawdopodobieństwem 1/200. Czy warto grać?

Wariancja: to, czego wartość oczekiwana nie mówi

Sama E(X)E(X) streszcza rozkład zbyt mocno. Porównajmy grę A z grą B, w której wygrana wynosi 44, 55 albo 66 zł z tymi samymi prawdopodobieństwami:

−10123456789101100,20,40,6wygrana (zł)p1/41/21/4gra B
Gra B na tej samej skali co gra A. Ta sama wartość oczekiwana 5 zł, ale słupki stoją ciasno przy środku zamiast rozjeżdżać się na całą oś.

Obie gry mają E(X)=5E(X) = 5, a różnią się wszystkim, co dla gracza istotne. Tę różnicę mierzy wariancja zmiennej:

σ2(X)=i(xiE(X))2pi\sigma^2(X) = \sum_{i} \big(x_i - E(X)\big)^2 p_i

W praktyce wygodniejszy jest wzór rachunkowy, wymagający jednego przejścia po tabeli zamiast dwóch:

σ2(X)=E(X2)(E(X))2\sigma^2(X) = E(X^2) - \big(E(X)\big)^2
Policz wariancję i odchylenie standardowe dla obu gier.

Uwaga na jednostki: wariancja jest w złotych do kwadratu, więc porównywalne z samą zmienną jest dopiero odchylenie standardowe σ(X)\sigma(X). To ta sama zależność, która wystąpiła przy miarach rozrzutu zbioru danych — z tą różnicą, że tam wagi były licznościami, a tu są prawdopodobieństwami.

Przydają się jeszcze dwie własności przekształceń liniowych:

E(aX+b)=aE(X)+b,σ2(aX+b)=a2σ2(X)E(aX + b) = a\,E(X) + b, \qquad \sigma^2(aX + b) = a^2 \sigma^2(X)

Przesunięcie bb nie zmienia rozrzutu — przesuwa cały rozkład, nie rozciąga go. Dlatego w drugim wzorze bb w ogóle nie występuje.

Dystrybuanta

Zamiast pytać „ile wynosi XX", często pytamy „czy XX nie przekroczy pewnej wartości". Odpowiada na to dystrybuanta:

F(x)=P(Xx)F(x) = P(X \leq x)

Dla zmiennej dyskretnej jest to funkcja schodkowa: stała między wartościami, skacząca o pip_i w każdej z nich. Dla gry A:

przedziałx<0x < 00x<50 \leq x < 55x<105 \leq x < 10x10x \geq 10
F(x)F(x)0014\tfrac{1}{4}34\tfrac{3}{4}11

Dystrybuanta jest zawsze niemalejąca, zaczyna od 00 i kończy na 11 — bo prawdopodobieństwa się dokładają, a nie odejmują. Wygodna jest tym, że prawdopodobieństwo przedziału czyta się z niej odejmowaniem: P(a<Xb)=F(b)F(a)P(a < X \leq b) = F(b) - F(a).

Wzmianka o zmiennej ciągłej

Wszystko powyżej zakładało, że wartości da się wypisać. Bywa inaczej: czas oczekiwania na autobus, wzrost losowo wybranej osoby czy błąd pomiaru przyjmują wszystkie wartości z przedziału.

Dla takiej zmiennej ciągłej prawdopodobieństwo trafienia w konkretną liczbę wynosi zero — nie dlatego, że to niemożliwe, tylko dlatego, że możliwości jest nieskończenie wiele. Sensowne jest pytanie o przedział, a odpowiada na nie pole pod wykresem funkcji gęstości:

P(aXb)=abf(x)dxP(a \leq X \leq b) = \int_a^b f(x)\,dx

Suma z tabeli zamienia się w całkę oznaczoną, a warunek „prawdopodobieństwa sumują się do jedynki" — w „pole pod całą gęstością wynosi jeden". Najważniejszą zmienną ciągłą, czyli rozkład normalny, poznasz dwie lekcje dalej.

Zadania

Zestaw pyta o dwie liczby, obie liczone z tej samej tabeli rozkładu. Prompt E(X)E(X) podaje pary „wartość → prawdopodobieństwo" i prosi o wartość oczekiwaną. Prompt σ2(X)\sigma^2(X) podaje taką samą tabelę i prosi o wariancję zmiennej — czyli o E(X2)(E(X))2E(X^2) - \big(E(X)\big)^2.

Odpowiedzią jest tu liczba, a nie ułamek nieskracalny: wartość oczekiwana nie jest prawdopodobieństwem, tylko wielkością w skali samej zmiennej. Prawdopodobieństwa w tabelach są tak dobrane, że wynik ma najwyżej dwa miejsca po przecinku.

Ćwiczenia

Rozwiąż zestaw zadań — trudność rośnie z każdym kolejnym. Na końcu zobaczysz wynik i listę pomyłek do powtórki.

Zadanie 1 z 8Punkty: 0
E(X): 0 → 3/5, 2 → 1/5, 3 → 1/5

Generator nie prosi o podanie rozkładu — odpowiedzią byłaby cała tabela wraz z nazwami zdarzeń, a nie liczba. Nie pyta też o dystrybuantę, z tego samego powodu. Oba masz rozpisane w przykładach powyżej.

Częste błędy

  • Mylenie E(X)E(X) ze średnią arytmetyczną danych — pierwsza jest własnością modelu i jest stała, druga własnością próby i zmienia się przy każdym powtórzeniu.
  • Oczekiwanie, że E(X)E(X) jest jedną z możliwych wartości3,53{,}5 oczka nie wypada nigdy, a mimo to jest poprawną wartością oczekiwaną.
  • Tabela, w której prawdopodobieństwa nie sumują się do jedynki — to nie jest rozkład i każda liczona z niego wielkość jest bez sensu. Warto zaczynać od tej kontroli.
  • Podnoszenie do kwadratu prawdopodobieństw zamiast wartości — we wzorze E(X2)E(X^2) kwadraty biorą xix_i, a wagi pip_i zostają nietknięte.
  • Porównywanie wariancji z samą zmienną — wariancja jest w jednostkach do kwadratu; porównywalne jest dopiero odchylenie standardowe.
  • Wniosek, że dwie zmienne o tej samej E(X)E(X) są równoważne — gry A i B mają tę samą wartość oczekiwaną i dwadzieścia pięć razy różną wariancję, co dla grającego jest całą różnicą.

Karta wzorów

Temat: Zmienna losowa

  • Zmienna losowa

    X:ΩRX: \Omega \to \mathbb{R}

    funkcja przypisująca liczbę wynikowi doświadczenia

  • Rozkład zmiennej

    ipi=1,pi=P(X=xi)\sum_{i} p_i = 1, \qquad p_i = P(X = x_i)

    tabela wartości i ich prawdopodobieństw

  • Wartość oczekiwana

    E(X)=ixipiE(X) = \sum_{i} x_i \, p_i

    średnia ważona prawdopodobieństwami

  • Wariancja zmiennej

    σ2(X)=i(xiE(X))2pi\sigma^2(X) = \sum_{i} \big(x_i - E(X)\big)^2 p_i

    średni kwadrat odchylenia od E(X)

  • Wzór rachunkowy

    σ2(X)=E(X2)(E(X))2\sigma^2(X) = E(X^2) - \big(E(X)\big)^2

    jeden przebieg tabeli zamiast dwóch

  • Odchylenie standardowe

    σ(X)=σ2(X)\sigma(X) = \sqrt{\sigma^2(X)}

    w tej samej jednostce co X

  • Dystrybuanta

    F(x)=P(Xx)F(x) = P(X \leq x)

    niemalejąca, od 0 do 1

  • Przekształcenie liniowe

    E(aX+b)=aE(X)+b,σ2(aX+b)=a2σ2(X)E(aX + b) = a\,E(X) + b, \qquad \sigma^2(aX + b) = a^2 \sigma^2(X)

    przesunięcie nie zmienia rozrzutu

−10123456789101100,20,40,6wygrana (zł)p1/41/21/4gra A
Rozkład wygranej w grze A: 0 zł, 5 zł albo 10 zł z prawdopodobieństwami 1/4, 1/2 i 1/4. Słupki sumują się do jedynki, a środek ciężkości tej tabeli — wartość oczekiwana — leży w 5 zł.
−10123456789101100,20,40,6wygrana (zł)p1/41/21/4gra B
Rozkład wygranej w grze B: 4 zł, 5 zł albo 6 zł z tymi samymi prawdopodobieństwami, narysowany na tej samej skali. Ta sama wartość oczekiwana 5 zł, wariancja dwadzieścia pięć razy mniejsza.

Najczęstsze pytania

Powiązane artykuły