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