Preskočiť na obsah
coden

Funkcia hustoty a pravdepodobnostná funkcia: prečo má spojitá premenná v bode pravdepodobnosť nula

Pravdepodobnostná funkcia, hustota pravdepodobnosti a plocha pod krivkou — prečo hustota nie je pravdepodobnosť, môže byť väčšia ako 1 a ako ju z dát odhadne histogram.

13 minút Matěj Žalmánek 2026-10-15

Séria Základy AI Diel 15 z 52
  1. Mapa umelej inteligencie: čo je čo a ako to do seba zapadá
  2. Stručná história umelej inteligencie: od Turinga po jazykové modely
  3. Dáta, model, tréning, inferencia: základné pojmy na jednom príklade
  4. Ako čítať matematický zápis: symboly, indexy, sumy a kvantifikátory
  5. Množiny a relácie: stavebnica, z ktorej sa skladajú definície v AI
  6. Funkcie a zobrazenia: prečo je každý model v AI funkcia
  7. Logika a dôkazy: pravidlá, podľa ktorých sa v AI uvažuje
  8. Grafy a stromy: štruktúra, na ktorej stojí hľadanie aj neurónové siete
  9. Lineárna algebra pre strojové učenie: vektory, matice a prečo AI hlavne násobí
  10. Derivácie, gradient a reťazové pravidlo: ako model zistí, ktorým smerom sa učiť
  11. Gradientný zostup: ako sa model učí
  12. Čo je pravdepodobnosť: od hodu kockou ku Kolmogorovovým axiómam
  13. Podmienená pravdepodobnosť a Bayesova veta: ako zmeniť odhad, keď sa dozviete niečo nové
  14. Náhodná premenná a distribučná funkcia: jedna funkcia, ktorá o rozdelení povie všetko
  15. Funkcia hustoty a pravdepodobnostná funkcia: prečo má spojitá premenná v bode pravdepodobnosť nula
  16. Stredná hodnota, rozptyl a kovariancia: kde rozdelenie leží, aké je široké a čo s čím súvisí
  17. Prehľad pravdepodobnostných rozdelení: Bernoulliho, binomické, Poissonovo, rovnomerné, normálne a exponenciálne
  18. Zákon veľkých čísel a centrálna limitná veta: prečo sa priemery správajú tak predvídateľne
  19. Štatistika a odhady parametrov: maximálna vierohodnosť, MAP a intervaly spoľahlivosti
  20. Teória informácie: entropia, krížová entropia a KL divergencia — odkiaľ sa berú stratové funkcie
  21. Prehľadávanie stavového priestoru: BFS, DFS, UCS a A* s animáciami
  22. Hry a protivník: minimax, alfa-beta orezávanie a Monte Carlo tree search
  23. Splňovanie obmedzujúcich podmienok: backtracking, propagácia obmedzení a AC-3 na sudoku a mape krajov
  24. Logika a znalostné systémy: dopredné a spätné reťazenie, rezolúcia a expertné systémy
  25. Plánovanie: STRIPS, svet kociek, heuristiky z relaxácie a plánovanie v AI agentoch
  26. Lokálne hľadanie a simulované žíhanie: obchodný cestujúci cez krajské mestá
  27. Genetické algoritmy: populácia, kríženie, mutácia a výber na batohu na výlet
  28. Evolučné stratégie a genetické programovanie: od pravidla jednej pätiny po Keplerov zákon
  29. Optimalizácia mravčou kolóniou: ako feromónová stopa nájde krátku cestu
  30. Optimalizácia rojom častíc: zotrvačnosť, vlastná skúsenosť a skúsenosť roja
  31. Neuroevolúcia: evolúcia váh aj tvaru neurónových sietí a algoritmus NEAT
  32. Fuzzy množiny a funkcie príslušnosti: ako zapísať „trochu teplo“
  33. Fuzzy logika a pravdepodobnosť: prečo „trochu teplo“ neznamená „pravdepodobne teplo“
  34. Fuzzy logika: t-normy, implikácie a pravidlá AK–TAK
  35. Fuzzy riadenie: Mamdani a Sugeno od ventilátora po kyvadlo na vozíku
  36. Fuzzy zhlukovanie: algoritmus fuzzy c-means a body medzi zhlukmi
  37. Neuro-fuzzy systémy (ANFIS): fuzzy pravidlá, ktoré sa učia z dát
  38. Markovove reťazce: prechodová matica, stacionárne rozdelenie a PageRank
  39. Skryté Markovove modely: Viterbiho algoritmus a nepoctivé kasíno
  40. Bayesovské siete: prečo je tráva mokrá a čo z toho vyplýva pre postrekovač
  41. Kalmanov filter: ako z nepresných meraní odhadnúť, kde auto je a ako rýchlo ide
  42. Druhy strojového učenia: s učiteľom, bez učiteľa, s čiastočným učiteľom a ďalšie
  43. Prompty, kontext a RAG: ako modelu dať správne informácie
  44. Variačné autoenkodéry (VAE): nové obrázky zo spojitého latentného priestoru
  45. Generatívne adversariálne siete (GAN): falšovateľ proti detektívovi
  46. Difúzne modely: ako sa z čistého šumu stane obrázok
  47. Markovove rozhodovacie procesy: stavy, akcie, odmeny a Bellmanove rovnice
  48. Q-learning a DQN: učenie hodnoty akcií od tabuľky po neurónovú sieť
  49. Policy gradient a PPO: učenie priamo stratégie
  50. Ako funguje AI agent: jazykový model v slučke s nástrojmi
  51. Programovanie s AI agentmi: zadanie, testy a kontrola
  52. Vyhodnocovanie a bezpečnosť AI systémov: ako merať a ako strážiť

