Přeskočit na obsah
coden

Střední hodnota, rozptyl a kovariance: kde rozdělení leží, jak je široké a co s čím souvisí

Střední hodnota jako těžiště, linearita, rozptyl a směrodatná odchylka, Čebyševova nerovnost, kovariance a korelace — a proč nulová kovariance ještě neznamená nezávislost.

14 minut Matěj Žalmánek 2026-10-15

Série Základy AI Díl 16 z 83
  1. Mapa umělé inteligence: co je co a jak to do sebe zapadá
  2. Stručná historie umělé inteligence: od Turinga po jazykové modely
  3. Data, model, trénink, inference: základní pojmy na jednom příkladu
  4. Jak číst matematický zápis: symboly, indexy, sumy a kvantifikátory
  5. Množiny a relace: stavebnice, ze které se skládají definice v AI
  6. Funkce a zobrazení: proč je každý model v AI funkce
  7. Logika a důkazy: pravidla, podle kterých se v AI uvažuje
  8. Grafy a stromy: struktura, na které stojí hledání i neuronové sítě
  9. Lineární algebra pro strojové učení: vektory, matice a proč AI hlavně násobí
  10. Derivace, gradient a řetízkové pravidlo: jak model zjistí, kterým směrem se učit
  11. Gradientní sestup: jak se model učí
  12. Co je pravděpodobnost: od hodu kostkou ke Kolmogorovovým axiomům
  13. Podmíněná pravděpodobnost a Bayesova věta: jak změnit odhad, když se dozvíte něco nového
  14. Náhodná veličina a distribuční funkce: jedna funkce, která o rozdělení řekne všechno
  15. Funkce hustoty a pravděpodobnostní funkce: proč má spojitá veličina v bodě pravděpodobnost nula
  16. Střední hodnota, rozptyl a kovariance: kde rozdělení leží, jak je široké a co s čím souvisí
  17. Přehled pravděpodobnostních rozdělení: Bernoulli, binomické, Poisson, rovnoměrné, normální a exponenciální
  18. Zákon velkých čísel a centrální limitní věta: proč se průměry chovají tak předvídatelně
  19. Statistika a odhady parametrů: maximální věrohodnost, MAP a intervaly spolehlivosti
  20. Teorie informace: entropie, křížová entropie a KL divergence — odkud se berou ztrátové funkce
  21. Prohledávání stavového prostoru: BFS, DFS, UCS a A* s animacemi
  22. Hry a protivník: minimax, alfa-beta ořezávání a Monte Carlo tree search
  23. Splňování omezujících podmínek: backtracking, propagace omezení a AC-3 na sudoku a mapě krajů
  24. Logika a znalostní systémy: dopředné a zpětné řetězení, rezoluce a expertní systémy
  25. Plánování: STRIPS, svět kostek, heuristiky z relaxace a plánování v AI agentech
  26. Lokální hledání a simulované žíhání: obchodní cestující přes krajská města
  27. Genetické algoritmy: populace, křížení, mutace a výběr na batohu na výlet
  28. Evoluční strategie a genetické programování: od pravidla jedné pětiny po Keplerův zákon
  29. Optimalizace mravenčí kolonií: jak feromonová stopa najde krátkou cestu
  30. Optimalizace hejnem částic: setrvačnost, vlastní zkušenost a zkušenost hejna
  31. Neuroevoluce: evoluce vah i tvaru neuronových sítí a algoritmus NEAT
  32. Fuzzy množiny a funkce příslušnosti: jak zapsat „trochu teplo“
  33. Fuzzy logika a pravděpodobnost: proč „trochu teplo“ neznamená „pravděpodobně teplo“
  34. Fuzzy logika: t-normy, implikace a pravidla KDYŽ–PAK
  35. Fuzzy řízení: Mamdani a Sugeno od ventilátoru po kyvadlo na vozíku
  36. Fuzzy shlukování: algoritmus fuzzy c-means a body mezi shluky
  37. Neuro-fuzzy systémy (ANFIS): fuzzy pravidla, která se učí z dat
  38. Markovovy řetězce: přechodová matice, stacionární rozdělení a PageRank
  39. Skryté Markovovy modely: Viterbiho algoritmus a nepoctivé kasino
  40. Bayesovské sítě: proč je tráva mokrá a co z toho plyne pro postřikovač
  41. Kalmanův filtr: jak z nepřesných měření odhadnout, kde auto je a jak rychle jede
  42. Druhy strojového učení: s učitelem, bez učitele, s částečným učitelem a další
  43. Přeučení a regularizace: proč model, který zná data nazpaměť, selže
  44. Vyhodnocení modelu: křížová validace, F1 a křivka ROC
  45. Shlukování a snížení dimenze: k-means, PCA a t-SNE
  46. Klasifikace a regrese: dvě úlohy učení s učitelem
  47. Lineární regrese: nejmenší čtverce geometricky i maticově
  48. Polynomiální a nelineární regrese: jak přímku ohnout
  49. Logistická regrese: sigmoida, šance a křížová entropie
  50. Vícetřídní klasifikace a softmax: od dvou tříd k mnoha
  51. Naivní Bayesův klasifikátor: spamový filtr, který funguje lépe, než by měl
  52. Ztrátové funkce: co přesně model minimalizuje
  53. Rozhodovací stromy a k nejbližších sousedů: modely bez vzorců
  54. Ansámbly: náhodný les a gradient boosting
  55. SVM a jádrový trik: nejširší hranice a jak ji ohnout
  56. Signál a vzorkování: proč platí Nyquistův teorém
  57. Fourierova transformace: rozklad signálu na frekvence
  58. Filtry a konvoluce: vyhlazení, hrany a odstranění šumu
  59. Zvuk a řeč: spektrogram, melová škála a MFCC
  60. Obraz jako signál: pixely, barvy, konvoluce a hrany
  61. Perceptron a umělý neuron: nejjednodušší neuronová síť a co neumí
  62. Neuronová síť od nuly: zpětné šíření chyby krok za krokem
  63. Aktivace, inicializace a normalizace: drobnosti, na kterých hluboká síť stojí
  64. Konvoluční sítě: jak síť poznává obrázky
  65. Rekurentní sítě: paměť pro posloupnosti a mizející gradient
  66. LSTM a GRU: brány, které rozhodují, co si síť zapamatuje
  67. Autoenkodéry: síť, která se učí data stlačit a zase rozbalit
  68. Embeddingy a word2vec: jak se slovo stane vektorem
  69. Mechanismus pozornosti: jak se model naučí, kam se dívat
  70. Transformery a velké jazykové modely: architektura krok za krokem
  71. Tokenizace: jak se text rozseká na kousky pro model
  72. Předtrénování a škálování: proč větší model a víc dat funguje
  73. Doladění a RLHF: jak se z doplňovače textu stane asistent
  74. Prompty, kontext a RAG: jak modelu dát správné informace
  75. Variační autoenkodéry (VAE): nové obrázky ze spojitého latentního prostoru
  76. Generativní adversariální sítě (GAN): padělatel proti detektivovi
  77. Difúzní modely: jak se z čistého šumu stane obrázek
  78. Markovovy rozhodovací procesy: stavy, akce, odměny a Bellmanovy rovnice
  79. Q-learning a DQN: učení hodnoty akcí od tabulky po neuronovou síť
  80. Policy gradient a PPO: učení přímo strategie
  81. Jak funguje AI agent: jazykový model ve smyčce s nástroji
  82. Programování s AI agenty: zadání, testy a kontrola
  83. Vyhodnocení a bezpečnost AI systémů: jak měřit a jak hlídat

