4 / 11

Teorema di Bayes e indipendenza tra eventi

Teorema di Bayes con esempi applicativi (test medico, scelta di corsi), definizione di indipendenza tra eventi multipli e regola del prodotto delle probabilità.

Ho 3 carte:

  • b/b che ha entrambi i lati bianchi
  • r/r che ha entrambi i lati rossi
  • b/r che ha un lato bianco e uno rosso

Allora pesco una carta a caso e guardo uno dei lati che è bianco. Qual è la probabilità che anche l'altro lato sia bianco?

Cioè:

  • {b/b}\{b/b\} evento in cui estraggo una carta con due lati bianchi
  • {b}\{b\} evento in cui estraggo una carta con almeno un lato bianco.

Allora voglio calcolare P(b/bb)P(b/b|b).

Notiamo che è ovvio sapere che P(bb/b)=1P(b | b/b) = 1, perché se estraggo la carta b/b allora sicuramente vedo un lato bianco. Allora:

P(b/bb)=P(b/bb)P(b)=P(b/b)bP(b)P(b/b)P(b/b)=P(b/bb)P(b/b)P(b/b)P(b)=P(bb/b)1/31/2=123=23\begin{aligned} P(b/b|b) &= \frac{P(b/b \cap b)}{P(b)} \\ &= \frac{P(b/b)\cap b}{P(b)} \cdot \frac{P(b/b)}{P(b/b)} \\ &= \frac{P(b/b\cap b)}{P(b/b)} \cdot \frac{P(b/b)}{P(b)} \\ &= P(b|b/b) \cdot \frac{1/3}{1/2} \\ &= 1 \cdot \frac{2}{3} = \frac{2}{3} \end{aligned}

Ora modifichiamo il problema: Se invece di prendere la carta da un mazzo di 3 carate, le mettessi tutte su un tavolo con un lato visibile e ne avessi scelta una di queste con il lato bianco visibile, qual è la probabilità che anche l'altro lato sia bianco?

In questo caso, avendo tutte le carte sul tavolo, ho due casi possibili per le facciate correntemente visibili:

  1. R R B
  2. R B B

Mi possono accadere entrambe con probabilità: P(1)=P(2)=1/2P(1) = P(2) = 1/2.

La domanda che voglio rispondere è la stessa: P(b/bb)P(b/b|b). Applichiamo la probabilità totale:

P(b/bb)=P(b/bb1)P(1b)+P(b/bb2)P(2b)=112+1212=34\begin{aligned} P(b/b|b) &= P(b/b|b\cap1)P(1|b) + P(b/b|b\cap2)P(2|b) \\ &= 1 \cdot \frac{1}{2} + \frac{1}{2} \cdot \frac{1}{2} = \frac{3}{4} \end{aligned}

Nel primo caso c'è solo una carta bianca visibile e le altre due rosse quindi per forza deve essere bianca da entrambi i lati Nel secondo caso, ci sono due carte bianche quindi l'altra faccia potrebbe essere in un caso rossa e l'altro bianca, quindi la probabilità è 1/2.

Teorema di Bayes

Se A e B sono eventi, con P(A)>0P(A) > 0 e P(B)>0P(B) > 0, allora:

P(AB)=P(BA)P(A)P(B)P(A|B) = \frac{P(B|A)P(A)}{P(B)}

Esempio

Immagina di sottoporti a un test medico per una malattia rara che colpisce solo 1 persona su 1000 (P(Malattia)=0.001P(Malattia) = 0.001).

Sappiamo che il test è molto sensibile: se hai la malattia, risulta positivo il 99% delle volte (P(PositivoMalattia)=0.99P(Positivo|Malattia) = 0.99). Tuttavia, il test produce anche dei falsi positivi: se sei completamente sano, risulta positivo in circa il 5% dei casi (P(PositivoSano)=0.05P(Positivo|Sano) = 0.05).

Fai il test e risulta positivo. Qual è la probabilità che tu sia effettivamente malato? In altre parole, vogliamo calcolare P(MalattiaPositivo)P(Malattia|Positivo).

L'intuito potrebbe ingannarci e farci credere che, dato il 99% di efficacia del test sui malati, siamo quasi certamente infetti. Il teorema di Bayes ci impone però di considerare la probabilità a priori (il fatto che la malattia sia rarissima).

Applicando la formula:

P(MalattiaPositivo)=P(PositivoMalattia)P(Malattia)P(Positivo)P(Malattia|Positivo) = \frac{P(Positivo|Malattia) \cdot P(Malattia)}{P(Positivo)}

Calcoliamo il denominatore P(Positivo)P(Positivo), ovvero la probabilità totale che un test risulti positivo (sia per i malati che per i sani):

P(Positivo)=(0.990.001)+(0.050.999)0.00099+0.04995=0.05094P(Positivo) = (0.99 \cdot 0.001) + (0.05 \cdot 0.999) \approx 0.00099 + 0.04995 = 0.05094

Quindi:

P(MalattiaPositivo)=0.000990.050940.0194P(Malattia|Positivo) = \frac{0.00099}{0.05094} \approx 0.0194

Sorprendentemente, pur avendo un test positivo in mano, hai meno del 2% di probabilità di essere realmente malato. Questo perché la malattia è talmente rara che la stragrande maggioranza dei test positivi sono, in realtà, falsi positivi provenienti da persone sane. Bayes ci insegna ad aggiornare le nostre convinzioni tenendo conto di quanto un evento è comune o raro di partenza.

Esercizio

Una moneta viene lanciata 2 volte, quale è la probabilità che esca testa in entrambi i lanci, sapendo che:

  1. il primo lancio è testa P(TTTx)P(TT|Tx)
  2. almeno un lancio è testa P(TTTxxT)P(TT|Tx \cup xT)

Dove il nostro Ω\Omega è formato da 4 eventi equiprobabili:

  • A1={TT}A_1 = \{TT\}
  • A2={TC}A_2 = \{TC\}
  • A3={CT}A_3 = \{CT\}
  • A4={CC}A_4 = \{CC\}

la x intende che non sappiamo se è testa o croce, quindi:

  • Tx=A1A2Tx = A_1 \cup A_2: il primo lancio è testa, secondo testa o croce
  • xT=A3A1xT = A_3 \cup A_1: il secondo lancio è testa, primo testa o croce

sostituendo al problema iniziale:

  1. P(TTTx)=P(A1A1A2)P(TT|Tx) = P(A_1|A_1 \cup A_2)
  2. P(TTTxxT)=P(A1A1A2A3)=P(A1A4c)P(TT|Tx \cup xT) = P(A_1|A_1 \cup A_2 \cup A_3) = P(A_1 | A_4^c)

Nel primo caso applichiamo bayes:

P(A1A1A2)=P(A1A2A1)P(A1)P(A1A2)=11/41/2=1/2P(A_1|A_1 \cup A_2) = P(A_1 \cup A_2 | A_1) \cdot \frac{P(A_1)}{P(A_1 \cup A_2)} = 1 \cdot \frac{1/4}{1/2} = 1/2

Nel secondo caso:

P(A1A4c)=P(A1A4c)P(A4c)=P(A1)P(A1)+P(A2)+P(A3)=1/43/4=1/3P(A_1 | A_4^c) = \frac{P(A_1 \cap A_4^c)}{P(A_4^c)} = \frac{P(A_1)}{P(A_1) + P(A_2) + P(A_3)} = \frac{1/4}{3/4} = 1/3

Osservazioni

  • Dati AA e BB eventi indipendenti, allora P(AB)=P(A)P(B)P(A\cap B) = P(A) \cdot P(B)
  • Dati AA e BB eventi indipendenti, allora P(AB)=P(A)P(A|B) = P(A). Dato che sono indipendenti, sapere che BB è accaduto non ci fornisce alcuna informazione su AA, quindi la probabilità di AA rimane invariata.
  • Se AA e BB sono due possibili esiti nello stesso spazio di probabilità, ed entrambi hanno probabilità non nulla, allora se P(AB)=0P(A\cap B) = 0, allora AA e BB NON sono indipendenti. Infatti se P(AB)=0P(A\cap B) = 0, e fossero indipendenti, allora P(AB)=P(A)P(B)P(A\cap B) = P(A) \cdot P(B), ma questo è impossibile se entrambi hanno probabilità non nulla. Quindi, se due eventi non possono verificarsi contemporaneamente (cioè sono mutuamente esclusivi), allora non possono essere indipendenti. Se uno dei due ha probabilità nulla, allora sono indipendenti.

Esercizio