Celá osnova série Hotovo: 52 z 83 tém

Čo si odnesiete
  • Čo je pravdepodobnostná funkcia diskrétnej premennej a ako súvisí s distribučnou funkciou.
  • Čo je hustota, prečo je pravdepodobnosť intervalu plocha pod ňou a prečo hustota môže byť väčšia ako 1.
  • Ako z dát odhadnúť hustotu histogramom a na čo slúžia hustoty v umelej inteligencii.

Čo sa hodí vedieť vopred: Náhodná premenná a distribučná funkcia: jedna funkcia, ktorá o rozdelení povie všetko, Derivácie, gradient a reťazové pravidlo: ako model zistí, ktorým smerom sa učiť

Diel o náhodnej premennej sme skončili paradoxom: to, že na autobus budete čakať presne 3 minúty, má pravdepodobnosť nula — a predsa nejakú presnú dobu čakať budete. Distribučná funkcia s tým počítať vie, ale nie je veľmi názorná. Tento diel ukáže dve funkcie, ktoré sa kreslia a čítajú oveľa ľahšie: pravdepodobnostnú funkciu pre diskrétne premenné a hustotu pre spojité.

1. Pravdepodobnostná funkcia

Definícia — pravdepodobnostná funkcia

Nech (Ω, 𝓕, P) je pravdepodobnostný priestor a X : Ω → ℝ diskrétna náhodná premenná (pozri diel Náhodná premenná), t. j. existuje najviac spočítateľná množina navzájom rôznych hodnôt H = {x1, x2, …} ⊆ ℝ s P(X ∈ H) = 1. Pravdepodobnostná funkcia premennej X je funkcia pX : ℝ → [0, 1],

pX(x) = P(X = x) = P({ω ∈ Ω : X(ω) = x})   pre x ∈ ℝ.

Hodnota pX(x) je teda pravdepodobnosť, že X vyjde presne x; pre x ∉ H je nulová.

Laicky povedané

Ku každej hodnote napíše, s akou pravdepodobnosťou padne. Pre súčet dvoch kociek je p(7) = 6/36 ≈ 0,167, p(2) = 1/36 a p(7,5) = 0 — súčet 7,5 padnúť nemôže. Po anglicky sa jej hovorí probability mass function, skrátene PMF: pravdepodobnosť je „hmota“ sústredená do jednotlivých bodov.

00,10,2 2 3 4 5 6 7 8 9 10 11 12 súčet k (oká) pravdepodobnosť p(k) (0–1)
Pravdepodobnostná funkcia súčtu dvoch kociek
Veta — vlastnosti pravdepodobnostnej funkcie

Nech X je diskrétna náhodná premenná, H = {x1, x2, …} ⊆ ℝ najviac spočítateľná množina navzájom rôznych hodnôt s P(X ∈ H) = 1, p : ℝ → [0, 1] pravdepodobnostná funkcia a F : ℝ → [0, 1] distribučná funkcia premennej X. Potom

  1. p(x) ≥ 0 pre každé x ∈ ℝ;
  2. ∑i p(xi) = 1;
  3. F(x) = ∑i : xi ≤ x p(xi) pre každé x ∈ ℝ.

Body 2 a 3 vyplývajú z aditivity: javy {X = xi} sú navzájom nezlučiteľné a ich zjednotenie je {X ∈ H}, resp. {X ∈ H, X ≤ x}, ktoré sa od {X ≤ x} líši len o jav s pravdepodobnosťou 0.

