Statistica e stima dei parametri: massima verosimiglianza, MAP e intervalli di confidenza
Campione, stimatore, correttezza e perché la varianza campionaria divide per n − 1, il metodo della massima verosimiglianza, la stima bayesiana MAP e gli intervalli di confidenza — e perché addestrare una rete neurale significa stimare parametri.
Serie Basi di IA Parte 19 di 52
- Mappa dell’intelligenza artificiale: che cos’è ciascuna cosa e come si incastra tutto
- Breve storia dell’intelligenza artificiale: da Turing ai modelli linguistici
- Dati, modello, addestramento, inferenza: i concetti di base su un unico esempio
- Come leggere la notazione matematica: simboli, indici, sommatorie e quantificatori
- Insiemi e relazioni: i mattoncini con cui si costruiscono le definizioni nell’IA
- Funzioni e applicazioni: perché ogni modello di IA è una funzione
- Logica e dimostrazioni: le regole con cui si ragiona nell’IA
- Grafi e alberi: la struttura su cui poggiano la ricerca e le reti neurali
- Algebra lineare per l’apprendimento automatico: vettori, matrici e perché l’IA soprattutto moltiplica
- Derivate, gradiente e regola della catena: come un modello capisce in quale direzione imparare
- Discesa del gradiente: come impara un modello
- Che cos’è la probabilità: dal lancio di un dado agli assiomi di Kolmogorov
- Probabilità condizionata e teorema di Bayes: come cambiare una stima quando scopri qualcosa di nuovo
- Variabile aleatoria e funzione di ripartizione: una sola funzione che dice tutto sulla distribuzione
- Funzione di densità e funzione di massa: perché una variabile continua ha probabilità zero in un punto
- Valore atteso, varianza e covarianza: dove si trova una distribuzione, quanto è larga e che cosa è legato a che cosa
- Panoramica delle distribuzioni di probabilità: Bernoulli, binomiale, Poisson, uniforme, normale ed esponenziale
- Legge dei grandi numeri e teorema del limite centrale: perché le medie si comportano in modo così prevedibile
- Statistica e stima dei parametri: massima verosimiglianza, MAP e intervalli di confidenza
- Teoria dell’informazione: entropia, entropia incrociata e divergenza KL — da dove vengono le funzioni di perdita
- Ricerca nello spazio degli stati: BFS, DFS, UCS e A* con animazioni
- Giochi e avversari: minimax, potatura alfa-beta e Monte Carlo tree search
- Soddisfacimento di vincoli: backtracking, propagazione dei vincoli e AC-3 su un sudoku e una mappa delle regioni ceche
- Logica e sistemi basati sulla conoscenza: concatenazione in avanti e all’indietro, risoluzione e sistemi esperti
- Pianificazione: STRIPS, il mondo dei blocchi, euristiche dal rilassamento e pianificazione negli agenti IA
- Ricerca locale e ricottura simulata: un commesso viaggiatore tra i capoluoghi di regione cechi
- Algoritmi genetici: popolazione, incrocio, mutazione e selezione su uno zaino per l’escursione
- Strategie evolutive e programmazione genetica: dalla regola di un quinto alla legge di Keplero
- Ottimizzazione con colonie di formiche: come una traccia di feromone trova un percorso breve
- Ottimizzazione a sciame di particelle: inerzia, esperienza propria ed esperienza dello sciame
- Neuroevoluzione: evoluzione dei pesi e della forma delle reti neurali e l’algoritmo NEAT
- Insiemi fuzzy e funzioni di appartenenza: come scrivere «un po’ caldo»
- Logica fuzzy e probabilità: perché «un po’ caldo» non significa «probabilmente caldo»
- Logica fuzzy: t-norme, implicazioni e regole SE–ALLORA
- Controllo fuzzy: Mamdani e Sugeno, dal ventilatore al pendolo inverso su un carrello
- Clustering fuzzy: l’algoritmo fuzzy c-means e i punti tra i cluster
- Sistemi neuro-fuzzy (ANFIS): regole fuzzy che imparano dai dati
- Catene di Markov: matrice di transizione, distribuzione stazionaria e PageRank
- Modelli di Markov nascosti: l’algoritmo di Viterbi e il casinò disonesto
- Reti bayesiane: perché l’erba è bagnata e che cosa significa per l’irrigatore
- Filtro di Kalman: come stimare da misure imprecise dove si trova un’auto e a che velocità va
- Tipi di apprendimento automatico: supervisionato, non supervisionato, semi-supervisionato e altri
- Prompt, contesto e RAG: come dare a un modello le informazioni giuste
- Autoencoder variazionali (VAE): immagini nuove da uno spazio latente continuo
- Reti generative avversarie (GAN): il falsario contro il detective
- Modelli di diffusione: come il rumore puro diventa un’immagine
- Processi decisionali di Markov: stati, azioni, ricompense ed equazioni di Bellman
- Q-learning e DQN: imparare il valore delle azioni, dalla tabella alla rete neurale
- Policy gradient e PPO: imparare direttamente la politica
- Come funziona un agente IA: un modello linguistico in un ciclo con strumenti
- Programmare con agenti IA: incarico, test e controllo
- Valutazione e sicurezza dei sistemi di IA: come misurare e come sorvegliare
Programma completo della serie 52 di 83 argomenti pubblicati
- Che cos’è uno stimatore, quando è corretto e consistente e perché la varianza campionaria divide per n − 1.
- Come funzionano la massima verosimiglianza e la stima MAP — e perché sono l’addestramento e la regolarizzazione.
- Che cosa significa esattamente un intervallo di confidenza al 95 % e qual è l’errore comune nel leggerlo.
Utile sapere prima: Legge dei grandi numeri e teorema del limite centrale: perché le medie si comportano in modo così prevedibile
La teoria della probabilità va dal modello ai dati: se sai che la moneta è equa, calcoli che cosa uscirà. La statistica fa il percorso inverso — dai dati al modello: sono uscite 6 croci su 20, com’è fatta questa moneta? Ed è esattamente ciò che fa l’apprendimento automatico. Il modello ha dei parametri, i dati sono ciò che è uscito e l’addestramento è una stima dei parametri. Questa parte mostra che cos’è un buono stimatore, da dove viene il metodo della massima verosimiglianza, su cui poggia quasi ogni funzione di perdita, e come dire onestamente quanto possiamo fidarci di una stima.
Riprende la parte Legge dei grandi numeri e teorema del limite centrale.
1. Campione, statistica, stimatore
Definizione — campione casuale
Un campione casuale di ampiezza n da una distribuzione con densità (o funzione di probabilità) f(x; θ) è formato da n variabili indipendenti e identicamente distribuite X1, …, Xn con questa distribuzione. Il parametro θ è ignoto.
I dati a nostra disposizione li consideriamo estratti a caso da una certa distribuzione. Il tipo di distribuzione lo scegliamo noi (per esempio una Bernoulli — spam sì o no), ma il suo parametro (quale parte della posta è spam) non lo conosciamo e vogliamo ricavarlo dai dati.
Definizione — statistica e stimatore
Una statistica è una qualsiasi funzione del campione T = T(X1, …, Xn) che non dipende dal parametro ignoto θ. La statistica che usiamo per stimare θ si chiama stimatore e si indica con θ̂.
Una statistica è un qualsiasi calcolo fatto sui dati: la media, il massimo, la mediana, la quota di spam. Uno stimatore è una statistica con un compito preciso: centrare il parametro ignoto. Siccome i dati sono casuali, è casuale anche lo stimatore: con altri dati verrebbe diverso. Per questo ci chiediamo come si comportano gli stimatori in media e man mano che i dati aumentano.
2. Un buono stimatore: corretto e consistente
Definizione — stimatore corretto
Uno stimatore θ̂ è corretto (o non distorto) se E[θ̂] = θ per ogni valore di θ. La differenza E[θ̂] − θ si chiama distorsione (bias).
Uno stimatore corretto nelle singole prove sbaglia ora da una parte ora dall’altra, ma in media centra il bersaglio. La media campionaria è uno stimatore corretto del valore atteso — l’abbiamo visto nella parte precedente.
Definizione — stimatore consistente
Uno stimatore θ̂n è consistente se per ogni ε > 0 vale limn→∞ P(|θ̂n − θ| ≥ ε) = 0.
Con abbastanza dati, lo stimatore centra il bersaglio con tutta la precisione che vuoi. La media campionaria è consistente — è esattamente la legge dei grandi numeri.
La correttezza non è scontata. Per stimare la varianza viene spontaneo usare lo scarto quadratico medio dalla media, (1/n) ∑ (Xi − X̄)². Ma questo stimatore sottostima sistematicamente la varianza: misuriamo gli scarti da X̄, che si è adattato ai dati, invece che dal vero μ.
Teorema — varianza campionaria corretta
Per un campione casuale con varianza σ² vale E[(1/n) ∑ (Xi − X̄)²] = ((n − 1)/n) σ². Per questo lo stimatore corretto della varianza è la varianza campionaria corretta
s² = 1/(n − 1) · ∑i=1n (Xi − X̄)².
Bisogna dividere per n − 1, non per n (correzione di Bessel). Un «grado di libertà» se ne va per stimare la media. Nei campioni piccoli la differenza si sente — con cinque valori, dividere per cinque sbaglia di un quinto —, con un milione di valori non conta nulla.
Verifichiamolo: abbiamo lanciato 2 000 volte cinque dadi e ogni volta abbiamo calcolato entrambe le stime della varianza. La varianza vera di un dado è 35/12 ≈ 2,92. La media delle stime con divisione per n si è fermata a 2,38 (punta a 4/5 · 2,92 ≈ 2,33), quella con divisione per n − 1 a 2,97:
3. Metodo della massima verosimiglianza
Definizione — funzione di verosimiglianza
Per i dati osservati x1, …, xn, la funzione di verosimiglianza è la funzione del parametro
L(θ) = ∏i=1n f(xi; θ).
I dati sono fissi; facciamo variare il parametro e ci chiediamo: «se il parametro fosse θ, quanto sarebbero probabili proprio questi dati?». Per osservazioni indipendenti le probabilità si moltiplicano. Attenzione: L(θ) non è una probabilità del parametro — sommata su θ non dà necessariamente 1.
Definizione — stimatore di massima verosimiglianza
Lo stimatore di massima verosimiglianza (in inglese MLE) è θ̂ = argmaxθ L(θ). Poiché il logaritmo è crescente, è la stessa cosa di argmaxθ ∑i ln f(xi; θ).
Scegli il parametro con cui i dati osservati sono più probabili. Il logaritmo trasforma il prodotto in una somma — è più comodo da calcolare e il computer non si ritrova con numeri piccolissimi arrotondati a zero.
Su 20 e-mail, 6 erano spam. Per una quota di spam p si ha L(p) = p6(1 − p)14. La derivata del logaritmo, 6/p − 14/(1 − p), si annulla per p̂ = 6/20 = 0,3 — esattamente la quota di spam nei dati. In generale, con k successi su n si ottiene p̂ = k/n. Con 60 spam su 200 la stima è la stessa, ma il picco è molto più stretto: i dati dicono di più sul parametro.
Per l’attesa di un’e-mail con distribuzione esponenziale, il metodo della massima verosimiglianza dà λ̂ = 1 / X̄. Da cinque attese misurate di 1,2; 3,5; 0,4; 2,1 e 2,8 minuti la media è di 2 minuti e λ̂ = 0,5 al minuto.
4. Stima bayesiana MAP
Definizione — stimatore del massimo a posteriori
Se π(θ) è la densità a priori del parametro, lo stimatore del massimo a posteriori (MAP) è θ̂ = argmaxθ L(θ) · π(θ).
È il teorema di Bayes della parte sulla probabilità condizionata applicato ai parametri: a posteriori ∝ verosimiglianza × a priori. Il MAP sceglie il parametro più probabile dopo aver tenuto conto di ciò che pensavamo in partenza. Se in anticipo diamo poco peso agli estremi (densità a priori proporzionale a p(1 − p)), la quota di spam risulta 7/22 ≈ 0,318 — la stima viene tirata un po’ verso il centro. Con l’aumentare dei dati, l’influenza della distribuzione a priori svanisce.
5. Intervallo di confidenza
Definizione — intervallo di confidenza
L’intervallo (T1, T2) con statistiche T1 ≤ T2 è un intervallo di confidenza per θ di livello 1 − α se P(T1 < θ < T2) = 1 − α per ogni θ.
Invece di un solo numero dai un intervallo e, insieme, una garanzia: il procedimento con cui lo calcoli centra il parametro vero nel 95 % dei casi. A essere casuali sono gli estremi dell’intervallo, non il parametro.
Per il valore atteso e un n abbastanza grande, il teorema del limite centrale fornisce un intervallo approssimato al 95 %: X̄ ± 1,96 · s/√n. Abbiamo lanciato quaranta volte 40 dadi e da ogni campione abbiamo calcolato un intervallo. La linea verticale è il valore atteso vero, 3,5:
Il valore vero è contenuto in 38 intervalli su 40 — il 95 %. Su 20 000 ripetizioni è stato il 94,1 %, un po’ sotto il 95 % promesso, perché σ lo stimiamo dai dati e n = 40 non è ancora l’infinito.
Un errore frequente: «con probabilità del 95 % μ sta nell’intervallo (3,2; 3,9)». In un intervallo concreto, μ o c’è o non c’è. Quel 95 % è una proprietà del procedimento — quanti intervalli, su molte ripetizioni, centrano il bersaglio. Un’affermazione sulla probabilità del parametro può darla solo l’approccio bayesiano.
6. La stima nell’intelligenza artificiale
- L’addestramento è una stima di massima verosimiglianza. L’entropia incrociata nella classificazione è il logaritmo negativo della verosimiglianza diviso per il numero di esempi, e l’errore quadratico medio è la stessa cosa per dati con rumore normale. Minimizzare la perdita = massimizzare la verosimiglianza.
- La regolarizzazione è MAP. La penalità sui pesi λ‖w‖² (L2, weight decay) corrisponde a una distribuzione a priori normale dei pesi, la L1 a una di Laplace. L’idea a priori «i pesi devono essere piccoli» frena l’overfitting.
- Bias e varianza del modello. Un modello è uno stimatore; i modelli complessi hanno poco bias ma molta varianza (overfitting), quelli semplici il contrario. La regolarizzazione accetta un po’ più di bias in cambio di meno varianza.
- Insieme di test e intervallo di confidenza. L’accuratezza sui dati di test è una stima dell’accuratezza reale; senza un intervallo di confidenza non puoi sapere se la differenza tra due modelli è reale. La batch normalization stima media e varianza da ogni batch, e i modelli linguistici sono addestrati con la massima verosimiglianza: massimizzano la verosimiglianza della parola successiva.
7. Riepilogo
- La statistica va dai dati al modello: uno stimatore è un calcolo sui dati che deve centrare un parametro ignoto.
- Un buono stimatore è corretto (in media mira giusto) e consistente (con più dati diventa più preciso); la varianza campionaria corretta divide per n − 1.
- Il metodo della massima verosimiglianza sceglie il parametro con cui i dati sono più probabili; per una quota è k/n.
- Il MAP aggiunge un’idea a priori; nell’apprendimento automatico è la regolarizzazione.
- Un intervallo di confidenza al 95 % nasce da un procedimento che centra il parametro vero nel 95 % dei casi.
La vedi diversamente? Scrivi a support@coden-app.com.