Celá osnova série 83 z 83 témat hotovo

Co si odnesete
  • Co je střední hodnota, proč je to těžiště rozdělení a proč se střední hodnoty vždy sčítají.
  • Co měří rozptyl a směrodatná odchylka, kdy se rozptyly sčítají a co zaručuje Čebyševova nerovnost.
  • Co je kovariance a korelace, jak je vidět v datech a proč nulová kovariance neznamená nezávislost.

Co se hodí znát předem: Funkce hustoty a pravděpodobnostní funkce: proč má spojitá veličina v bodě pravděpodobnost nula

Rozdělení náhodné veličiny říká všechno — ale „všechno“ je často moc. Chceme-li dvě rozdělení rychle porovnat, stačí pár čísel: kde rozdělení leží, jak je široké a jestli spolu dvě veličiny rostou, nebo jdou proti sobě. Těm číslům se říká střední hodnota, rozptyl a kovariance a potkáte je v každém modelu strojového učení — od ztrátové funkce po optimalizátor Adam.

Navazuje na díly Náhodná veličina a Funkce hustoty.

1. Střední hodnota

Definice — střední hodnota diskrétní veličiny

Nechť X je diskrétní náhodná veličina, H = {x1, x2, …} ⊆ ℝ nejvýš spočetná množina navzájem různých hodnot s P(X ∈ H) = 1 a p : ℝ → [0, 1] pravděpodobnostní funkce veličiny X (viz díly Náhodná veličina a Funkce hustoty). Řekneme, že X má konečnou střední hodnotu, právě když ∑i |xi| · p(xi) < ∞ (řada konverguje absolutně). Pak je střední hodnota veličiny X číslo

E[X] = ∑i xi · p(xi).

Absolutní konvergence zaručuje, že součet nezávisí na pořadí hodnot xi. Značí se také μ nebo μX. Střední hodnota je tedy průměr hodnot vážený jejich pravděpodobnostmi.

Laicky řečeno

Vážený průměr možných hodnot, kde vahou je pravděpodobnost. U jedné kostky 1 · 1/6 + 2 · 1/6 + … + 6 · 1/6 = 1 + 2 + … + 6/6 = 3,5 (barvy podle vzorce níže). Střední hodnota přitom nemusí být hodnota, která může padnout — 3,5 na kostce nepadne nikdy. Je to „průměr, ke kterému se dlouhodobě dopracujete“.

