Salta al contenuto
Note per Studenti Esercizio 56 · approssimazione normale di binomiali e Poisson

Esercizio 56approssimazione normale di binomiali e Poisson

Esame
In questa pagina 9

Testo.

  1. (II compitino, 4 febbraio 2015, seguito degli esercizi 49 e 50.) Per il cromosoma mutato (B(20 000,10−4)B(20\,000, 10^{-4})), è ragionevole approssimare P(X≤2)P(X \le 2) con una variabile normale? Di quali parametri?
  2. (Appello 4, 2014/15.) Il partito VivaVenezia ha avuto il 18%18\% dei voti; un exit poll intervista 10001000 elettori che votano indipendentemente, e XX è il numero di quelli che hanno votato VivaVenezia. Con un'approssimazione normale, stimare la probabilità che la società dichiari per il partito una percentuale sbagliata di almeno un punto percentuale.
  3. Approssimando X∼B(1000,0.002)X \sim B(1000, 0.002) con una normale, fornire un valore approssimato di P(X≤3)P(X \le 3).
  4. Un docente dichiara erroneamente sufficiente una prova incerta con probabilità 15\frac15. Su 3737 prove incerte, probabilità che vengano promosse ingiustamente almeno 55 persone? Approssimare con una Poisson e con una normale.
  5. Una macchina fotografica ordinata online è difettosa con probabilità 3%3\%. (a) Probabilità che su 10001000 spedite ve ne siano al più 2020 difettose? Quale variabile si considera? (b) Approssimarla con un'altra variabile discreta; (c) con una continua.
  6. Le chiamate a un centralino in un'ora sono una Poisson di media 250250. (a) Probabilità di più di 260260 chiamate? (b) Approssimarla con la normale.
  7. Il 35%35\% dell'elettorato è a favore di Pinco Pallino; in una sezione votano 200200 persone scelte a caso e XX conta quelle a suo favore. (a) Probabilità che X≥75X \ge 75; (b) approssimazione con il TCL.

Teoria: Teorema centrale del limite e correzione di continuitàSe X₁, X₂, … sono i.i.d. con media μ e varianza σ², la somma standardizzata (X₁ + ⋯ + Xₙ − nμ)/√(nσ²) tende in distribuzione a N(0,1): P(Sₙ ≤ a) ≈ Φ((a − nμ)/√(nσ²)) per n grande, qualunque sia la legge delle Xᵢ (servono indipendenza e stessa legge). Per la media campionaria: media μ, varianza σ²/n. Casi particolari: B(n,p) ≈ N(np, np(1−p)) per n grande; Po(λ) ≈ N(λ, λ) per λ ≥ 50. Per variabili a valori interi si usa la correzione di continuità: P(X ≤ k) ≈ P(Y ≤ k + 0.5), P(X = k) ≈ P(k − 0.5 ≤ Y ≤ k + 0.5).Teorema centrale del limite e correzione di continuità →, Variabile di Poisson e approssimazione della binomialeX ~ Po(λ), λ > 0, assume i valori 0, 1, 2, … con P(X = k) = e^{−λ} λᵏ/k!; λ è il numero medio di eventi. Nasce come limite della binomiale: se Xₙ ~ B(n, λ/n) allora P(Xₙ = k) → e^{−λ}λᵏ/k!. In pratica B(n,p) ≈ Po(np) se n è grande e p piccolo (regole del corso: n ≥ 20, p ≤ 0.05, p ≤ 10/n). Modella conteggi di eventi rari in un intervallo di tempo o spazio: telefonate, arrivi, terremoti, difetti, vincite.Variabile di Poisson e approssimazione della binomiale →.


Promemoria: B(n,p)≈N(np,np(1−p))B(n, p) \approx N(np, np(1 - p)) se npnp e n(1−p)n(1 - p) sono grandi; Po⁡(λ)≈N(λ,λ)\operatorname{Po}(\lambda) \approx N(\lambda, \lambda) per λ≥50\lambda \ge 50; correzione di continuità per variabili intere. Valori esatti e approssimati calcolati con Python.

1. Il cromosoma: la normale non va bene

np=2np = 2 e np(1−p)≈2np(1 - p) \approx 2: la media è troppo piccola. La binomiale è fortemente asimmetrica (vale 00, 11, 22 con probabilità 0.1350.135, 0.2710.271, 0.2710.271, e non può essere negativa), mentre una N(2,2)N(2, 2) darebbe probabilità non trascurabile ai valori negativi. Non è ragionevole: qui va usata la Poisson (nn grande, pp piccolo, npnp moderato).

2. L'exit poll