Laicky povedané

Stĺpce nie sú záporné a spolu dajú jednotku. Distribučná funkcia z minulého dielu je len priebežný súčet stĺpcov zľava: F(7) = 1/36 + 2/36 + … + 6/36 = 21/36. Naopak p(x) je výška schodu, ktorý F v bode x urobí.

Príklad — pravdepodobnostná funkcia súčtu dvoch kociek jedným vzorcom

Všetkých 36 dvojíc (a, b) je rovnako možných. Súčet k dávajú dvojice (a, k − a), kde 1 ≤ a ≤ 6 a zároveň 1 ≤ k − a ≤ 6, teda max(1, k − 6) ≤ a ≤ min(6, k − 1). Takých a je

min(6, k − 1) − max(1, k − 6) + 1 = 6 − |k − 7|,   k = 2, 3, …, 12,

p(k) = 6 − |k − 7|/36:   p(2) = 6 − 5/36 = 1/36,   p(7) = 6 − 0/36 = 6/36,   p(10) = 6 − 3/36 = 3/36.

Kontrola druhej vlastnosti: 1 + 2 + 3 + 4 + 5 + 6 + 5 + 4 + 3 + 2 + 1 = 36, takže ∑ p(k) = 36/36 = 1.

Stĺpce rastú o 1/36 až k sedmičke a potom rovnako klesajú — preto je graf trojuholník.

2. Prečo spojitá premenná nemá pravdepodobnosť v bode

Vezmime spojitú premennú: dobu X v minútach, kým vám príde ďalší e-mail, keď v priemere chodí jeden za 2 minúty. Takéto čakanie má exponenciálne rozdelenie s distribučnou funkciou F(x) = 1 − e−0,5x pre x ≥ 0 (prečo práve takú, ukáže diel o rozdeleniach). Pravdepodobnostná funkcia tu nepomôže — v každom bode by bola nula. Skúsme teda namiesto bodu interval [2, 2 + h] a zužovať ho:

Definícia — exponenciálne rozdelenie

Nech λ ∈ ℝ, λ > 0, a X je náhodná premenná s distribučnou funkciou F : ℝ → [0, 1]. X má exponenciálne rozdelenie s parametrom (intenzitou) λ práve vtedy, keď pre každé x ∈ ℝ

F(x) = 1 − e−λx   pre x ≥ 0   a   F(x) = 0   pre x < 0.

Číslo e = limn→∞ (1 + 1/n)n ≈ 2,718 je základ prirodzeného logaritmu. Premenná teda nadobúda len nezáporné hodnoty a pravdepodobnosť, že prekročí x ≥ 0, je 1 − F(x) = e−λx.

Laicky povedané

λ je, koľko udalostí príde v priemere za jednotku času — pri e-mailoch jeden za 2 minúty, teda λ = 0,5 za minútu. Pravdepodobnosť, že budete čakať dlhšie než x, je 1 − F(x) = e−λx: každá ďalšia minúta ju vynásobí rovnakým číslom e−0,5 ≈ 0,61.

00,250,5 0123456 čakanie x (min) hustota f(x) (1/min)
Interval [2, 2 + h] sa zužuje
šírka hP(2 < X ≤ 2 + h)P / h
20,23250,1163
10,14470,1447
0,50,08140,1627
0,10,01790,1794
0,010,00180,1835
→ 0→ 0→ f(2) ≈ 0,1839

Pravdepodobnosť intervalu ide k nule — 0,2325, 0,1447, …, 0,0018 —, a preto v jedinom bode vyjde nula. Ale pravdepodobnosť na jednotku dĺžky, teda P/h, sa ustáli: 0,1163, 0,1447, 0,1627, 0,1794, 0,1835 … a blíži sa k 0,1839. Toto číslo už o bode 2 niečo hovorí: koľko pravdepodobnosti pripadá na minútu okolo druhej minúty. Je to hustota.

f(x) = limh → 0+ P(x < X ≤ x + h)/h

f(x)
hustota v bode x — číslo, ktoré má o jednom bode niečo povedať, hoci jeho pravdepodobnosť je nula
limh → 0+
zužovanie intervalu — kým má interval nenulovú šírku, opisuje celý úsek, nie bod x; v limite zostane len bod
P(x < X ≤ x + h)
pravdepodobnosť úzkeho intervalu za bodom — sama ide k nule, preto ju nemôžeme vziať priamo
h
šírka intervalu — delením z pravdepodobnosti urobíme pravdepodobnosť na jednotku dĺžky, a tá sa k nule neblíži
Rovnaká stavba ako derivácia: F(x + h) − F(x)/h — zmena distribučnej funkcie delená dĺžkou kroku.
Príklad — prvý riadok tabuľky a limita