E[X] = ∑i xi · p(xi)

E[X]
střední hodnota — jedno číslo, které má zastoupit celé rozdělení: kde „v průměru“ leží
∑i
přes všechny možné hodnoty — žádnou nevynecháme, každá k průměru přispěje
xi
hodnota — to, co se průměruje
p(xi)
váha hodnoty — častější hodnoty mají táhnout průměr víc; váhy dají dohromady 1, takže nemusíme na konci nic dělit
Při N pokusech padne xi zhruba N · p(xi)-krát; obyčejný průměr výsledků je pak ∑ xi · Np(xi)/N — přesně tenhle vzorec.

Fyzikálně je střední hodnota těžiště. Postavte pravděpodobnosti jako závaží na páku: podpěra ji udrží ve vodorovné poloze jen ve střední hodnotě. Na obrázku je větší ze dvou čísel na kostkách — maximum. Velká čísla jsou častější (šestka je maximem v 11 z 36 případů, jednička jen v jednom), takže těžiště leží vpravo od středu: E[max] = 161/36 ≈ 4,47.

1/36 1 3/36 2 5/36 3 7/36 4 9/36 5 11/36 6 E ≈ 4,47 maximum dvou kostek (oka) pravděpodobnost P (0–1)
Střední hodnota je těžiště: páka se vyrovná jen v bodě 4,47
Příklad — střední hodnota maxima dvou kostek a rovnováha páky

Maximum je nejvýš k právě tehdy, když jsou obě kostky nejvýš k: to je k² dvojic. Přesně k je tedy v k² − (k − 1)² = 2k − 1 dvojicích.

p(1), …, p(6) = 1/36, 3/36, 5/36, 7/36, 9/36, 11/36

E[max] = 1 · 1/36 + 2 · 3/36 + 3 · 5/36 + 4 · 7/36 + 5 · 9/36 + 6 · 11/36 = 1 + 6 + 15 + 28 + 45 + 66/36 = 161/36 ≈ 4,47

(barvy podle vzorce výše)

Páka: moment kolem bodu c je ∑ (k − c) · p(k) = E[max] − c, a ten je nulový právě pro c = E[max].

Definice — střední hodnota spojité veličiny

Nechť X je náhodná veličina s hustotou f : ℝ → [0, ∞). Řekneme, že X má konečnou střední hodnotu, právě když ∫−∞∞ |x| · f(x) dx < ∞ (integrál konverguje absolutně). Pak je střední hodnota veličiny X číslo

E[X] = ∫−∞∞ x · f(x) dx.

Je to tentýž vážený průměr jako v diskrétním případě, jen součet přes hodnoty nahradí integrál a pravděpodobnost p(xi) „pravděpodobnost na úseku“ f(x) dx.

Laicky řečeno

Totéž, jen se místo sčítání přes body integruje přes plochu pod hustotou. Pro čekání na e-mail z minulého dílu (hustota 0,5 · e−0,5x) vyjde E[X] = 2 minuty — přesně ten průměr, který jsme do příkladu vložili.

Příklad — střední doba čekání na e-mail

Integrál součinu spočítáme metodou per partes, ∫ u v′ = [u v] − ∫ u′ v, s u = x a v′ = λ e−λx, tedy v = −e−λx:

E[X] = ∫0∞ x · λ e−λx dx = [−x e−λx]0∞ + ∫0∞ e−λx dx

= (0 − 0) + [−e−λx/λ]0∞ = 0 + 1/λ = 1/λ

λ = 0,5   →   E[X] = 1/0,5 = 2 minuty

Člen x e−λx jde pro x → ∞ k nule, protože exponenciála klesá rychleji, než x roste.

2. Linearita střední hodnoty

Věta — linearita střední hodnoty

Nechť (Ω, 𝓕, P) je pravděpodobnostní prostor, X, Y : Ω → ℝ náhodné veličiny s konečnou střední hodnotou a a, b ∈ ℝ. Pak veličiny aX + b a X + Y (definované po výsledcích: (X + Y)(ω) = X(ω) + Y(ω)) mají konečnou střední hodnotu a platí

E[aX + b] = a · E[X] + b,    E[X + Y] = E[X] + E[Y],

a to i když X a Y nejsou nezávislé.

Důkaz (pro konečný prostor výsledků Ω s 𝓕 = 2Ω, jako jsou dvě kostky). Označme P(ω) = P({ω}). Střední hodnotu lze sčítat přes výsledky místo přes hodnoty: E[X] = ∑ω ∈ Ω X(ω) P(ω) (výsledky se stejnou hodnotou X(ω) = x dají dohromady x · p(x)). Pak

E[X + Y] = ∑ω (X(ω) + Y(ω)) P(ω) = ∑ω X(ω) P(ω) + ∑ω Y(ω) P(ω) = E[X] + E[Y].