X∼B(1000,0.18)X \sim B(1000, 0.18): media 180180, varianza 1000⋅0.18⋅0.82=147.61000\cdot 0.18\cdot 0.82 = 147.6, deviazione standard 12.1512.15. La percentuale stimata è X1000\frac{X}{1000}; sbagliare di almeno un punto significa ∣X1000−0.18∣≥0.01\left|\frac{X}{1000} - 0.18\right| \ge 0.01, cioè ∣X−180∣≥10|X - 180| \ge 10, cioè X≤170X \le 170 oppure X≥190X \ge 190.

  • Senza correzione: 2(1−Φ(1012.15))=2(1−Φ(0.82))≈0.412\left(1 - \Phi\left(\frac{10}{12.15}\right)\right) = 2(1 - \Phi(0.82)) \approx 0.41.
  • Con correzione (X≤170→Y≤170.5X \le 170 \to Y \le 170.5, X≥190→Y≥189.5X \ge 190 \to Y \ge 189.5): 2(1−Φ(9.512.15))=2(1−Φ(0.78))≈0.4342\left(1 - \Phi\left(\frac{9.5}{12.15}\right)\right) = 2(1 - \Phi(0.78)) \approx 0.434, praticamente uguale al valore esatto (0.4340.434).

Con 10001000 intervistati l'errore di un punto percentuale capita quasi una volta su due.

3. B(1000,0.002)B(1000, 0.002) e la normale

Media 22, varianza 1000⋅0.002⋅0.998=1.9961000\cdot 0.002\cdot 0.998 = 1.996: P(X≤3)≈Φ(3.5−21.996)=Φ(1.06)≈0.856P(X \le 3) \approx \Phi\left(\frac{3.5 - 2}{\sqrt{1.996}}\right) = \Phi(1.06) \approx 0.856 Il valore esatto è 0.8570.857: qui il risultato è buono, ma con np=2np = 2 la normale è in generale poco affidabile (vedi il punto 1); la Poisson Po⁡(2)\operatorname{Po}(2) dà P(Y≤3)=e−2(1+2+2+43)≈0.857P(Y \le 3) = e^{-2}\left(1 + 2 + 2 + \frac43\right) \approx 0.857.

4. Le promozioni ingiuste

X∼B(37,15)X \sim B\left(37, \frac15\right): esatto P(X≥5)≈0.888P(X \ge 5) \approx 0.888.

  • Poisson con λ=37⋅0.2=7.4\lambda = 37\cdot 0.2 = 7.4: 1−∑k=04e−7.47.4kk!≈0.8601 - \sum_{k=0}^{4}e^{-7.4}\frac{7.4^k}{k!} \approx 0.860 (non eccellente: p=0.2p = 0.2 non è piccolo).
  • Normale con media 7.47.4 e varianza 37⋅0.2⋅0.8=5.9237\cdot 0.2\cdot 0.8 = 5.92: P(X≥5)≈1−Φ(4.5−7.45.92)=Φ(1.19)≈0.883P(X \ge 5) \approx 1 - \Phi\left(\frac{4.5 - 7.4}{\sqrt{5.92}}\right) = \Phi(1.19) \approx 0.883 (migliore).

5. Le macchine fotografiche

(a) X∼B(1000,0.03)X \sim B(1000, 0.03): P(X≤20)=∑k=020(1000k)0.03k 0.971000−k≈0.0333P(X \le 20) = \sum_{k=0}^{20}\binom{1000}{k}0.03^k\,0.97^{1000-k} \approx 0.0333.

(b) Poisson Po⁡(30)\operatorname{Po}(30) (nn grande, p=0.03p = 0.03 piccolo): ∑k=020e−3030kk!≈0.0353\sum_{k=0}^{20}e^{-30}\frac{30^k}{k!} \approx 0.0353.

(c) Normale N(30,29.1)N(30, 29.1): Φ(20.5−3029.1)=Φ(−1.76)≈0.039\Phi\left(\frac{20.5 - 30}{\sqrt{29.1}}\right) = \Phi(-1.76) \approx 0.039. Nelle code lontane (qui 1.81.8 deviazioni standard sotto la media) l'approssimazione normale è meno precisa.

6. Il centralino

(a) X∼Po⁡(250)X \sim \operatorname{Po}(250): P(X>260)=1−∑k=0260e−250250kk!≈0.251P(X > 260) = 1 - \sum_{k=0}^{260}e^{-250}\frac{250^k}{k!} \approx 0.251 (calcolo al computer).

(b) λ=250≥50\lambda = 250 \ge 50, quindi X≈N(250,250)X \approx N(250, 250); con la correzione, X>260  ⟺  X≥261→Y≥260.5X > 260 \iff X \ge 261 \to Y \ge 260.5: P(X>260)≈1−Φ(10.5250)=1−Φ(0.66)≈0.253P(X > 260) \approx 1 - \Phi\left(\frac{10.5}{\sqrt{250}}\right) = 1 - \Phi(0.66) \approx 0.253