P(2 < X ≤ 3) = F(3) − F(2) = (1 − e−1,5) − (1 − e−1) = e−1 − e−1,5 ≈ 0,3679 − 0,2231 = 0,1447

h = 1,   takže P/h = 0,1447/1 = 0,1447

V limite h → 0+ je to derivácia F v bode 2: (1 − e−0,5x)′ = 0,5 · e−0,5x,   f(2) = 0,5 · e−1 ≈ 0,1839

(farby podľa vzorca vyššie)

Veta — spojitá premenná má v každom bode pravdepodobnosť nula

Nech X je náhodná premenná, ktorej distribučná funkcia F : ℝ → [0, 1] je spojitá v každom bode ℝ. Potom

∀x ∈ ℝ :   P(X = x) = 0.

Dôkaz. Pre každé h > 0 je jav {X = x} časťou javu {x − h < X ≤ x}, takže 0 ≤ P(X = x) ≤ F(x) − F(x − h). Pravá strana ide pre h → 0 zo spojitosti F k nule, a preto P(X = x) = 0. ∎

3. Plocha pod krivkou: určitý integrál

Definícia — určitý integrál

Nech a, b ∈ ℝ, a < b, a f : [a, b] → ℝ je funkcia. Pre n ∈ ℕ, n ≥ 1, rozdeľme [a, b] na n dielikov [ti−1, ti], i = 1, …, n, kde

Δx = b − a/n,   ti = a + i · Δx   (i = 0, 1, …, n),

a v každom dieliku zvoľme ľubovoľný bod xi ∈ [ti−1, ti]. Funkcia f je (riemannovsky) integrovateľná na [a, b] práve vtedy, keď existuje číslo I ∈ ℝ také, že pri každej voľbe bodov xi (pre každé n zvlášť) platí

limn→∞ ∑i=1n f(xi) · Δx = I.

Číslo I je potom určitý integrál funkcie f od a do b a označuje sa ∫ab f(x) dx. Integrál je teda limita súčtov plôch obdĺžnikov, ak existuje a nezávisí od toho, kde v dielikoch výšku meriame. (Obvyklá Riemannova definícia pripúšťa aj nerovnako široké dieliky, ktorých najväčšia šírka ide k nule; dáva rovnaký pojem.) Pre x ∈ ℝ a funkciu integrovateľnú na každom [a, x] sa definuje nevlastný integrál ∫−∞x f(t) dt = lima→−∞ ∫ax f(t) dt, ak limita existuje; obdobne ∫a∞ a ∫−∞∞.

Laicky povedané

Integrál je plocha pod grafom funkcie medzi a a b. Rozkrájate ju na úzke obdĺžniky „výška krát šírka“, sčítate ich a čím užšie obdĺžniky, tým presnejší výsledok. Je to opačná operácia k derivácii: derivácia hovorí, ako rýchlo niečo pribúda, integrál sčíta, koľko toho pribudlo celkovo.

∫ab f(x) dx = limn→∞ ∑i=1n f(xi) · Δx

∫ab f(x) dx
plocha pod grafom medzi a a b — to, čo chceme; znak ∫ je natiahnuté S ako „suma“
limn→∞
čím viac dielikov, tým užšie obdĺžniky a tým menej prečnievajú cez krivku alebo pod ňou chýbajú; presná plocha je to, k čomu sa súčty blížia
∑i=1n
sčítanie všetkých dielikov — plocha celku je súčet plôch kúskov
f(xi)
výška obdĺžnika — funkcia na úzkom dieliku takmer nemení hodnotu, takže ju nahradíme jednou hodnotou vnútri dielika
Δx
šírka dielika — výška krát šírka je plocha obdĺžnika; v zápise integrálu z nej zostane dx
Integrál je limita súčtu plôch obdĺžnikov „výška krát šírka“.

Definícia je zároveň návod, ako integrál spočítať počítačom — obdĺžniková metóda berie výšku v strede každého dielika:

def integral(f, a, b, n):
    dx = (b - a) / n
    plocha = 0
    for i in range(n):
        x = a + (i + 0.5) * dx    # stred i-teho dielika
        plocha += f(x) * dx       # výška krát šírka
    return plocha
Príklad — plocha pod hustotou čakania medzi 1 a 3 štyrmi obdĺžnikmi

f(x) = 0,5 · e−0,5x,   a = 1,   b = 3,   n = 4,   Δx = 3 − 1/4 = 0,5