Nikde jsme nepotřebovali, jak spolu X a Y souvisí — obě se vyhodnotí na stejném výsledku ω. Rovnost pro aX + b plyne stejně z ∑ P(ω) = 1. Obecný případ se dokazuje stejnou úvahou s integrály. ∎

Laicky řečeno

Průměr součtu je součet průměrů. Součet dvou kostek má střední hodnotu 3,5 + 3,5 = 7 — bez jakékoli tabulky. A když hodíte jen jednou a výsledek zdvojnásobíte, vyjde 2 · 3,5 = 7 taky. Linearita je nejužitečnější vlastnost v celé pravděpodobnosti: složité veličiny rozložíte na jednoduché kousky a průměry sečtete.

Příklad — kolik šestek padne na dvou kostkách

Pro jev A ∈ 𝓕 je indikátor jevu A náhodná veličina 𝟙A : Ω → {0, 1}, 𝟙A(ω) = 1 pro ω ∈ A a 𝟙A(ω) = 0 pro ω ∉ A. Počet šestek N rozložíme na dva indikátory: I1 = 𝟙A1 pro A1 = „šestka na první kostce“ a I2 = 𝟙A2 pro A2 = „šestka na druhé kostce“. Pak N = I1 + I2.

E[I1] = 1 · 1/6 + 0 · 5/6 = 1/6,   E[I2] = 1/6   (barvy podle vzorce střední hodnoty výše)

E[N] = E[I1] + E[I2] = 1/6 + 1/6 = 1/3

Kontrola výčtem: dvě šestky v 1 z 36 případů, jedna v 10 z 36   →   2 · 1 + 1 · 10/36 = 12/36 = 1/3

Střední hodnota indikátoru je pravděpodobnost jeho jevu — díky tomu se linearitou dá spočítat průměrný počet čehokoli.

3. Rozptyl a směrodatná odchylka

Definice — rozptyl

Nechť X je náhodná veličina s konečnou střední hodnotou μ = E[X]. Pokud má náhodná veličina (X − μ)² konečnou střední hodnotu (to nastane, právě když E[X²] < ∞), je rozptyl veličiny X číslo

Var(X) = E[(X − μ)²]

a říkáme, že X má konečný rozptyl. Pro diskrétní X je E[(X − μ)²] = ∑i (xi − μ)² · p(xi), pro X s hustotou f je to ∫−∞∞ (x − μ)² · f(x) dx. Rozptyl je tedy střední hodnota čtverce odchylky od střední hodnoty; vždy Var(X) ≥ 0.

Laicky řečeno

Průměrná čtvercová vzdálenost od průměru. Když hodnoty padají těsně kolem střední hodnoty, je rozptyl malý; když lítají daleko, velký. Čtverec zajistí, že odchylky nahoru a dolů se nevyruší, a velké odchylky potrestá víc než malé.

Var(X) = E[(X − μ)²]

Var(X)
rozptyl — jedno číslo pro šířku rozdělení
X − μ
odchylka od středu — šířku měříme od střední hodnoty, takže posunutí celého rozdělení rozptyl nezmění
²
čtverec — bez něj by se odchylky nahoru a dolů vyrušily a vyšla by vždy nula (E[X − μ] = μ − μ = 0); navíc velkou odchylku započte víc než dvě poloviční
E
průměr přes celé rozdělení — časté odchylky váží víc než vzácné, stejně jako u střední hodnoty
Proč čtverec, a ne absolutní hodnota? Čtverec se dá derivovat a rozepsat, a díky tomu platí pravidla níže — rozptyly nezávislých veličin se sčítají.
Věta — výpočetní tvar rozptylu

Nechť X je náhodná veličina s konečným rozptylem. Pak

Var(X) = E[X²] − (E[X])².

Důkaz. Označme μ = E[X]. Roznásobíme a použijeme linearitu, μ je konstanta: E[(X − μ)²] = E[X² − 2μX + μ²] = E[X²] − 2μ · E[X] + μ² = E[X²] − 2μ² + μ² = E[X²] − μ². ∎

Laicky řečeno

„Průměr čtverců minus čtverec průměru.“ U jedné kostky je E[X²] = 1 + 4 + 9 + 16 + 25 + 36/6 = 91/6, takže Var(X) = 91/6 − 3,5² = 35/12 ≈ 2,92. Plyne to z linearity: stačí roznásobit (X − μ)².

Definice — směrodatná odchylka

Nechť X je náhodná veličina s konečným rozptylem Var(X) ≥ 0. Směrodatná odchylka veličiny X je nezáporné číslo

σ = σX = √Var(X).

Platí tedy σ ≥ 0 a σ² = Var(X); σ má stejné jednotky jako X.

Laicky řečeno

Rozptyl má divné jednotky — u čekání v minutách vyjde v minutách na druhou. Odmocnina to vrátí zpět: σ je „typická vzdálenost od průměru“ ve stejných jednotkách jako data. U jedné kostky σ ≈ 1,71, u čekání na e-mail je rozptyl 4 a σ = 2 minuty.