Data una moneta, calcolare:

  1. Se si eseguono 5 lanci, quale è lo spazio degli eventi? Quale è la probabilità di ottenere 5 volte testa?
  2. Quale è la probabilità che al quinto lancio esca testa, sapendo che nei primi 4 lanci è uscita solo testa?
  3. Quale è la probabilità che al quinto lancio esca testa, sapendo che nei primi 4 lanci è uscita solo croce?

1: Lo spazio di un lancio è Ω={T,C}\Omega = \{T, C\}, quindi per 5 lanci lo spazio è Ω×Ω×Ω×Ω×Ω\Omega \times \Omega \times \Omega \times \Omega \times \Omega, ovvero l'insieme di tutte le possibili combinazioni di 5 lanci, che sono 25=322^5 = 32. La probabilità di ottenere 5 volte testa è P(T1×T2×T3×T4×T5)=P(T1T2T3T4T5)P(T_1 \times T_2 \times T_3 \times T_4 \times T_5) = P(T_1 \cap T_2 \cap T_3 \cap T_4 \cap T_5). Dato che i lanci sono indipendenti, questa probabilità è P(T1)P(T2)P(T3)P(T4)P(T5)=(1/2)5=1/32P(T_1) \cdot P(T_2) \cdot P(T_3) \cdot P(T_4) \cdot P(T_5) = (1/2)^5 = 1/32.

2: La probabilità dell'evento è data da

P(T5T1T2T3T4)=P(T5T1T2T3T4)P(XT1T2T3T4)==P(T5T1T2T3T4)P(T5T1T2T3T4)+P(CT1T2T3T4)==(1/2)5(1/2)5+(1/2)5=1/321/16=1/2\begin{aligned} P(T_5|T_1 \cap T_2 \cap T_3 \cap T_4) &= \frac{P(T_5 \cap T_1 \cap T_2 \cap T_3 \cap T_4)}{P(X \cap T_1 \cap T_2 \cap T_3 \cap T_4)} = \\ &= \frac{P(T_5 \cap T_1 \cap T_2 \cap T_3 \cap T_4)}{P(T_5 \cap T_1 \cap T_2 \cap T_3 \cap T_4) + P(C \cap T_1 \cap T_2 \cap T_3 \cap T_4)} = \\ &= \frac{(1/2)^5}{(1/2)^5 + (1/2)^5} = \frac{1/32}{1/16} = 1/2 \end{aligned}

3: Poiché gli eventi (i lanci) sono indipendenti, sapere che nei primi 4 lanci è uscita croce non altera la probabilità del quinto lancio. Analogo al caso precedente.

P(T5C1C2C3C4)=P(T5C1C2C3C4)P(XC1C2C3C4)==P(T5C1C2C3C4)P(T5C1C2C3C4)+P(CC1C2C3C4)==(1/2)5(1/2)5+(1/2)5=1/321/16=1/2\begin{aligned} P(T_5|C_1 \cap C_2 \cap C_3 \cap C_4) &= \frac{P(T_5 \cap C_1 \cap C_2 \cap C_3 \cap C_4)}{P(X \cap C_1 \cap C_2 \cap C_3 \cap C_4)} = \\ &= \frac{P(T_5 \cap C_1 \cap C_2 \cap C_3 \cap C_4)}{P(T_5 \cap C_1 \cap C_2 \cap C_3 \cap C_4) + P(C \cap C_1 \cap C_2 \cap C_3 \cap C_4)} = \\ &= \frac{(1/2)^5}{(1/2)^5 + (1/2)^5} = \frac{1/32}{1/16} = 1/2 \end{aligned}

Osservazione 1

3 Eventi AA, BB e CC sono indipendenti se e solo se:

P(ABC)=P(A)P(B)P(C)P(A\cap B \cap C) = P(A) \cdot P(B) \cdot P(C)

Osservazione 2

Date delle coppie, se ogni coppia di evento AA e BB è indipendente, non è detto che anche la tripla AA, BB e CC (o qualsiasi n-upla) sia indipendente.

Per esempio:

  • Ω={1,2,3,4}\Omega = \{1, 2, 3, 4\} equiprobabili
  • P(1)=P(2)=P(3)=P(4)=1/4P(1) = P(2) = P(3) = P(4) = 1/4
  • Gli eventi:
    • A={12}A = \{1 \cup 2\}
    • B={13}B = \{1 \cup 3\}
    • C={14}C = \{1 \cup 4\}

Sono indipendenti a coppie? Per dirlo dobbiamo verificare che per ogni coppia Ai,AjA_i, A_j (con iji \neq j) valga P(AiAj)=P(Ai)P(Aj)P(A_i \cap A_j) = P(A_i) \cdot P(A_j).