stredy dielikov xi = 1,25; 1,75; 2,25; 2,75   →   výšky f(xi) = 0,2676; 0,2084; 0,1623; 0,1264

(0,2676 + 0,2084 + 0,1623 + 0,1264) · 0,5 = 0,7647 · 0,5 ≈ 0,3824

s 10 dielikmi 0,3832, so 100 dielikmi 0,3834 — presná hodnota (spočítame ju nižšie) je 0,3834   (farby podľa vzorca vyššie)

Veta — Newtonova–Leibnizova formula

Nech a, b ∈ ℝ, a < b, f : [a, b] → ℝ je spojitá funkcia a F : [a, b] → ℝ je funkcia, pre ktorú F′(x) = f(x) pre každé x ∈ [a, b] (v krajných bodoch jednostranná derivácia); taká F sa nazýva primitívna funkcia k f. Potom je f integrovateľná na [a, b] a

∫ab f(x) dx = F(b) − F(a).

Laicky povedané

Keď poznáte funkciu, ktorej derivácia je f, nemusíte sčítať obdĺžniky: plocha je rozdiel jej hodnôt na koncoch. Rýchlosť auta integrovaná cez čas dá prejdenú dráhu — a tú prečítate ako rozdiel dvoch stavov tachometra.

4. Funkcia hustoty

Definícia — hustota pravdepodobnosti

Nech X je náhodná premenná s distribučnou funkciou FX : ℝ → [0, 1] a f : ℝ → [0, ∞) funkcia integrovateľná na každom ohraničenom intervale. f je hustota premennej X práve vtedy, keď

∀x ∈ ℝ :   FX(x) = ∫−∞x f(t) dt,

kde vpravo je nevlastný integrál z definície vyššie. Premenná, ktorá hustotu má, sa nazýva absolútne spojitá; jej distribučná funkcia je spojitá, takže je spojitá aj v zmysle dielu Náhodná premenná. Hustota nie je určená jednoznačne: ak zmeníme f v konečne mnohých bodoch, integrály sa nezmenia.

Laicky povedané

Hustota je krivka, pod ktorou leží celá pravdepodobnosť. Plocha pod ňou naľavo od x je pravdepodobnosť, že premenná vyjde najviac x. Kde je krivka vysoko, tam hodnoty padajú často; kde je nízko, zriedka. Pri čakaní na e-mail je najvyššie pri nule a potom klesá — krátke čakania sú časté, dlhé zriedkavé.

Veta — pravdepodobnosť ako plocha

Nech X je náhodná premenná s distribučnou funkciou F a hustotou f. Potom

  1. pre všetky a, b ∈ ℝ, a < b, platí P(a < X ≤ b) = ∫ab f(x) dx = F(b) − F(a);
  2. ∫−∞∞ f(x) dx = 1;
  3. ak je f spojitá v bode x ∈ ℝ, tak je F v bode x diferencovateľná a F′(x) = f(x).

Dôkaz. 1: F(b) − F(a) = ∫−∞b f − ∫−∞a f = ∫ab f a podľa vety o pravdepodobnosti intervalu je to P(a < X ≤ b). 2: ∫−∞∞ f = limx→∞ F(x) = 1. 3: F(x + h) − F(x)/h = 1/h ∫xx+h f(t) dt je priemer f na intervale dĺžky |h|, a ten ide zo spojitosti f v x k f(x). ∎

Laicky povedané

Pravdepodobnosť intervalu je plocha pod hustotou nad týmto intervalom a celá plocha je 1. Hustota je derivácia distribučnej funkcie — presne to, k čomu sa blížil pomer P/h v tabuľke. Pre čakanie na e-mail vyjde f(x) = 0,5 · e−0,5x.

Príklad — hustota exponenciálneho rozdelenia a plocha pod ňou

Pre x > 0:   f(x) = F′(x) = (1 − e−λx)′ = 0 − (−λ) · e−λx = λ · e−λx;   pre x < 0 je f(x) = 0.

Celková plocha: ∫0∞ λ e−λx dx = limb→∞ (F(b) − F(0)) = limb→∞ (1 − e−λb) − 0 = 1

Pre λ = 0,5:   f(0) = 0,5,   f(2) ≈ 0,1839,   f(4) ≈ 0,0677 — každé 2 minúty hustota klesne na e−1 ≈ 37 %.

To, že e-mail príde medzi prvou a treťou minútou, je plocha pod krivkou medzi 1 a 3. Podľa Newtonovej–Leibnizovej formuly ju nemusíme sčítať po obdĺžnikoch: F(3) − F(1) = (1 − e−1,5) − (1 − e−0,5) = e−0,5 − e−1,5 ≈ 0,6065 − 0,2231 ≈ 0,383.