Z dat místo z rozdělení

Rozdělení často neznáme, máme jen naměřená čísla x1, …, xn. Pak každé z nich dostane váhu 1/n — jako by to bylo rozdělení, ve kterém je každé měření stejně pravděpodobné:

def prumer_a_rozptyl(data):
    n = len(data)
    mu = sum(data) / n                          # E[X] s vahami 1/n
    var = sum((x - mu) ** 2 for x in data) / n  # E[(X − μ)²] s vahami 1/n
    return mu, var, var ** 0.5                  # μ, rozptyl, σ
Příklad — průměr a rozptyl osmi čísel

data 2, 4, 4, 4, 5, 5, 7, 9;   μ = 40/8 = 5

odchylky xi − μ: −3, −1, −1, −1, 0, 0, 2, 4;   čtverce 9, 1, 1, 1, 0, 0, 4, 16

Var = 9 + 1 + 1 + 1 + 0 + 0 + 4 + 16/8 = 32/8 = 4,   σ = √4 = 2

(barvy podle rozboru rozptylu výše: průměr E s vahami 1/8 je dělení osmi)

Když data jsou jen vzorek z většího celku a odhadujeme rozptyl celku, dělí se často n − 1 místo n (tzv. výběrový rozptyl, tady 32/7 ≈ 4,57): odchylky měříme od průměru vzorku, který k datům sedí líp než skutečná μ, a dělení n by rozptyl v průměru podhodnotilo.

4. Rozptyl součtu: proč se vyplatí hodit dvakrát

Věta — rozptyl lineární funkce a součtu

Nechť X, Y jsou náhodné veličiny na témž pravděpodobnostním prostoru s konečným rozptylem a a, b ∈ ℝ. Pak

  1. Var(aX + b) = a² · Var(X);
  2. jsou-li X a Y nezávislé (definice v části Kovariance a korelace), pak Var(X + Y) = Var(X) + Var(Y).

Důkaz první části. Z linearity je E[aX + b] = aμ + b, takže odchylka od středu je (aX + b) − (aμ + b) = a(X − μ) — b se odečte. Pak Var(aX + b) = E[a²(X − μ)²] = a² Var(X). Druhou část dokážeme v kapitole o kovarianci. ∎

Laicky řečeno

Posunutí o b šířku nezmění; natažení a-krát ji zvětší a²-krát. A rozptyly nezávislých veličin se sčítají — na rozdíl od směrodatných odchylek.

Porovnejme dvě hry se stejnou střední hodnotou 7. V první hodíte dvěma kostkami a sečtete je: Var = 35/12 + 35/12 = 35/6 ≈ 5,83, σ ≈ 2,42. Ve druhé hodíte jednou a výsledek zdvojnásobíte: Var = 2² · 35/12 = 35/3 ≈ 11,67, σ ≈ 3,42. Průměr je stejný, ale druhá hra je o dost divočejší — u dvou kostek se smůla na jedné často vyrovná štěstím na druhé.

součet dvou kostek 24681012 μ = 7, σ ≈ 2,42 hodnota (oka) pravděpodobnost P (0–1) jedna kostka × 2 24681012 μ = 7, σ ≈ 3,42 hodnota (oka) pravděpodobnost P (0–1)
Stejný průměr, různý rozptyl: zelená úsečka je μ ± σ

Na stejném principu stojí diverzifikace investic i to, proč průměr z více měření kolísá méně než jedno měření. Kolik přesně méně, ukáže příští díl o zákonu velkých čísel.

5. Čebyševova nerovnost

Věta — Čebyševova nerovnost

Nechť X je náhodná veličina s konečným rozptylem, μ = E[X] a σ = √Var(X) > 0. Pak pro každé k ∈ ℝ, k > 0, platí

P(|X − μ| ≥ kσ) ≤ 1/k²,

kde {|X − μ| ≥ kσ} = {ω ∈ Ω : |X(ω) − μ| ≥ kσ}.

Laicky řečeno

Daleko od průměru moc pravděpodobnosti být nemůže — a platí to pro každé rozdělení. Víc než dvě směrodatné odchylky od průměru padne nejvýš čtvrtina hodnot, víc než tři nejvýš devítina. U součtu dvou kostek je 2σ ≈ 4,83 a tak daleko od sedmičky jsou jen součty 2 a 12: skutečně 2/36 ≈ 0,056, hluboko pod zaručenou mezí 1/2² = 0,25 (barvy podle vzorce níže).

P(|X − μ| ≥ kσ) ≤ 1/k²

