“
La probabilità di un evento, conoscendo il verificarsi di un altro, si ricalcola con una formula semplice e
potente: \(P(B|A) = P(A \cap B)/P(A)\). È il modo matematico di dire “data l’informazione, aggiorno le mie aspettative”. Da questa semplice
formula discende quasi tutta la statistica bayesiana, che oggi pervade scienza, medicina e intelligenza
artificiale.”
— Pierre-Simon Laplace, Théorie analytique des probabilités (1812) — prima trattazione sistematica della probabilità condizionata e della formula di disintegrazione
____________________________________________________________________________________
Nel capitolo precedente (122) abbiamo introdotto, intuitivamente, l’idea di eventi dipendenti e indipendenti, e il calcolo delle probabilità composte tramite diagrammi ad albero. In questo capitolo formalizziamo:
Sono i fondamenti del calcolo delle probabilità moderno, da cui nel prossimo capitolo (124) dedurremo lo schema di Bernoulli, e poi (125) il teorema di Bayes.
Definizione 123.1 — Probabilità condizionata
Senso intuitivo. Sapere che \(A\) si è verificato restringe lo spazio campionario da \(\Omega \) ad \(A\). La nuova “misura” di \(B\) è la frazione di \(A\) che cade in \(B\): \(|A \cap B|/|A|\), che generalizzando i casi finiti diventa \(P(A \cap B)/P(A)\).
Proprietà della probabilità condizionata:
La funzione \(B \mapsto P(B|A)\) è essa stessa una probabilità sullo spazio \(A\) (verifica gli assiomi di Kolmogorov, con \(A\) come nuovo “universo”). Da qui possiamo applicare tutte le formule probabilistiche, sostituendo \(P\) con \(P(\cdot | A)\). Ad esempio:
Lancio un dado. Probabilità di ottenere “\(2\)” sapendo che è uscito un numero pari?
\(A\) = “pari” = \(\{2, 4, 6\}\), \(B\) = “è \(2\)” = \(\{2\}\).
\(P(A) = 3/6 = 1/2\). \(P(A \cap B) = P(\{2\}) = 1/6\).
\(P(B|A) = (1/6)/(1/2) = 1/3\).
Senza informazione, \(P(B) = 1/6\). Con l’informazione che è pari, sale a \(1/3\).
Teorema della probabilità composta: dati due eventi \(A, B\) con \(P(A) > 0\),
Equivalentemente, con \(P(B) > 0\): \(P(A \cap B) = P(B) \cdot P(A|B)\).
Dimostrazione. Dalla definizione di probabilità condizionata: \(P(B|A) = P(A \cap B)/P(A)\). Moltiplicando entrambi i membri per \(P(A)\): \(P(A \cap B) = P(A) P(B|A)\).
Significato. “La probabilità che \(A\) e \(B\) si verifichino entrambi è la probabilità che \(A\) si verifichi, moltiplicata per la probabilità che \(B\) si verifichi dato che \(A\) si è verificato.” È esattamente la regola del prodotto lungo un percorso del diagramma ad albero (cap. 122).
Probabilità composta per \(n\) eventi (regola della catena):
Esempio. Estraggo \(3\) carte da un mazzo da \(52\) senza ripristino. Probabilità di estrarre \(3\) assi:
Definizione 123.2 — Partizione dello spazio campionario
Esempi di partizione:
Teorema della probabilità totale (formula di disintegrazione): sia \(\{A_1, A_2, \ldots , A_n\}\) una partizione di \(\Omega \) (con tutti i \(P(A_i) > 0\)) e sia \(B\) un evento qualunque. Allora
Dimostrazione. Poiché gli \(A_i\) partizionano \(\Omega \), anche le intersezioni \(B \cap A_i\) partizionano \(B\) (cioè \(B = \bigcup _i (B \cap A_i)\) con unione disgiunta). Quindi:
usando il teorema della probabilità composta nell’ultimo passaggio.
Caso particolare: per \(n = 2\), partizione \(\{A, A^c\}\):
È la forma a due rami (più frequente nelle applicazioni).
La probabilità di un evento \(B\) è la media pesata delle probabilità condizionate \(P(B|A_i)\), con pesi pari alle probabilità \(P(A_i)\) delle “alternative” della partizione. Equivalentemente: si “disintegra” la probabilità di \(B\) in tante porzioni, una per ogni \(A_i\), e si sommano.
Una fabbrica ha due macchine: \(M_1\) produce \(60\%\) dei pezzi, \(M_2\) il \(40\%\). La frazione di pezzi difettosi è \(2\%\) per \(M_1\), \(5\%\) per \(M_2\). Probabilità che un pezzo a caso sia difettoso?
Partizione: \(A_1 = \) “proviene da \(M_1\)”, \(A_2 = \) “da \(M_2\)”. \(P(A_1) = 0{,}6\), \(P(A_2) = 0{,}4\). \(P(D|A_1) = 0{,}02\), \(P(D|A_2) = 0{,}05\).
Probabilità totale: \(P(D) = 0{,}6 \cdot 0{,}02 + 0{,}4 \cdot 0{,}05 = 0{,}012 + 0{,}020 = 0{,}032 = 3{,}2\%\).
Definizione 123.3 — Eventi indipendenti (definizione formale)
Equivalenze (per \(P(A), P(B) > 0\)): \(A, B\) indipendenti se e solo se vale una qualunque di queste tre condizioni:
Dimostrazione equivalenza 1 \(\leftrightarrow \) 2. Da \(P(A \cap B) = P(A) P(B)\) e \(P(B|A) = P(A \cap B)/P(A) = P(B)\). Viceversa, \(P(B|A) = P(B)\) implica \(P(A \cap B)/P(A) = P(B)\), quindi \(P(A \cap B) = P(A) P(B)\).
Errore comune. “Indipendenti” e “incompatibili” sono concetti diversi:
In generale, eventi con probabilità non nulla incompatibili non sono indipendenti: se sai che \(A\) si è verificato, \(B\) certamente non si verifica (\(P(B|A) = 0 \ne P(B)\)), quindi \(A\) ti dà informazione su \(B\).
Definizione 123.4 — Indipendenza di \(n\) eventi
Attenzione. “Indipendenza a coppie” (cioè \(P(A_i \cap A_j) = P(A_i) P(A_j)\) per ogni \(i, j\)) è più debole dell’indipendenza collettiva: ci possono essere eventi a coppie indipendenti ma non collettivamente indipendenti.
Esempio 123.1 — Probabilità condizionata da un albero
Un’urna contiene \(3\) palline rosse e \(7\) blu. Estraggo \(2\) palline senza ripristino. Calcolare \(P(R_2 | R_1)\), dove \(R_i\) = “\(i\)-esima estratta è rossa”.
Calcolo. \(P(R_1) = 3/10\), \(P(R_1 \cap R_2) = (3/10)(2/9) = 6/90 = 1/15\).
Condizionata: \(P(R_2 | R_1) = P(R_1 \cap R_2)/P(R_1) = (1/15)/(3/10) = 10/45 = 2/9\).
Verifica intuitiva: dopo aver estratto una rossa, restano \(2\) rosse su \(9\), quindi \(P(R_2|R_1) = 2/9\). \(\checkmark\)
Esempio 123.2 — Probabilità composta per tre eventi
Estraggo \(3\) carte da un mazzo da \(52\) senza ripristino. Calcolare la probabilità che siano tutte di cuori.
Catena: \(P(H_1 \cap H_2 \cap H_3) = P(H_1) P(H_2|H_1) P(H_3|H_1 \cap H_2)\).
Numeri: \(P(H_1) = 13/52 = 1/4\). Dopo una di cuori, restano \(12\) su \(51\): \(P(H_2|H_1) = 12/51\). Dopo due, \(11\) su \(50\): \(P(H_3| \ldots ) = 11/50\).
Prodotto: \((13/52)(12/51)(11/50) = 1716/132600 = 11/850 \approx 0{,}0129\).
Esempio 123.3 — Probabilità totale: due macchine
Riprendo l’esempio della fabbrica. \(M_1\) produce \(60\%\) con difetto \(2\%\), \(M_2\) produce \(40\%\) con difetto \(5\%\). Probabilità che un pezzo difettoso provenga da \(M_1\)?
Strategia. Cerco \(P(M_1 | D)\). Per Bayes (anticipo cap. 125):
Probabilità totale: \(P(D) = 0{,}6 \cdot 0{,}02 + 0{,}4 \cdot 0{,}05 = 0{,}032\).
Bayes: \(P(M_1 | D) = (0{,}02 \cdot 0{,}6)/0{,}032 = 0{,}012/0{,}032 = 0{,}375\).
Interpretazione. Anche se \(M_1\) produce di più, è \(M_2\) ad avere più alta frazione di difetti, quindi la maggior parte dei difetti viene da \(M_2\) (\(62{,}5\%\)).
Esempio 123.4 — Eventi indipendenti: lancio di dadi
Lancio due dadi. Eventi:
Sono indipendenti?
Calcolo. \(P(A) = 3/6 = 1/2\). \(P(B) = \) casi con somma pari/\(36\) \(= 18/36 = 1/2\).
\(P(A \cap B) = \) “primo dispari E somma pari” = “primo dispari, secondo dispari”. (Dispari + dispari = pari.) \(P(A \cap B) = (3/6)(3/6) = 9/36 = 1/4\).
Verifica: \(P(A) \cdot P(B) = (1/2)(1/2) = 1/4 = P(A \cap B)\). \(\checkmark\)Indipendenti.
Esempio 123.5 — Verificare l’indipendenza (controesempio)
Lancio una moneta tre volte. Eventi:
Verificare se \(A, B, C\) sono a coppie indipendenti e collettivamente indipendenti.
Probabilità singole: \(P(A) = P(B) = P(C) = 1/2\).
Coppie: \(P(A \cap B) = 1/4 = P(A) P(B)\). \(\checkmark\)\(P(A \cap C) = \) “T al 1° e numero pari di T”. Casi: TTT? no (3 teste, dispari); TTC? no (2 teste, pari, sì); TCT? no (2 teste, sì); TCC? no (1 testa, dispari). Quindi TTC, TCT: \(2/8 = 1/4 = P(A)P(C)\). \(\checkmark\)\(P(B \cap C)\): simmetria, \(1/4\). \(\checkmark\)
Indipendenza collettiva: \(P(A \cap B \cap C) = P(\)TTT con numero pari di T\()\). TTT ha \(3\) teste, dispari, dunque \(A \cap B \cap C\) contiene solo TT\(x\) con numero pari di T, cioè TTC. \(P = 1/8\). Ma \(P(A) P(B) P(C) = 1/8\). \(\checkmark\)
Conclusione: anche collettivamente indipendenti. (In altri esempi può non essere così!)
Esempio 123.6 — Lampadina rossa o verde
In una scatola ci sono \(10\) lampadine: \(6\) rosse e \(4\) verdi. Tra le rosse, \(2\) sono difettose; tra le verdi, \(1\) è difettosa. Estraggo una lampadina a caso. Probabilità che sia difettosa?
Partizione: \(R\) = rossa, \(V\) = verde. \(P(R) = 6/10 = 0{,}6\), \(P(V) = 4/10 = 0{,}4\).
Condizionate: \(P(D|R) = 2/6 = 1/3\), \(P(D|V) = 1/4\).
Probabilità totale: \(P(D) = 0{,}6 \cdot (1/3) + 0{,}4 \cdot (1/4) = 0{,}2 + 0{,}1 = 0{,}3\).
Esempio 123.7 — Sondaggio elettorale
Un sondaggio campione di \(1000\) persone: \(60\%\) ha votato \(X\), \(40\%\) ha votato \(Y\). Tra chi ha votato \(X\), \(70\%\) è soddisfatto. Tra chi ha votato \(Y\), \(30\%\) è soddisfatto. Probabilità che una persona qualunque del campione sia soddisfatta?
Probabilità totale: \(P(S) = 0{,}6 \cdot 0{,}7 + 0{,}4 \cdot 0{,}3 = 0{,}42 + 0{,}12 = 0{,}54 = 54\%\).
Esercizio 123.1
Sono indipendenti? Calcolare \(P(A), P(B), P(A \cap B)\) e verificare.
Verificare se sono a coppie indipendenti e collettivamente indipendenti.
Riepilogo
Probabilità condizionata: per \(P(A) > 0\),
Si legge “\(B\) dato \(A\)”. È una nuova probabilità sullo spazio \(A\).
Teorema della probabilità totale (disintegrazione): per partizione \(\{A_i\}\),
Forma più semplice (con \(\{A, A^c\}\)): \(P(B) = P(A) P(B|A) + P(A^c) P(B|A^c)\).