Calcoliamo le probabilità dei singoli eventi:

  • P(A)=2/4=1/2P(A) = 2/4 = 1/2
  • P(B)=2/4=1/2P(B) = 2/4 = 1/2
  • P(C)=2/4=1/2P(C) = 2/4 = 1/2

Calcoliamo le probabilità delle intersezioni a due a due:

  • AB={1}    P(AB)=1/4A \cap B = \{1\} \implies P(A \cap B) = 1/4 e P(A)P(B)=(1/2)(1/2)=1/4P(A)P(B) = (1/2)(1/2) = 1/4. Quindi AA e BB sono indipendenti.
  • AC={1}    P(AC)=1/4A \cap C = \{1\} \implies P(A \cap C) = 1/4 e P(A)P(C)=(1/2)(1/2)=1/4P(A)P(C) = (1/2)(1/2) = 1/4. Quindi AA e CC sono indipendenti.
  • BC={1}    P(BC)=1/4B \cap C = \{1\} \implies P(B \cap C) = 1/4 e P(B)P(C)=(1/2)(1/2)=1/4P(B)P(C) = (1/2)(1/2) = 1/4. Quindi BB e CC sono indipendenti.

Quindi sì, sono indipendenti a coppie.

Sono indipendenti a triple? Per essere indipendenti a triple (o globalmente), deve valere anche la condizione:

P(ABC)=P(A)P(B)P(C)P(A \cap B \cap C) = P(A) \cdot P(B) \cdot P(C)

Calcoliamo l'intersezione a tre: ABC={1}    P(ABC)=1/4A \cap B \cap C = \{1\} \implies P(A \cap B \cap C) = 1/4

Calcoliamo il prodotto delle tre probabilità: P(A)P(B)P(C)=(1/2)(1/2)(1/2)=1/8P(A) \cdot P(B) \cdot P(C) = (1/2) \cdot (1/2) \cdot (1/2) = 1/8

Poiché 1/41/81/4 \neq 1/8, concludiamo che NON sono indipendenti a triple. Questo esempio dimostra che l'indipendenza a coppie non implica l'indipendenza globale.

Esercizio

Uno studente sta svolgendo un esame da consegnare entro un ora.

La probabilità che lo studente consegni e termini l'esame entro xx ore è di x/2x/2 per 0x10 \leq x \leq 1.

Sapendo che lo studente sta ancora lavorando dopo 3/43/4 di ora, qual è la probabilità che userà l'inter ora?

  • ExE_x è l'evento "lo studente consegna entro xx ore"
  • FF è l'evento "lo studente impiega tutta l'ora"
  • P(FE3/4c)P(F|E_{3/4}^c) è la domanda del problema.

Allora:

P(FE3/4c)=P(FE3/4c)P(E3/4c)=P(F)P(E3/4c)\begin{aligned} P(F|E_{3/4}^c) &= \frac{P(F \cap E_{3/4}^c)}{P(E_{3/4}^c)} \\ &= \frac{P(F)}{P(E_{3/4}^c)} \end{aligned}

Abbiamo sostituito P(FE3/4c)P(F \cap E_{3/4}^c) con P(F)P(F) perché se lo studente impiega tutta l'ora, allora sicuramente non consegna entro 3/43/4 di ora. Allora possiamo applicare il complemento:

P(F)P(E3/4c)=P(F)1P(E3/4)\begin{aligned} \frac{P(F)}{P(E_{3/4}^c)} &= \frac{P(F)}{1 - P(E_{3/4})} \end{aligned}

Allora calcoliamo le probabilità dei nostri eventi:

  • P(F)=x2P(F) = \frac{x}{2} che per x=1x = 1 è P(F)=1/2P(F) = 1/2
  • P(E3/4)=342=38P(E_{3/4}) = \frac{\frac{3}{4}}{2} = \frac{3}{8}

Mettendo tutto insieme otteniamo:

P(FE3/4c)=1/213/8=1/25/8=45\begin{aligned} P(F|E_{3/4}^c) &= \frac{1/2}{1 - 3/8} = \frac{1/2}{5/8} = \frac{4}{5} \end{aligned}

Esercizio

