Esercizio 6malattia XXX, Bayes e approssimazione di Poisson

Testo (I parziale, 17 novembre 2025, esercizio 2, Tema A; punti 9 = 3+3+3). Si ritiene che una certa malattia XXX colpisca lo $0{,}4\%$ della popolazione maschile e lo $0{,}5\%$ della popolazione femminile. Si sceglie un campione casuale di 2000 persone, di cui metà maschi e metà femmine.

(a) Qual è la probabilità che un individuo scelto a caso nel campione sia malato?

(b) Sapendo che l'individuo scelto è malato, qual è la probabilità che sia un maschio?

(c) Qual è la probabilità che il numero di femmine malate nel campione sia maggiore di 5? Approssimare tale valore usando l'approssimazione di Poisson. Lasciare indicate le somme ottenute.

Impostazione: dare un nome agli eventi

L'esperimento dei punti (a) e (b) è: "si sceglie a caso un individuo del campione". Gli eventi utili sono

• $M$ = "l'individuo è maschio", $F$ = "l'individuo è femmina";
• $A$ = "l'individuo è malato".

$M$ e $F$ sono complementari ($F = M^c$): formano una partizioneDue eventi disgiunti che insieme coprono tutto Ω: ogni individuo è maschio o femmina, mai entrambi.Spazio campionario ed eventi → di $\Omega$. È la struttura tipica "in due fasi" (prima il sesso, poi la malattia) delle probabilità totaliP(A) = P(A ∣ M) P(M) + P(A ∣ F) P(F).Formula delle probabilità totali e formula di Bayes →.

Traduzione dei dati.

• "Metà maschi e metà femmine", individuo scelto a caso: $\mathbb{P}(M) = \mathbb{P}(F) = \frac{1000}{2000} = 0{,}5$.
• "Lo $0{,}4\%$ della popolazione maschile" è una percentuale calcolata tra i maschi: è una probabilità condizionata, $\mathbb{P}(A \mid M) = 0{,}004$. Allo stesso modo $\mathbb{P}(A \mid F) = 0{,}005$.

Attenzione: $0{,}004$ non è $\mathbb{P}(A \cap M)$ ("maschio e malato"), che è più piccola.

(a) Probabilità che l'individuo sia malato

Non sappiamo il sesso dell'individuo, ma sappiamo la probabilità di malattia per ciascun sesso: è esattamente la situazione della formula delle probabilità totaliSe gli A_i formano una partizione, P(B) è la media delle P(B ∣ A_i) pesata con le P(A_i).Formula delle probabilità totali e formula di Bayes → con la partizione $\{M, F\}$:

$$\mathbb{P}(A) = \mathbb{P}(A \mid M)\, \mathbb{P}(M) + \mathbb{P}(A \mid F)\, \mathbb{P}(F) = 0{,}004 \cdot 0{,}5 + 0{,}005 \cdot 0{,}5 = 0{,}002 + 0{,}0025 = 0{,}0045$$

$$\boxed{\mathbb{P}(A) = 0{,}0045 = 0{,}45\%}$$

Controllo: è una media pesata di $0{,}004$ e $0{,}005$ con pesi uguali, quindi deve stare esattamente a metà ✓. I due addendi hanno anche un significato: $0{,}002 = \mathbb{P}(A \cap M)$ e $0{,}0025 = \mathbb{P}(A \cap F)$ (regola del prodotto).

(b) Probabilità che un malato sia maschio

Ora si conosce l'"effetto" (malato) e si chiede la "causa" (sesso): si inverte il condizionamento, da $\mathbb{P}(A \mid M)$ a $\mathbb{P}(M \mid A)$. È il compito della formula di BayesP(M ∣ A) = P(A ∣ M) P(M) / P(A).Formula delle probabilità totali e formula di Bayes →:

$$\mathbb{P}(M \mid A) = \frac{\mathbb{P}(A \mid M)\, \mathbb{P}(M)}{\mathbb{P}(A)} = \frac{0{,}004 \cdot 0{,}5}{0{,}0045} = \frac{0{,}002}{0{,}0045} = \frac{4}{9}$$

$$\boxed{\mathbb{P}(M \mid A) = \frac49 \approx 0{,}444}$$