00,250,5 012345678910 čakanie x (min) hustota f(x) (1/min)
Hustota čakania na e-mail a plocha medzi 1 a 3 minútami P(1 < X ≤ 3): 0,383

5. Hustota nie je pravdepodobnosť

Hodnota hustoty v bode nie je pravdepodobnosť — je to pravdepodobnosť na jednotku. Pravdepodobnosť dostanete, až keď ju vynásobíte šírkou intervalu: pre úzky interval platí P(x < X ≤ x + h) ≈ f(x) · h. Pre x = 2 a h = 0,1 je to 0,1839 · 0,1 ≈ 0,0184, v skutočnosti 0,0179 (farby podľa vzorca nižšie).

P(x < X ≤ x + h) ≈ f(x) · h

P(x < X ≤ x + h)
skutočná pravdepodobnosť — presne je to plocha pod hustotou nad úzkym intervalom
f(x)
hustota — pravdepodobnosť na jednotku dĺžky; na úzkom intervale sa takmer nemení, takže ju berieme ako výšku obdĺžnika
h
šírka intervalu — až vynásobením dĺžkou sa z „na jednotku“ stane pravdepodobnosť; preto môže byť f veľká a súčin napriek tomu malý
Plocha pod krivkou nahradená jedným obdĺžnikom; chyba mizne rýchlejšie než h, keď je f spojitá v x.

Preto môže byť hustota väčšia ako 1. Keď prídete na zastávku, kam autobus chodí každú pol minútu, čakanie je rovnomerné na intervale [0; 0,5] a hustota je tam všade 2 — plocha obdĺžnika 0,5 × 2 je stále 1. Je to ako hustota obyvateľstva: 2 000 obyvateľov na km² nehovorí, koľko ľudí býva v jednom dome — počet dostanete, až keď hustotu vynásobíte plochou.

Definícia — rovnomerné rozdelenie

Nech a, b ∈ ℝ, a < b, a X je náhodná premenná. X má rovnomerné rozdelenie na intervale [a, b] práve vtedy, keď má hustotu f : ℝ → [0, ∞),

f(x) = 1/b − a   pre x ∈ [a, b],    f(x) = 0   pre x ∉ [a, b].

Integráciou dostaneme F(x) = x − a/b − a pre x ∈ [a, b], 0 vľavo a 1 vpravo — tá istá distribučná funkcia ako v definícii z dielu Náhodná premenná, obe definície sa teda zhodujú.

Laicky povedané

Všetky hodnoty v intervale sú „rovnako husté“ a celá jednotka sa rozprestrie na dĺžku b − a. Pre autobus po pol minúte je f = 1/0,5 − 0 = 2; pre autobus po 10 minútach by to bolo 1/10 = 0,1. Čím kratší interval, tým vyššia hustota — plocha musí zostať 1.

6. Odhad hustoty z dát: histogram

Definícia — histogram ako odhad hustoty

Nech n ∈ ℕ, n ≥ 1, x1, …, xn ∈ ℝ sú dáta, w > 0 šírka koša a c ∈ ℝ počiatok. Koše sú intervaly

Bk = [c + kw, c + (k + 1)w),   k ∈ ℤ;

navzájom sa neprekrývajú a pokrývajú celé ℝ, takže každé x ∈ ℝ leží v práve jednom koši B(x). Pre kôš B označme nB = |{i ∈ {1, …, n} : xi ∈ B}| počet dát v ňom. Histogram (odhad hustoty) je funkcia f̂ : ℝ → [0, ∞),

f̂(x) = nB(x)/n · w   pre x ∈ ℝ.

Na každom koši je teda histogram konštantný a rovný podielu dát v koši delenému šírkou koša.

Laicky povedané

Spočítate, koľko meraní padlo do každého koša, a vydelíte to počtom meraní aj šírkou koša — potom je celková plocha stĺpcov presne 1 a histogram je odhad hustoty. Je to to isté ako pomer P/h, len namiesto skutočnej pravdepodobnosti použijete podiel dát.

f̂(x) = nB/n · w

f̂(x)
odhad hustoty v bode x — strieška hovorí „spočítané z dát“, nie skutočná f
nB
počet meraní v koši B, do ktorého x padne — čím viac dát v okolí x, tým je tam hustota vyššia
n
počet všetkých meraní — z počtu urobí podiel nB/n, odhad pravdepodobnosti koša; inak by viac dát dalo vyššie stĺpce
w
šírka koša — z pravdepodobnosti koša urobí pravdepodobnosť na jednotku dĺžky, rovnako ako h v pomere P/h; bez nej by širšie koše dávali vyššie stĺpce
Plocha stĺpca je w · nB/n · w = nB/n, a keďže každé meranie leží v práve jednom koši, dajú stĺpce spolu ∑B nB/n = n/n = 1.