Susanna è indecisa se frequentare il corso di francese o quello di chimica. Se susanna basa la sua decisione sull'esito di una moneta non truccata, e se ha la probabilità di 12\frac{1}{2} di prendere più di 27 se frequenta francese e probabilità di 23\frac{2}{3} di prendere più di 27 se frequenta chimica, con quale probabilità prenderà più di 27 all'esame di chimica? Con quale probabilità prenderà più di 27?

  • EfE_f è l'evento "Susanna frequenta francese"
  • EcE_c è l'evento "Susanna frequenta chimica"
  • AA è l'evento "Susanna prende più di 27"

Il problema chiede 2 cose:

  • P(EcA)P(E_c \cap A): la probabilità che Susanna prenda più di 27 all'esame di chimica, ovvero la probabilità che frequenti chimica e prenda più di 27
  • P(A)P(A): la probabilità che Susanna prenda più di 27, indipendentemente dal corso che frequenta

Allora:

  • P(Ef)=P(Ec)=1/2P(E_f) = P(E_c) = 1/2 perché Susanna decide con una moneta non truccata
  • P(AEc)=2/3P(A|E_c) = 2/3 perché se frequenta chimica ha questa probabilità di prendere più di 27
  • P(AEf)=1/2P(A|E_f) = 1/2 perché se frequenta francese ha questa probabilità di prendere più di 27 P(EcA)=P(AEc)P(Ec)=(2/3)(1/2)=1/3P(E_c \cap A) = P(A|E_c) \cdot P(E_c) = (2/3) \cdot (1/2) = 1/3
  • EcEf=ΩE_c \cup E_f = \Omega perché Susanna deve scegliere uno dei due corsi, quindi possiamo applicare la probabilità totale: P(A)=P(AEc)P(Ec)+P(AEf)P(Ef)=(2/3)(1/2)+(1/2)(1/2)=7/12P(A) = P(A|E_c) \cdot P(E_c) + P(A|E_f) \cdot P(E_f) = (2/3) \cdot (1/2) + (1/2) \cdot (1/2) = 7/12.

Regola del prodotto

Se AA e BB sono eventi, allora:

P(AB)=P(AB)P(B)P(A \cap B) = P(A|B) \cdot P(B)

Ma possiamo generalizzarlo a n eventi:

P(A1A2...An)=P(A1)P(A2A1)P(A3A1A2)...P(AnA1A2...An1)P(A_1 \cap A_2 \cap ... \cap A_n) = P(A_1) \cdot P(A_2|A_1) \cdot P(A_3|A_1 \cap A_2) \cdot ... \cdot P(A_n|A_1 \cap A_2 \cap ... \cap A_{n-1})

Esercizio

In un esame a risposta multipla uno studente può conoscere la risposta o tentare a caso. Se pp è la probabilità che lo studente conosca la risposta, e mm il numero di possibili risposte, sapendo che lo studente ha risposto correttamente a una domanda, qual è la probabilità che conoscesse la risposta?

  • pp probabilità che lo studente conosca la risposta
  • 1p1-p probabilità che lo studente non conosca la risposta e quindi tenti a caso
  • mm numero di possibili risposte
  • cc è l'evento "lo studente risponde correttamente"
  • kk è l'evento "lo studente conosce la risposta"

Allora dobbiamo calcolare P(kc)P(k|c), ovvero la probabilità che lo studente conoscesse la risposta dato che ha risposto correttamente.

Applichiamo Bayes:

P(kc)=P(ck)P(k)P(c)=1pP(c)\begin{aligned} P(k|c) &= \frac{P(c|k) \cdot P(k)}{P(c)} \\ &= \frac{1 \cdot p}{P(c)}\\ \end{aligned}

Dobbiamo calcolare P(c)P(c), ovvero la probabilità che lo studente risponda correttamente. Questa può accadere in due modi:

  1. Lo studente conosce la risposta e risponde correttamente (con probabilità pp)
  2. Lo studente non conosce la risposta e tenta a caso

Allora usiamo la probabilità totale per calcolare P(c)P(c):

P(c)=P(ck)P(k)+P(ckc)P(kc)=1p+1m(1p)=p+1pm=pm+1pm=1+p(m1)m\begin{aligned} P(c) &= P(c|k) \cdot P(k) + P(c|k^c) \cdot P(k^c) \\ &= 1 \cdot p + \frac{1}{m} \cdot (1-p) \\ &= p + \frac{1-p}{m} \\ &= \frac{pm + 1 - p}{m} \\ &= \frac{1 + p(m-1)}{m} \end{aligned}