Perché meno di $\frac12$. A priori maschi e femmine sono ugualmente probabili, ma le femmine si ammalano di più: tra i malati sono quindi in maggioranza ($\mathbb{P}(F \mid A) = \frac{0{,}0025}{0{,}0045} = \frac59$). Le due probabilità a posteriori sommano a $1$ ✓.

(c) Più di 5 femmine malate

La variabile aleatoria e la sua legge

Ora l'esperimento è diverso: si guarda l'intero campione. Sia $X$ = numero di femmine malate tra le 1000 femmine del campione.

Ogni femmina è malata con probabilità $p = 0{,}005$, e (implicitamente) le femmine si ammalano indipendentemente l'una dall'altra. Contare i "successi" (malate) in $n = 1000$ prove indipendenti con la stessa probabilità $p$ è esattamente lo schema del modello binomialeSe X conta i successi in n prove indipendenti con probabilità p, allora P(X = k) = (n su k) p^k (1−p)^(n−k).Bernoulli, binomiale e uniforme discreta →:

$$X \sim \text{Bin}(1000;\ 0{,}005), \qquad \mathbb{P}(X = k) = \binom{1000}{k} (0{,}005)^k (0{,}995)^{1000-k}, \quad k = 0, \dots, 1000$$

Passare al complementare

"Maggiore di 5" significa $X \in \{6, 7, \dots, 1000\}$: sono 995 termini. Il complementare $\{X \le 5\}$ ne ha solo 6, quindi

$$\mathbb{P}(X > 5) = 1 - \mathbb{P}(X \le 5) = 1 - \sum_{k=0}^{5} \binom{1000}{k} (0{,}005)^k (0{,}995)^{1000-k}$$

L'approssimazione di Poisson

I coefficienti $\binom{1000}{k}$ e le potenze $0{,}995^{1000-k}$ sono scomodi a mano. Quando $n$ è grande, $p$ è piccolo e $\lambda = np$ è moderato, la binomiale si approssima con la distribuzione di PoissonPoi(λ): P(Y = k) = e^(−λ) λ^k / k!, per k = 0, 1, 2, ...; approssima Bin(n, p) con λ = np se n è grande e p piccolo.Distribuzione di Poisson → di parametro $\lambda = np$. Qui

$$n = 1000 \ \text{(grande)}, \qquad p = 0{,}005 \ \text{(piccolo)}, \qquad \lambda = np = 1000 \cdot 0{,}005 = 5$$

Le condizioni sono soddisfatte. Con $Y \sim \text{Poi}(5)$:

$$\mathbb{P}(X > 5) \simeq \mathbb{P}(Y > 5) = 1 - \mathbb{P}(Y \le 5) = \boxed{1 - e^{-5} \sum_{k=0}^{5} \frac{5^k}{k!}}$$

Il valore numerico (calcolato con Python)

$$e^{-5} \sum_{k=0}^5 \frac{5^k}{k!} = e^{-5}\left(1 + 5 + \frac{25}{2} + \frac{125}{6} + \frac{625}{24} + \frac{3125}{120}\right) \approx 0{,}6160 \implies \mathbb{P}(Y > 5) \approx 0{,}3840$$

Il valore binomiale esatto è $\mathbb{P}(X > 5) \approx 0{,}38404$: l'approssimazione è eccellente (differenza dell'ordine di $10^{-7}$ in questo caso particolare; in generale l'errore è dell'ordine di $p$).

Plausibilità: la media di $X$ è $np = 5$; "più di 5" è "più della media", e per una distribuzione quasi simmetrica attorno a 5 ci si aspetta poco meno di metà della probabilità (il valore 5 stesso si prende circa il $17{,}5\%$). $0{,}384$ è coerente.

Errori comuni

• Leggere $0{,}004$ come $\mathbb{P}(A \cap M)$ invece che come $\mathbb{P}(A \mid M)$.
• In (b) rispondere $0{,}004$ o $0{,}5$: Bayes combina entrambi.
• In (c) usare $n = 2000$ o $p = 0{,}0045$: la domanda riguarda solo le femmine, quindi $n = 1000$ e $p = 0{,}005$.
• Confondere "maggiore di 5" con "almeno 5": $\mathbb{P}(X > 5) = 1 - \mathbb{P}(X \le 5)$, la somma arriva a $k = 5$ incluso.
• Usare la Poisson con $\lambda = p$ invece di $\lambda = np$.