Dáta k histogramu sme si vyrobili. Počítač vie losovať len rovnomerne z intervalu [0, 1); exponenciálne čakanie z toho dostaneme takýmto prevodom:

Veta — metóda inverznej transformácie

Nech λ ∈ ℝ, λ > 0, a U : Ω → [0, 1) je náhodná premenná s rovnomerným rozdelením na [0, 1] (hodnota 1 má pravdepodobnosť 0, takže ju smieme vynechať). Potom náhodná premenná

X = −ln(1 − U)/λ

má exponenciálne rozdelenie s parametrom λ.

Dôkaz. Keďže 1 − U ∈ (0, 1], je ln(1 − U) ≤ 0 a X ≥ 0, takže pre x < 0 je P(X ≤ x) = 0. Pre x ≥ 0 je −ln(1 − U)/λ ≤ x ⇔ 1 − U ≥ e−λx ⇔ U ≤ 1 − e−λx, takže P(X ≤ x) = P(U ≤ 1 − e−λx) = 1 − e−λx = F(x). ∎

Laicky povedané

Vylosované číslo u beriete ako „hladinu“ distribučnej funkcie a hľadáte, kde ju F dosiahne: x = F−1(u). Ak padne u = 0,5, je to −ln 0,5/0,5 ≈ 1,39 minúty — medián čakania. Rovnako sa dá losovať z každého rozdelenia, ktorého distribučnú funkciu vieme obrátiť.

def histogram(n, lam=0.5, w=0.5, konec=10):
    pocty = [0] * int(konec / w)
    for _ in range(n):
        u = nahodne()                    # rovnomerne z [0, 1)
        x = -log(1 - u) / lam            # inverzná transformácia
        if x < konec:
            pocty[int(x / w)] += 1       # do ktorého koša padne
    return [c / (n * w) for c in pocty]  # výšky: plocha stĺpcov je 1

Takto sme nasimulovali 100, 1 000 a 10 000 čakaní na e-mail (s košom 0,5 minúty) a cez histogramy nakreslili skutočnú hustotu:

00,250,5 012345678910 čakanie x (min) hustota (1/min)
Histogramy zo 100, 1 000 a 10 000 čakaní a skutočná hustota n = 100, najväčšia odchýlka: 0,098 n = 1 000, najväčšia odchýlka: 0,017 n = 10 000, najväčšia odchýlka: 0,008

Zo 100 čakaní sa stĺpec od skutočnosti líši až o 0,098, z 1 000 o 0,017 a z 10 000 o 0,008. Podiel čakaní do jednej minúty vyšiel 0,440, 0,397 a 0,390 — skutočná pravdepodobnosť je F(1) ≈ 0,393. S pribúdajúcimi dátami sa histogram primkne k hustote, rovnako ako sa empirická distribučná funkcia primkla k skutočnej.

Príklad — prvý stĺpec histogramu

Kôš [0; 0,5), w = 0,5. Zo 100 čakaní ich doň padlo 27:   f̂ = 27/100 · 0,5 = 0,54

Z 10 000 čakaní ich padlo 2 172:   f̂ = 2 172/10 000 · 0,5 ≈ 0,434   (farby podľa vzorca vyššie)

Histogram odhaduje priemernú hustotu na koši: F(0,5) − F(0)/0,5 = 1 − e−0,25/0,5 ≈ 0,442

Zo 100 čakaní je prvý stĺpec o 0,098 vyšší, než má byť — to je tá najväčšia odchýlka z textu vyššie.

7. Diskrétne a spojité vedľa seba

Diskrétna premennáSpojitá premenná
Opisuje jupravdepodobnostná funkcia p(x)hustota f(x)
P(X = x)p(x)0
Hodnotymedzi 0 a 1≥ 0, môžu byť aj väčšie ako 1
Spolu∑ p(xi) = 1∫ f(x) dx = 1
P(a < X ≤ b)súčet stĺpcovplocha pod krivkou
Distribučná funkciaschodisko, skoky o p(x)spojitá, sklon f(x)