P
pravděpodobnost „ulétlé“ hodnoty — tu chceme shora omezit, aniž bychom znali tvar rozdělení
|X − μ|
vzdálenost od střední hodnoty — absolutní hodnota, protože nás zajímají úlety na obě strany
kσ
práh měřený ve směrodatných odchylkách — tak platí jedna nerovnost pro všechna rozdělení a všechny jednotky; „daleko“ znamená daleko vzhledem k typické vzdálenosti
1/k²
mez — čtverec pochází z rozptylu (průměrného čtverce odchylky): každá hodnota aspoň kσ od středu přispěje do rozptylu aspoň k²σ² krát svou pravděpodobnost; kdyby jich dohromady bylo víc než 1/k², vyšel by rozptyl větší než σ²
Pro k ≤ 1 je mez ≥ 1 a nic neříká; užitečná je až pro k > 1.
Věta — Markovova nerovnost a důkaz Čebyševovy

Markovova nerovnost: Nechť Y je náhodná veličina s konečnou střední hodnotou, která je nezáporná (Y(ω) ≥ 0 pro každé ω ∈ Ω), a a ∈ ℝ, a > 0. Pak

P(Y ≥ a) ≤ E[Y]/a.

Důkaz. Vždy platí Y ≥ a · I, kde I = 𝟙{Y ≥ a} je indikátor jevu {Y ≥ a} (1, když Y ≥ a, jinak 0) (pro Y ≥ a je to Y ≥ a, jinak Y ≥ 0). Střední hodnota nerovnost zachová: E[Y] ≥ a · E[I] = a · P(Y ≥ a).

Čebyšev. Použijeme Markova na Y = (X − μ)² ≥ 0 a a = k²σ² > 0:

P(|X − μ| ≥ kσ) = P((X − μ)² ≥ k²σ²) ≤ E[(X − μ)²]/k²σ² = σ²/k²σ² = 1/k². ∎

Jak volná mez je, ukazuje součet dvou kostek (σ = √35/6 ≈ 2,42) pro několik k (barvy podle vzorce výše):

kpráh kσsoučty tak daleko od 7skutečná Pmez 1/k²
12,422, 3, 4, 10, 11, 1212/36 ≈ 0,3331
1,53,622, 3, 11, 126/36 ≈ 0,1670,444
24,832, 122/36 ≈ 0,0560,25
37,25žádné00,111

U spojitého čekání na e-mail (μ = σ = 2 minuty) je pro k = 2 práh 4 minuty; dolů tak daleko čekání nedosáhne, nahoru znamená čekat aspoň 6 minut: P(X ≥ 6) = e−3 ≈ 0,050, opět hluboko pod 0,25. Mez je hrubá, protože musí platit i pro nejhorší možné rozdělení.

6. Kovariance a korelace

Definice — kovariance

Nechť X, Y : Ω → ℝ jsou náhodné veličiny na témž pravděpodobnostním prostoru (Ω, 𝓕, P) s konečným rozptylem a μX = E[X], μY = E[Y]. Kovariance veličin X a Y je číslo

Cov(X, Y) = E[(X − μX)(Y − μY)].

Tato střední hodnota je konečná, protože |(X − μX)(Y − μY)| ≤ (X − μX)² + (Y − μY)²/2. Kovariance je tedy střední hodnota součinu odchylek obou veličin od jejich středních hodnot.

Laicky řečeno

Když je X nad svým průměrem, bývá nad průměrem i Y? Pak je součin odchylek většinou kladný a kovariance kladná. Když jdou proti sobě, je záporná. Na obrázku je každý součin obdélník od průsečíku průměrů k bodu: zelené kladné, oranžové záporné. Kovariance je jejich průměrná plocha se znaménkem.

Definice — korelační koeficient

Nechť X, Y jsou náhodné veličiny na témž pravděpodobnostním prostoru s konečným rozptylem a se směrodatnými odchylkami σX > 0, σY > 0. Korelační koeficient (korelace) veličin X a Y je číslo

ρ(X, Y) = Cov(X, Y)/σX σY.

Vždy platí −1 ≤ ρ(X, Y) ≤ 1 (plyne z Cauchyovy–Schwarzovy nerovnosti |Cov(X, Y)| ≤ σX σY).

Laicky řečeno

Kovariance závisí na jednotkách — v centimetrech vyjde 10 000× větší než v metrech. Korelace ji vydělí oběma směrodatnými odchylkami a dostane číslo bez jednotek: 1 znamená, že body leží na rostoucí přímce, −1 na klesající, 0 žádnou lineární souvislost.

první kostka a součet Cov = 2,92 ρ = 0,71 první kostka (oka) součet (oka) první a druhá kostka Cov = 0,00 ρ = 0,00 první kostka (oka) druhá kostka (oka) první kostka a druhá − první Cov = −2,92 ρ = −0,71 první kostka (oka) druhá − první (oka)
Kovariance jako průměr obdélníků: kladná, nulová, záporná

První kostka a součet spolu rostou: Cov = 35/12 ≈ 2,92 a ρ ≈ 0,71. První a druhá kostka jsou nezávislé a obdélníky se přesně vyruší: Cov = 0. A čím víc padne na první kostce, tím menší je rozdíl „druhá minus první“: Cov ≈ −2,92, ρ ≈ −0,71.