7. Pinco Pallino

(a) X∼B(200,0.35)X \sim B(200, 0.35): P(X≥75)=∑k=75200(200k)0.35k 0.65200−k≈0.251P(X \ge 75) = \sum_{k=75}^{200}\binom{200}{k}0.35^k\,0.65^{200-k} \approx 0.251.

(b) Media 7070, varianza 200⋅0.35⋅0.65=45.5200\cdot 0.35\cdot 0.65 = 45.5, deviazione standard 6.756.75: P(X≥75)≈1−Φ(74.5−706.75)=1−Φ(0.67)≈0.252P(X \ge 75) \approx 1 - \Phi\left(\frac{74.5 - 70}{6.75}\right) = 1 - \Phi(0.67) \approx 0.252


Riepilogo: quale approssimazione

Situazione Approssimazione
B(n,p)B(n, p), nn grande, pp piccolo, npnp moderato Po⁡(np)\operatorname{Po}(np)
B(n,p)B(n, p), npnp e n(1−p)n(1-p) grandi N(np,np(1−p))N(np, np(1-p)) con correzione
Po⁡(λ)\operatorname{Po}(\lambda), λ≥50\lambda \ge 50 N(λ,λ)N(\lambda, \lambda) con correzione

Errori comuni

  • Usare la normale con npnp piccolo.
  • Sbagliare il verso della correzione (per X≥190X \ge 190 si usa 189.5189.5, per X≤170X \le 170 si usa 170.5170.5).

Versione ripasso

Testo.

  1. (II compitino, 4 febbraio 2015, seguito degli esercizi 49 e 50.) Per il cromosoma mutato (B(20 000,10−4)B(20\,000, 10^{-4})), è ragionevole approssimare P(X≤2)P(X \le 2) con una variabile normale? Di quali parametri?
  2. (Appello 4, 2014/15.) Il partito VivaVenezia ha avuto il 18%18\% dei voti; un exit poll intervista 10001000 elettori che votano indipendentemente, e XX è il numero di quelli che hanno votato VivaVenezia. Con un'approssimazione normale, stimare la probabilità che la società dichiari per il partito una percentuale sbagliata di almeno un punto percentuale.
  3. Approssimando X∼B(1000,0.002)X \sim B(1000, 0.002) con una normale, fornire un valore approssimato di P(X≤3)P(X \le 3).
  4. Un docente dichiara erroneamente sufficiente una prova incerta con probabilità 15\frac15. Su 3737 prove incerte, probabilità che vengano promosse ingiustamente almeno 55 persone? Approssimare con una Poisson e con una normale.
  5. Una macchina fotografica ordinata online è difettosa con probabilità 3%3\%. (a) Probabilità che su 10001000 spedite ve ne siano al più 2020 difettose? Quale variabile si considera? (b) Approssimarla con un'altra variabile discreta; (c) con una continua.
  6. Le chiamate a un centralino in un'ora sono una Poisson di media 250250. (a) Probabilità di più di 260260 chiamate? (b) Approssimarla con la normale.
  7. Il 35%35\% dell'elettorato è a favore di Pinco Pallino; in una sezione votano 200200 persone scelte a caso e XX conta quelle a suo favore. (a) Probabilità che X≥75X \ge 75; (b) approssimazione con il TCL.

Teoria: Teorema centrale del limite e correzione di continuitàSe X₁, X₂, … sono i.i.d. con media μ e varianza σ², la somma standardizzata (X₁ + ⋯ + Xₙ − nμ)/√(nσ²) tende in distribuzione a N(0,1): P(Sₙ ≤ a) ≈ Φ((a − nμ)/√(nσ²)) per n grande, qualunque sia la legge delle Xᵢ (servono indipendenza e stessa legge). Per la media campionaria: media μ, varianza σ²/n. Casi particolari: B(n,p) ≈ N(np, np(1−p)) per n grande; Po(λ) ≈ N(λ, λ) per λ ≥ 50. Per variabili a valori interi si usa la correzione di continuità: P(X ≤ k) ≈ P(Y ≤ k + 0.5), P(X = k) ≈ P(k − 0.5 ≤ Y ≤ k + 0.5).Teorema centrale del limite e correzione di continuità →, Variabile di Poisson e approssimazione della binomialeX ~ Po(λ), λ > 0, assume i valori 0, 1, 2, … con P(X = k) = e^{−λ} λᵏ/k!; λ è il numero medio di eventi. Nasce come limite della binomiale: se Xₙ ~ B(n, λ/n) allora P(Xₙ = k) → e^{−λ}λᵏ/k!. In pratica B(n,p) ≈ Po(np) se n è grande e p piccolo (regole del corso: n ≥ 20, p ≤ 0.05, p ≤ 10/n). Modella conteggi di eventi rari in un intervallo di tempo o spazio: telefonate, arrivi, terremoti, difetti, vincite.Variabile di Poisson e approssimazione della binomiale →.