8. Hustoty a pravdepodobnostné funkcie v umelej inteligencii

  • Klasifikátor aj jazykový model vracajú pravdepodobnostnú funkciu. Posledná vrstva (softmax) dá ku každej triede alebo každému slovu číslo medzi 0 a 1 a všetky spolu dajú 1 — presne vlastnosti pravdepodobnostnej funkcie.
  • Vierohodnosť. Model sa učí tak, aby trénovacie dáta mali čo najväčšiu pravdepodobnosť (pri diskrétnych dátach) alebo hustotu (pri spojitých). Záporný logaritmus tohto čísla je bežná stratová funkcia.
  • Hľadanie anomálií. Model odhadne hustotu bežnej prevádzky; transakcia alebo meranie v mieste s veľmi malou hustotou je podozrivé.
  • Generatívne modely sa učia hustotu obrázkov, zvukov či textov — a nový obrázok je vzorka z tejto naučenej hustoty.

Vierohodnosť na príklade

Definícia — vierohodnosť a záporná logaritmická vierohodnosť

Nech Θ je množina parametrov a pre každé θ ∈ Θ je fθ : ℝ → [0, ∞) hustota (alebo pravdepodobnostná funkcia) nejakého rozdelenia — model. Nech n ≥ 1 a x1, …, xn ∈ ℝ sú dáta, ktoré považujeme za hodnoty nezávislých náhodných premenných X1, …, Xn s rovnakým rozdelením. Vierohodnosť je funkcia L : Θ → [0, ∞),

L(θ) = ∏i=1n fθ(xi) = fθ(x1) · fθ(x2) · … · fθ(xn).

Pre θ s L(θ) > 0 je záporná logaritmická vierohodnosť (NLL)

NLL(θ) = −ln L(θ) = −∑i=1n ln fθ(xi).

Odhad metódou maximálnej vierohodnosti je θ̂ ∈ Θ, pre ktoré L(θ̂) ≥ L(θ) pre všetky θ ∈ Θ (ak také existuje). Keďže −ln je klesajúca, je to to isté ako NLL(θ̂) ≤ NLL(θ) pre všetky θ s L(θ) > 0 — metóda volí parametre, ktoré NLL minimalizujú.

Laicky povedané

Hustoty jednotlivých meraní sa násobia, pretože dáta sú nezávislé. Súčin mnohých malých čísel je v počítači nepoužiteľne malý, logaritmus z neho urobí súčet a mínus z maximalizácie minimalizáciu — aby to bola strata, ktorú zmenšuje gradientný zostup.

Namerali sme tri čakania na e-mail: 1, 2 a 3 minúty. Akú intenzitu λ majú? Pre exponenciálnu hustotu je ln(λ e−λx) = ln λ − λx, takže

NLL(λ) = −n · ln λ + λ · ∑i xi

NLL(λ)
strata — čím menšia, tým vyššiu hustotu dáva model nameraným dátam
−n · ln λ
výška hustoty pri nule je λ — tento člen tlačí λ nahor, pretože vyššia krivka dá dátam vyššiu hustotu
λ · ∑i xi
pokuta za dlhé čakania — s veľkým λ hustota rýchlo padá a namerané dlhšie čakania by boli nepravdepodobné; tlačí λ nadol
Najlepšie λ je kompromis oboch členov: derivácia −n/λ + ∑ xi = 0 dáva λ̂ = n/∑ xi, prevrátený priemer dát.
Príklad — strata pre tri čakania 1, 2 a 3 minúty

n = 3,   ∑ xi = 6

NLL(0,25) = −3 · ln 0,25 + 0,25 · 6 ≈ 4,159 + 1,5 = 5,659

NLL(0,5) = −3 · ln 0,5 + 0,5 · 6 ≈ 2,079 + 3 = 5,079

NLL(1) = −3 · ln 1 + 1 · 6 = 0 + 6 = 6,000   (farby podľa vzorca vyššie)

minimum: λ̂ = 3/6 = 0,5 — jeden e-mail za 2 minúty, presne priemer dát

9. Zhrnutie

  • Diskrétnu premennú opisuje pravdepodobnostná funkcia p(x) = P(X = x); stĺpce dajú spolu 1.
  • Spojitá premenná má v každom bode pravdepodobnosť 0; zmysel má pravdepodobnosť na jednotku dĺžky — hustota.
  • Hustota f je derivácia distribučnej funkcie a pravdepodobnosť intervalu je plocha pod ňou — integrál.
  • Hustota nie je pravdepodobnosť a môže byť väčšia ako 1; celková plocha pod ňou je vždy 1.
  • Histogram vydelený počtom dát a šírkou koša odhaduje hustotu z dát.

Máte na to iný názor? Napíšte na support@coden-app.com.

Agent na vašom počítači, ovládaný z telefónu