Věta — výpočetní tvar kovariance

Nechť X, Y jsou náhodné veličiny na témž pravděpodobnostním prostoru s konečným rozptylem. Pak

Cov(X, Y) = E[XY] − E[X] · E[Y].

Důkaz. Stejně jako u rozptylu roznásobíme: E[XY − μYX − μXY + μXμY] = E[XY] − μYμX − μXμY + μXμY. Pro Y = X je to výpočetní tvar rozptylu: Cov(X, X) = Var(X). ∎

Příklad — první kostka a součet

X = první kostka, Y = druhá, S = X + Y;   E[XS] = E[X²] + E[XY] = 91/6 + 3,5 · 3,5 = 91/6 + 49/4 = 329/12

Cov(X, S) = 329/12 − 3,5 · 7 = 329/12 − 294/12 = 35/12 ≈ 2,92

ρ = 35/12/√35/12 · √35/6 = √35/12/√35/6 = √1/2 ≈ 0,71

E[XY] = E[X] · E[Y] jsme směli použít, protože kostky jsou nezávislé (věta níže).

Definice — nezávislé náhodné veličiny

Nechť X, Y : Ω → ℝ jsou náhodné veličiny na témž pravděpodobnostním prostoru (Ω, 𝓕, P) s distribučními funkcemi FX, FY. X a Y jsou nezávislé, právě když

∀x, y ∈ ℝ :   P(X ≤ x, Y ≤ y) = FX(x) · FY(y),

kde {X ≤ x, Y ≤ y} = {X ≤ x} ∩ {Y ≤ y}; jevy {X ≤ x} a {Y ≤ y} jsou tedy nezávislé pro každou dvojici x, y (nezávislé jevy viz díl Podmíněná pravděpodobnost). Pro diskrétní X a Y je to ekvivalentní s podmínkou

∀x, y ∈ ℝ :   P(X = x, Y = y) = P(X = x) · P(Y = y).

Laicky řečeno

Když se dozvíte hodnotu jedné veličiny, nic se nezmění na šancích té druhé. Pro dvě kostky: P(první 6, druhá 6) = 1/36 = 1/6 · 1/6. Pro součet a rozdíl (první − druhá) to neplatí: součet 12 a zároveň rozdíl 1 padnout nemůže, pravděpodobnost je 0, ale P(součet 12) · P(rozdíl 1) = 1/36 · 5/36 ≠ 0.

Věta — rozptyl součtu obecně

Nechť X, Y jsou náhodné veličiny na témž pravděpodobnostním prostoru s konečným rozptylem. Pak

Var(X + Y) = Var(X) + Var(Y) + 2 · Cov(X, Y).

Jsou-li navíc X a Y nezávislé, je E[XY] = E[X] · E[Y], tedy Cov(X, Y) = 0 a Var(X + Y) = Var(X) + Var(Y).

Důkaz. Odchylka součtu od středu je (X − μX) + (Y − μY); její čtverec je (X − μX)² + (Y − μY)² + 2(X − μX)(Y − μY) a střední hodnota z linearity dá první rovnost. U nezávislých diskrétních veličin je E[XY] = ∑x,y xy · P(X = x) P(Y = y) = (∑x x P(X = x)) · (∑y y P(Y = y)) = E[X] · E[Y], takže Cov = 0; pro veličiny s hustotou se postupuje stejně s integrály. ∎

Var(X + Y) = Var(X) + Var(Y) + 2 · Cov(X, Y)

Var(X + Y)
jak moc kolísá součet — třeba celé portfolio nebo součet dvou chyb měření
Var(X) + Var(Y)
kolísání každé části zvlášť — kdyby spolu nesouvisely, byl by to celý výsledek
Cov(X, Y)
souhra — když jdou obě veličiny nahoru zároveň, odchylky se sčítají a součet kolísá víc (kladná); když proti sobě, odchylky se ruší (záporná)
2
smíšený člen 2ab z (a + b)² — dvojice (X, Y) a (Y, X) přispějí stejně
Součet dvou kostek: 35/12 + 35/12 + 2 · 0 = 35/6. Jedna kostka × 2 je „součet kostky se sebou“: 35/12 + 35/12 + 2 · 35/12 = 35/3.
Laicky řečeno

Veličiny, které spolu rostou, rozptyl součtu zvětší; veličiny, které jdou proti sobě, ho zmenší — přesně to dělá zajištění v investicích. U nezávislých veličin člen s kovariancí zmizí a zbude pravidlo z předchozí kapitoly.

Pozor, obráceně to neplatí: nulová kovariance neznamená nezávislost. Rozdíl kostek D = první − druhá a jejich součet S mají kovarianci Var(první) − Var(druhá) = 0. Přesto jsou závislé: když padne součet 12, musí být rozdíl 0. Kovariance měří jen lineární souvislost.

Věta — kovariance je bilineární