Promemoria: B(n,p)≈N(np,np(1−p))B(n,p)\approx N(np,np(1-p)) se npnp e n(1−p)n(1-p) sono grandi; Po⁡(λ)≈N(λ,λ)\operatorname{Po}(\lambda)\approx N(\lambda,\lambda) per λ≥50\lambda\ge50; correzione di continuità.

  1. np=2np = 2: la binomiale (0.1350.135, 0.2710.271, 0.2710.271 in 0,1,20,1,2) è asimmetrica e non negativa, una N(2,2)N(2,2) darebbe massa ai negativi: non ragionevole, serve la Poisson.
  2. X∼B(1000,0.18)X\sim B(1000,0.18): media 180180, varianza 147.6147.6, σ=12.15\sigma = 12.15. Errore di almeno un punto: ∣X−180∣≥10|X - 180|\ge10, cioè X≤170X\le170 o X≥190X\ge190. Senza correzione 2(1−Φ(0.82))≈0.412(1 - \Phi(0.82))\approx0.41; con correzione 2(1−Φ(9.512.15))=2(1−Φ(0.78))≈0.4342\left(1 - \Phi\left(\frac{9.5}{12.15}\right)\right) = 2(1 - \Phi(0.78))\approx0.434 (come il valore esatto).
  3. Media 22, varianza 1.9961.996: Φ(3.5−21.996)=Φ(1.06)≈0.856\Phi\left(\frac{3.5 - 2}{\sqrt{1.996}}\right) = \Phi(1.06)\approx0.856 (esatto 0.8570.857; la Poisson Po⁡(2)\operatorname{Po}(2) dà 0.8570.857, ma con npnp piccolo la normale è poco affidabile).
  4. X∼B(37,15)X\sim B(37,\frac15), esatto P(X≥5)≈0.888P(X\ge5)\approx0.888. Poisson λ=7.4\lambda = 7.4: ≈0.860\approx0.860 (pp non piccolo). Normale (media 7.47.4, varianza 5.925.92): Φ(7.4−4.55.92)=Φ(1.19)≈0.883\Phi\left(\frac{7.4 - 4.5}{\sqrt{5.92}}\right) = \Phi(1.19)\approx0.883, migliore.
  5. (a) X∼B(1000,0.03)X\sim B(1000,0.03): P(X≤20)≈0.0333P(X\le20)\approx0.0333. (b) Po⁡(30)\operatorname{Po}(30): ≈0.0353\approx0.0353. (c) N(30,29.1)N(30,29.1): Φ(20.5−3029.1)=Φ(−1.76)≈0.039\Phi\left(\frac{20.5 - 30}{\sqrt{29.1}}\right) = \Phi(-1.76)\approx0.039 (nelle code lontane è meno precisa).
  6. (a) P(Po⁡(250)>260)≈0.251P(\operatorname{Po}(250) > 260)\approx0.251. (b) N(250,250)N(250,250), X≥261→Y≥260.5X\ge261\to Y\ge260.5: 1−Φ(10.5250)=1−Φ(0.66)≈0.2531 - \Phi\left(\frac{10.5}{\sqrt{250}}\right) = 1 - \Phi(0.66)\approx0.253.
  7. (a) X∼B(200,0.35)X\sim B(200,0.35): P(X≥75)≈0.251P(X\ge75)\approx0.251. (b) Media 7070, varianza 45.545.5: 1−Φ(74.5−706.75)=1−Φ(0.67)≈0.2521 - \Phi\left(\frac{74.5 - 70}{6.75}\right) = 1 - \Phi(0.67)\approx0.252.
Situazione Approssimazione
B(n,p)B(n,p), nn grande, pp piccolo, npnp moderato Po⁡(np)\operatorname{Po}(np)
B(n,p)B(n,p), npnp e n(1−p)n(1-p) grandi N(np,np(1−p))N(np,np(1-p)) con correzione
Po⁡(λ)\operatorname{Po}(\lambda), λ≥50\lambda\ge50 N(λ,λ)N(\lambda,\lambda) con correzione

Errori tipici: usare la normale con npnp piccolo; sbagliare il verso della correzione (X≥190→189.5X\ge190\to189.5, X≤170→170.5X\le170\to170.5).

Esercizi su questo argomento

Teoria collegata