03 · Probabilitat i estadística
Idea central: un LLM és un model probabilístic . Entrenar-lo és estimació de màxima versemblança. Avaluar-lo és inferència estadística. Si no domines aquest bloc, confondràs soroll amb progrés.
1. Fonaments# sense fórmules
Espai mostral , esdeveniments, probabilitat amb els axiomes de Kolmogorov. El que realment necessites:
Independència: . Independència condicional: — la base dels models gràfics i de la factorització autoregressiva.
La regla de la cadena de la probabilitat (el cor d'un LLM)
Factoritzador autoregressiu
Un LLM només modela p(xₜ | x₍<ₜ₎). Com es converteix això en la probabilitat d'una frase sencera?
p(frase) =
Exemples
- tokens
- 6
- vocabulari
- 6
- log p(frase)
- -3,634
- p(frase)
- 0,026
- pèrdua mitjana
- 0,606 nats
- perplexity
- 1,83
La probabilitat de la frase és el producte dels 6 condicionals: 0,026. Sumant els logaritmes: -3,633967, i exp(-3,634) = 0,026 — diferència 3,47·10⁻¹⁸. Es treballa amb logaritmes perquè el producte de centenars de números petits arriba a zero en coma flotant de seguida.
La pèrdua d'entrenament d'un LLM és exactament = 0,6057 nats, que són 0,874 bits per token. Un model que no sabés res donaria 1,792 nats (la uniforme sobre 6 paraules); aquest en dona 0,606, és a dir que ha guanyat 1,71 bits per token.
2. Teorema de Bayes# sense fórmules
Prior, versemblança i posterior
Quantes dades calen perquè el que observes s'imposi al que ja creies?
- prior (el que creies)
- versemblança (el que has vist)
- posterior (el que hauries de creure)
Taxa observada: 70,0 %
- prior
- Beta(1, 1), mitjana 0,500
- posterior
- Beta(15,0, 7,0)
- mitjana posterior
- 0,6818
- interval creïble 95 %
- [0,478, 0,854]
- màxima versemblança
- 0,7000
- pes del prior
- 9,1 %
No hi ha cap aproximació numèrica aquí: la beta és el prior conjugat de la binomial, de manera que el posterior té forma tancada. Amb Beta(1, 1) i 14 encerts de 20, surt Beta(15,0, 7,0), amb mitjana (a+k)/(a+b+n) = 0,681818.
La mitjana posterior és exactament una mitjana ponderada del prior i de les dades: 9,1 % × 0,500 + 90,9 % × 0,700 = 0,681818. El pes del prior és (a+b)/(a+b+n): el prior equival literalment a 2 observacions imaginàries.
Per comparar: l'interval de confiança freqüentista de Wilson dona [0,481, 0,855]. Amb n gran i prior feble tots dos convergeixen, tot i que volen dir coses diferents — l'interval creïble és una afirmació sobre el paràmetre, i el de confiança, sobre el procediment.
On apareix a l'IA moderna:
- Regularització L2 ≡ prior gaussià sobre els pesos (MAP).
- Model d'un LLM com a inferència bayesiana implícita sobre "quina distribució de text estic imitant?" — una de les explicacions teòriques de l'in-context learning.
- Calibratge i quantificació d'incertesa.
logprobsd'un model es poden combinar amb un prior de tasca per a classificació.
3. Variables aleatòries i distribucions# sense fórmules
| Distribució | Paràmetres | On apareix en IA |
|---|---|---|
| Bernoulli | classificació binària, dropout | |
| Categòrica | sortida d'un LLM sobre el vocabulari | |
| Binomial | comptes d'encerts en una eval | |
| Gaussiana | inicialització, soroll, VAE, difusió | |
| Gumbel | — | Gumbel-max trick per mostrejar categòriques |
| Laplace | — | prior L1, esparsitat |
| Beta / Dirichlet | — | priors conjugats, mescles de dades |
| Poisson | comptes d'esdeveniments rars | |
| Llei de potència (Zipf) | freqüència de paraules, mides de fitxers, scaling laws |
La distribució de Zipf del llenguatge (la freqüència de la -èssima paraula més comuna ) explica per què la tokenització subparaula funciona, per què la cua llarga és tan cara i per què sempre hi ha coneixement rar que el model no ha vist prou.
Esperança, variància, covariància
Propietats que fas servir constantment:
- Linealitat: (sempre, fins i tot amb dependència).
- ; si independents, .
Aplicació directa: per què l'atenció divideix per . Si independents, aleshores té variància . Sense escalar, els logits creixen amb , la softmax se satura i els gradients s'esvaeixen. Dividir per restaura variància 1. Aquest és l'únic misteri del "scaled" dot-product attention.
Per què l'atenció divideix per √dₖ
Què li passa al softmax quan els logits creixen amb la dimensió?
- logits observats
- N(0, 1,02) teòrica
Entropia 2,376 nats de 2,773 possibles (86 %). L'atenció reparteix el pes entre diverses claus, que és el que li permet aprendre.
Els models reals fan servir entre 64 i 128 per cap.
Desactiva-ho i puja la dimensió: aquest és tot el misteri del «scaled» dot-product attention.
La mateixa llavor dona sempre el mateix resultat.
- variància dels logits
- 65,10
- desviació mostrada
- 1,009
- entropia del softmax
- 2,376 nats
- en bits
- 3,428
- pes màxim
- 21,9 %
- max p(1−p)
- 0,17114
Si qᵢ, kᵢ ~ N(0,1) independents, llavors qᵀk = Σqᵢkᵢ és una suma de dₖ termes de variància 1, i per tant té variància dₖ = 64. Mesurat sobre 4 000 productes: 65,101, una desviació del 1,7 % respecte de la teoria (l'error de mostreig esperat és de 2,2 %).
Dividint per √dₖ = 8,00, la variància passa a 1,017 ≈ 1. Per això es divideix per l'arrel i no per dₖ: l'arrel és el que restaura la desviació a 1; dividir per dₖ aixafaria els logits fins a fer-los indistingibles.
4. Teoremes límit# sense fórmules
- Llei dels grans nombres: la mitjana mostral convergeix a .
- Teorema central del límit: per a gran. L'error estàndard cau com .
Simulador del teorema central del límit
Per què la mitjana de qualsevol cosa acaba sent normal? I amb quina n comença a ser-ho?
- mitjanes observades (4.000)
- N(μ, σ²/n) teòrica
- μ de la distribució base
Amb n = 1 veus la distribució base tal qual.
La mateixa llavor dona sempre el mateix resultat.
- rèpliques dibuixades
- 4.000
- mitjana empírica
- 1,0116
- μ teòrica
- 1,0000
- desviació empírica
- 0,4540
- σ/√n teòrica
- 0,4472
- asimetria
- 0,859
La normal dibuixada no s'ha ajustat a les dades: té mitjana μ = 1,0000 i desviació σ/√n = 1,000/√5 = 0,4472, calculades només a partir de la distribució base. L'histograma mesura 1,0116 i 0,4540: una desviació relativa del 1,53 %.
L'àrea de discrepància entre l'histograma i la corba val 24,9 % de l'àrea total. Amb n = 5 encara queda asimetria (0,86): l'asimetria de la distribució base és el que triga més a desaparèixer, i cau com 1/√n.
5. Estimació: MLE, MAP i bayesiana# sense fórmules
Màxima versemblança (MLE):
MAP afegeix un prior: . Un prior gaussià dona exactament el terme de weight decay L2.
Bayesiana completa: integrar sobre tots els . Inviable per a xarxes grans; s'aproxima amb ensembles, dropout en inferència (MC-dropout), o deep ensembles.
MLE contra MAP: d'on surt el weight decay
Amb poques dades, la màxima versemblança es dispara. Què hi fa un prior?
- prior
- posterior
- MLE
- valor real
Beta(6,0, 6,0): equival a 10 observacions imaginàries repartides.
La mateixa llavor dona sempre el mateix resultat.
- dades
- 8 de 10
- MLE
- 0,8000
- MAP
- 0,6500
- mitjana posterior
- 0,6364
- valor real
- 0,700
- error del MLE
- 0,1000
- error del MAP
- 0,0500
El MAP maximitza . Amb un prior Beta(6,0, 6,0), el segon terme val 5,0·[log θ + log(1−θ)]: una penalització que estira l'estimació cap a 0.5. És exactament la mateixa estructura que el weight decay, on un prior gaussià sobre els pesos dona el terme que els estira cap a zero.
Amb 10 observacions, el prior pesa 54,5 % i les dades la resta. A mesura que n creix, MAP i MLE convergeixen: el prior importa quan no tens dades, que és precisament quan el necessites.
Propietats dels estimadors
- Biaix: . Variància: dispersió entre mostres. MSE = biaix² + variància.
- Aquest compromís reapareix arreu: models grans vs petits, regularització, mida del batch.
6. El compromís biaix–variància (i com el trenca el deep learning)# sense fórmules
Error esperat = biaix² + variància + soroll irreductible.
- Alt biaix (underfitting): el model és massa simple.
- Alta variància (overfitting): el model memoritza el soroll.
La corba clàssica en forma d'U (complexitat vs error de test) falla en el règim modern: apareix el double descent. Quan el model supera el punt d'interpolació (pot ajustar perfectament l'entrenament), l'error de test torna a baixar. Moltes xarxes profundes viuen en aquest règim sobreparametritzat. Els LLMs, en canvi, solen tenir menys paràmetres que tokens d'entrenament (un 70B en veu uns 15 bilions) i s'entrenen sovint amb menys d'una època sobre les dades: no veuen cap exemple prou vegades per memoritzar-lo, i per això "overfitting" en el sentit clàssic gairebé no és el problema principal al preentrenament (sí que ho és al fine-tuning amb poques dades).
Biaix, variància i double descent
El model massa gran sobreajusta... fins que deixa de fer-ho. Passa el punt d'interpolació i mira-ho.
- error d’entrenament
- error de test
- soroll irreductible / funció real
El punt d'interpolació és el grau 11: a partir d'aquí el model pot passar per tots els punts.
La mateixa llavor dona sempre el mateix resultat.
- paràmetres
- 4
- error d’entrenament
- 0,218
- error de test
- 0,190
- soroll irreductible
- 0,063
- norma dels coeficients
- 1,498
- millor grau
- 7
Abans del pic, l'error de test toca fons al grau 7 (0,068): és el mínim de la U clàssica. Al voltant del punt d'interpolació —grau 11, on el model té tants paràmetres com dades— es dispara fins a 0,279 (grau 10), i després torna a baixar fins a 0,101 al grau 30. No arriba a millorar el mínim clàssic, però la corba en U del llibre de text no acaba a la dreta: continua.
L'error d'entrenament arriba a 0,218 amb grau 3. A partir del punt d'interpolació no passa de 4,67·10⁻²⁴: el model passa per tots els punts, i tot i així, passat el pic, l'error de test torna a millorar en afegir-hi paràmetres. Per això als LLMs, entrenats amb menys d'una època i molt sobreparametritzats, l'«overfitting» clàssic gairebé no és el problema principal al preentrenament.
Cap model pot baixar de la variància del soroll: 0,063. Si veus un error de test per sota d'aquest valor, és soroll de mostreig del conjunt de test.
7. Inferència estadística per avaluar models# sense fórmules
Aquesta secció és la que més t'estalviarà el ridícul professional.
Interval de confiança per a una proporció
Amb ítems i encerts (Wilson és més fiable que la normal per a petit o extrem):
Comparar dos models
Com que els avalues sobre els mateixos ítems, usa un test aparellat (McNemar per a encerts binaris), que té molta més potència que comparar dues proporcions independents:
on = ítems que A encerta i B falla, = a l'inrevés.
Barres d'error en avaluacions
El model A treu 71 % i el B, 68 %. És una diferència real o és soroll?
| B encerta | B falla | total | |
|---|---|---|---|
| A encerta | 250 | 34 (b) | 284 |
| A falla | 22 (c) | 94 | 116 |
| total | 272 | 128 | 400 |
Les caselles a i d — on els dos models coincideixen — no entren al test. Només compten les 56 discordances: aquesta és la raó per la qual el test aparellat té molta més potència que comparar dues proporcions independents.
Mou-lo i mira quan la diferència passa a ser detectable.
- model A
- 71,0 % [66,4 %, 75,2 %]
- model B
- 68,0 % [63,3 %, 72,4 %]
- diferència
- +3,00 punts
- χ² de McNemar
- 2,161
- valor p
- 0,1416
- potència
- 15 %
Amb p = 0,1416, un atzar com aquest passaria 14,2 vegades de cada 100 encara que els dos models fossin idèntics. Publicar "el model A guanya per 3,0 punts" amb aquestes dades seria un error.
Per detectar una diferència de 3,0 punts amb un 80 % de potència caldrien aproximadament 3.694 ítems si els models fossin independents. Ara en tens 400. El test aparellat en necessita menys — per això s'usa McNemar i no una comparació de proporcions.
Els intervals de Wilson se superposen, però això no és el criteri correcte: dos intervals poden superposar-se i la diferència aparellada ser significativa igualment, perquè el test aparellat aprofita que els dos models responen els mateixos ítems. Aquí el test exacte dona p = 0,1409.
Bootstrap (la navalla suïssa)
Quan la mètrica és complicada (BLEU, puntuació d'un jutge LLM, pass@k), remostreja:
import numpy as np
def bootstrap_ci(scores, n_boot=10_000, alpha=0.05):
scores = np.asarray(scores)
means = [np.mean(np.random.choice(scores, len(scores), replace=True))
for _ in range(n_boot)]
return np.percentile(means, [100*alpha/2, 100*(1-alpha/2)])
Bootstrap
Com poses una barra d'error a una mètrica per a la qual no tens cap fórmula?
Cada remostreig agafa 40 d'aquestes amb reemplaçament: algunes hi surten dues vegades i d'altres no hi surten. La variabilitat entre remostrejos és l'estimació de la variabilitat mostral.
- mitjanes dels remostrejos
- mitjana observada
- interval percentil 95 %
Tipus de puntuacions
Per sobre de ~2 000 l'interval ja no es mou gaire.
La mateixa llavor dona sempre el mateix resultat.
- mitjana observada
- 3,4750
- IC percentil 95 %
- [3,200, 3,750]
- amplada
- 0,5500
- error estàndard
- 0,1476
- remostrejos fets
- 3.000
Interval percentil del bootstrap: [3,2000, 3,7500], amplada 0,5500.
Interval normal analític x̄ ± 1,96·s/√n: [3,1858, 3,7642], amplada 0,5785. Difereixen un 4,9 %.
Aquesta coincidència és el motiu per confiar-hi. El bootstrap no aporta res quan tens fórmula; aporta que segueix funcionant quan la mètrica és BLEU, un jutge LLM o pass@k, on la fórmula no existeix.
Múltiples comparacions
Si compares 20 variants de prompt amb , esperes 1 fals positiu per pur atzar. Corregeix (Bonferroni, Benjamini–Hochberg) o, millor, valida la millor variant en un conjunt held-out que no hagis tocat.
Múltiples comparacions
Si proves 20 variants de prompt, quantes «guanyen» per pur atzar?
- fals positiu
- efecte real detectat
- efecte real que s'ha escapat
- nul·la certa, correctament no rebutjada
Taxa de falsos positius entre els rebutjats: Sense corregir: 100 %Bonferroni: —Benjamini–Hochberg: —
Amb 0, totes les hipòtesis nul·les són certes: tot el que surti és fals.
La mateixa llavor dona sempre el mateix resultat.
- nul·les certes
- 20
- falsos positius sense corregir
- 1
- amb Bonferroni
- 0
- amb BH
- 0
Amb 20 hipòtesis nul·les certes i α = 0,050, esperes 1,00 falsos positius per pur atzar. En aquesta simulació n'han sortit 1. Canvia la llavor unes quantes vegades: la mitjana s'hi acostarà.
Bonferroni divideix el llindar per k (0,050/20 = 0,0025) i controla la probabilitat de tenir-ne cap a ≤ α. BH és més permissiu: controla la proporció de falsos entre els rebutjats, i per això en detecta més de reals.
La solució més bona no és estadística sinó de procediment: valida la variant guanyadora en un conjunt held-out que no hagis tocat. Una sola comparació, cap correcció necessària.
Variància del propi model
Amb temperature > 0 la sortida és estocàstica. Reporta mitjana i desviació sobre execucions, i fixa llavors quan calgui reproducibilitat. Compte: en GPU, ni amb temperature=0 tens determinisme bit a bit garantit (ordre de reducció, batching dinàmic).
8. Mostreig# sense fórmules
- Mostreig de la categòrica: és el que fa un LLM a cada pas. Vegeu temperatura/top-p al fitxer 15.
- Gumbel-max trick: amb mostreja exactament de — útil per a versions diferenciables (Gumbel-softmax) en routers de MoE i models discrets.
- Importance sampling: reponderar mostres d'una distribució per estimar esperances d'una altra. És exactament el que fa la ràtio de PPO/GRPO.
- MCMC: menys central en LLMs, però clau per entendre models basats en energia i difusió.
Mostreig d'una categòrica
Temperatura, top-p, min-p: què fa cadascun a la distribució de sortida?
- abans del filtre
- després, renormalitzat
- freqüència de 2000 mostres
S'aplica sempre, abans del filtre.
La mateixa llavor dona sempre el mateix resultat.
- tokens supervivents
- 5 de 16
- entropia abans del filtre
- 1,677 nats
- entropia després
- 1,332 nats
- en bits
- 1,921
- opcions efectives
- 3,79
- p del token més probable
- 44,1 %
El truc de Gumbel-max diu que arg max(log pᵢ + gᵢ) amb g ~ Gumbel(0,1) mostreja exactament de p. Comprovat sobre 2000 mostres: la desviació màxima entre les freqüències observades i les probabilitats teòriques és 0,0135, i l'error de mostreig esperat amb aquesta mida és de l'ordre de 0,0224.
No és una curiositat: el mateix truc, amb una softmax en comptes d'un argmax (Gumbel-softmax), és el que permet fer diferenciable una tria discreta, i és el que hi ha darrere dels encaminadors dels models de mescla d'experts.
9. Comprova que ho has entès#
0 de 6 resoltes
- Escriu la factorització autoregressiva d'una frase de 5 tokens. Quantes distribucions condicionals hi ha?factorsMarge acceptat: ± 0,100
- Demostra que minimitzar l'entropia creuada equival a la màxima versemblança.
- Un model puntua 68 % i un altre 71 % en un benchmark de 400 ítems. Quin error estàndard té cada proporció, en punts percentuals?punts percentualsMarge acceptat: 6 %
- Troba una configuració on els dos models difereixin en 3 punts i la diferència sí que sigui significativa.
Objectiu:Al widget de barres d'error, deixa una diferència d'uns 3 punts entre els dos models i ajusta les discordances b i c perquè el valor p de McNemar baixi de 0,05.
Vés al widget ↑b = per defecte · c = per defecte · hitsA = per defecte · hitsB = per defecte · n = per defecte - Per què l'atenció divideix per i no per ?
- Què és el double descent i per què desmunta la intuïció clàssica de "model massa gran = overfitting"?
10. Recursos#
Blitzstein & Hwang, Introduction to Probability (i el curs Stat110 de Harvard, gratuït).
Mathematics for Machine Learning, capítol 6.
Murphy, Probabilistic Machine Learning: An Introduction — la referència moderna.
Belkin et al. (2019), Reconciling modern ML practice and the bias-variance trade-off (double descent).
Miller (2024), Adding Error Bars to Evals — per què les barres d'error són obligatòries en avaluació d'LLMs.