Nechť X, Y, Z jsou náhodné veličiny na témž pravděpodobnostním prostoru s konečným rozptylem a a, b ∈ ℝ. Pak

  1. Cov(X, Y) = Cov(Y, X) (symetrie);
  2. Cov(X, X) = Var(X);
  3. Cov(aX + bY, Z) = a · Cov(X, Z) + b · Cov(Y, Z) (linearita v prvním argumentu; se symetrií i ve druhém).

Plyne to z linearity střední hodnoty, stejně jako u rozptylu.

Příklad — proč mají rozdíl a součet kovarianci nula

Cov(X − Y, X + Y) = Cov(X, X) + Cov(X, Y) − Cov(Y, X) − Cov(Y, Y)

= Var(X) − Var(Y) = 35/12 − 35/12 = 0

Smíšené členy se odečtou díky symetrii a zbude rozdíl rozptylů — nula, protože obě kostky jsou stejné. Na závislosti to nic nemění.

7. Střední hodnota a rozptyl v umělé inteligenci

  • Ztrátová funkce. Střední kvadratická chyba (MSE) je střední hodnota čtverce chyby, E[(Y − ŷ)²]. Nejlepší konstantní odhad je právě střední hodnota Y a chyba, která zbude, je rozptyl.
  • Standardizace vstupů. Než se data pošlou do sítě, převedou se na z = x − μ/σ — střední hodnota 0, směrodatná odchylka 1. Učení pak běží rychleji a stabilněji; totéž uvnitř sítě dělá batch normalization.
  • Kovarianční matice a PCA. Kovariance všech dvojic příznaků tvoří matici a její vlastní vektory (viz díl o lineární algebře) ukazují směry, ve kterých data kolísají nejvíc.
  • Optimalizátor Adam (z dílu o gradientním sestupu) průběžně odhaduje střední hodnotu gradientu a střední hodnotu jeho čtverce — a krok dělí odmocninou z toho druhého. Kde gradient hodně kolísá, dělá opatrnější kroky.
  • Vychýlení a rozptyl modelu. Chybu modelu lze rozložit na systematickou chybu (vychýlení) a kolísání mezi trénovacími sadami (rozptyl). Jednoduché modely mají velké vychýlení, složité velký rozptyl.

Proč je nejlepší konstantní odhad střední hodnota

Věta — rozklad střední kvadratické chyby

Nechť Y je náhodná veličina s konečným rozptylem a c ∈ ℝ libovolné číslo. Pak

E[(Y − c)²] = Var(Y) + (E[Y] − c)².

Funkce c ↦ E[(Y − c)²] má proto jediné minimum v bodě c = E[Y] a jeho hodnota je Var(Y).

Důkaz. Var(Y − c) = Var(Y) (posunutí šířku nemění) a výpočetní tvar rozptylu pro Y − c dává E[(Y − c)²] = Var(Y − c) + (E[Y − c])² = Var(Y) + (E[Y] − c)². ∎

E[(Y − c)²] = Var(Y) + (E[Y] − c)²

E[(Y − c)²]
střední kvadratická chyba (MSE) odhadu c — ztráta, kterou model minimalizuje
Var(Y)
neodstranitelná chyba — v c vůbec není, takže ji žádný konstantní odhad nezmenší; je to šum samotných dat
(E[Y] − c)²
vychýlení odhadu na druhou — jediné, co volbou c ovlivníme; nulové je právě pro c = E[Y]
Odtud tvrzení z prvního bodu: MSE je nejmenší pro c = E[Y] a zbude rozptyl. Stejný rozklad v obecnější podobě je rozklad na vychýlení a rozptyl modelu.
Příklad — odhad hodu kostkou a standardizace

Odhad c = 3:   E[(Y − 3)²] = 4 + 1 + 0 + 1 + 4 + 9/6 = 19/6 ≈ 3,17

Podle rozkladu: 35/12 + (3,5 − 3)² = 35/12 + 1/4 = 38/12 = 19/6 ✓;   pro c = 3,5 jen 35/12 + 0 ≈ 2,92   (barvy podle vzorce výše)

Standardizace šestky: z = 6 − 3,5/1,71 ≈ 1,46 — šestka leží 1,46 směrodatné odchylky nad průměrem

8. Shrnutí

  • Střední hodnota je vážený průměr — těžiště rozdělení. Je lineární, i pro závislé veličiny.
  • Rozptyl je průměrná čtvercová odchylka od průměru, směrodatná odchylka jeho odmocnina ve stejných jednotkách jako data.
  • Rozptyly nezávislých veličin se sčítají; proto součet dvou kostek kolísá méně než jedna kostka × 2.
  • Čebyševova nerovnost: dál než kσ od průměru leží nejvýš 1/k² pravděpodobnosti.
  • Kovariance a korelace měří lineární souvislost dvou veličin; nulová kovariance ještě neznamená nezávislost.

Máte na to jiný názor? Napište na support@coden-app.com.

Agent na vašem počítači, ovládaný z telefonu