Distribucions de probabilitat

Més enllà del dau: models probabilístics

Formació d'Aprofundiment en Probabilitat per a Professorat

De l'Histograma a la Corba

Distribució de les alçades d'una població:

1. Mida de la mostra (N): Augmenta els individus per a tenir dades més representatives.

2. Amplada de l'interval (w): Fes més estretes les classes per afinar el detall i evitar la pèrdua d'informació.

Densitat de probabilitat

$f(x) \approx$ Freqüència relativa Amplada interval (w)
A mesura que la mostra tendeix a infinit ($N \to \infty$) i l'interval tendeix a zero ($w \to 0$), el perfil aspre de l'histograma s'adapta perfectament a una funció de densitat contínua i simètrica. Aquesta corba (en vermell) és la Campana de Gauss.

El Catàleg de Distribucions

Formes, equacions i aplicacions

La distribució normal o Gaussiana: $N(\mu, \sigma)$

La forma de la campana depèn de dos valors:

Desplaça el centre de la distribució.
Controla la concentració de les dades, l'amplada de la campana.
$$ f(x) = \frac{1}{\sigma\sqrt{2\pi}} e^{-\frac{1}{2}\left(\frac{x-\mu}{\sigma}\right)^2} $$

Reptes: El pes de les taronges

Paràmetres del lot

Probabilitat = 0.0000

Càlcul de Probabilitats

$X \sim N(\mu, \sigma)$

La probabilitat es defineix com l'àrea sota la corba de densitat.

Prem la fletxa DRETA (→) per descobrir les propietats.

I com calculem aquestes àrees?

1. El mur de la integració

La funció de densitat de la Normal no té una primitiva elemental. Això vol dir que no podem aplicar la famosa Regla de Barrow per calcular-ne l'àrea a mà:

$$ P(a \le X \le b) = \int_{a}^{b} \frac{1}{\sigma\sqrt{2\pi}} e^{-\frac{1}{2}\left(\frac{x-\mu}{\sigma}\right)^2} dx = \text{???} $$

2. La solució (i el seu defecte)

Com que no es pot resoldre analíticament, es va recórrer a mètodes numèrics per calcular les àrees i imprimir-les en taules de probabilitat.

Però atenció: tenim infinites combinacions de $\mu$ i $\sigma$...

Necessitaríem infinites taules per cobrir tots els casos!

3. La Clau Universal

La solució brillant va ser imprimir una única taula per a la distribució més senzilla de totes: la Normal Estàndard.

$$ Z \sim N(0,1) $$

Qualsevol altra distribució $X$ haurà de patir una translació i un escalat (tipificar-se) per poder entrar a la taula.

Lectura de la Taula $N(0,1)$

Cerca de la probabilitat acumulada $P(Z \le z)$


P( Z ≤ 1.25 ) = 0.8944

La Normal Tipificada: De X a Z

Com utilitzar la Taula Estàndard

Hi ha infinites distribucions normals, però només tenim una taula de probabilitats: la de la Normal Estàndard $N(0,1)$.

Per poder utilitzar-la amb qualsevol distribució $N(\mu, \sigma)$, hem de tipificar la variable.

Prem la fletxa DRETA (→) per descobrir els passos.

Original: $X \sim N(\mu, \sigma)$
Tipificada: $Z \sim N(0, 1)$

Exemple Pràctic: Despesa Turística

La despesa d'un turista a Barcelona segueix una distribució Normal amb $\mu = 180€$ i $\sigma = 30€$.

$X \sim N(180, 30)$

Prem la fletxa DRETA (→) per començar.

$z_1 = $ -4.00
$z_2 = $ 4.00
Probabilitat: 1.0000

El Teorema Central del Límit

La suma de variables independents

1. Un sol dau (Uniforme)

Si tirem un dau clàssic perfecte, tots els valors de l'1 al 6 tenen exactament la mateixa probabilitat de sortir ($1/6 \approx 16.6\%$). No hi ha cap mena de campana, només un bloc recte.

$$ X_1 \sim \text{Uniforme}(1, 6) $$
Pas 0 de 3

L'Espai $L^2$: La geometria de l'atzar

La independència com a perpendicularitat

1. Variables a l'espai $L^2$

En probabilitat, una variable aleatòria $X$ centrada és un vector geomètriquíssim. La seva desviació típica ($\sigma_X$) n'és exactament la longitud o norma, definida per l'arrel de la integral del seu quadrat:

$$ \|X\| = \sqrt{\int_{\Omega} X^2 dP} = \sqrt{E[X^2]} = \sigma_X = 4 $$
Pas 0 de 4

Per què la desviació típica creix tan a poc a poc?

El secret de la independència geomètrica

1. El parany lineal

La intuïció ens enganya: pensem que si sumem $n$ daus, el "desordre" o dispersió total s'hauria de multiplicar per $n$. Però això només passaria si els daus estiguessin completament sincronitzats!

$$ \sigma_{\text{total}} \neq n \cdot \sigma \implies \sigma_{\text{total}} = \sqrt{n} \cdot \sigma $$
Pas 0 de 3

La distribució binomial

Fins ara mesuràvem variables contínues (pes, alçada...).

Ara passem a l'estadística de comptar èxits.

1. L'Experiment de Bernoulli

És l'experiment aleatori més senzill que existeix: aquell on només hi ha dos resultats possibles i excloents.

Èxit: passa amb una probabilitat $p$
Fracàs: passa amb una probabilitat $q = 1 - p$
1. Exp. Bernoulli Només 2 resultats Prob.: $p$ i $q$

2. L'exemple i la mostra ($n$)

En unes eleccions, aturar a 1 persona per preguntar si ha votat la candidata és un exp. de Bernoulli (Ex: Èxit $\implies p=0.4$).

Si agafem una mostra de $n$ persones:
$x = \text{"nº de persones que l'han votat"}$ $x\sim B(n,p)$
1. Exp. Bernoulli Només 2 resultats Prob.: $p$ i $q$
2. La Mostra $n$ persones Variable $x$

3. La gran pregunta

Si $x$ és el nombre d'èxits...

Quina és la probabilitat d'aconseguir exactament $k$ èxits?
$P(x = k) = ?$

La Fórmula de la Binomial

Volem calcular la probabilitat d'obtenir exactament $k$ èxits en $n$ experiments.

$P(X = k) = ?$

Anem a construir la fórmula pas a pas.

1. Quantes combinacions possibles hi ha?

Els $k$ èxits poden aparèixer en qualsevol ordre al llarg dels $n$ experiments. Per comptar totes les maneres d'escollir $k$ posicions d'un total de $n$ (sense que n'importi l'ordre), fem servir el nombre combinatori:

$\binom{n}{k} = \frac{n!}{k!(n-k)!}$
1. Combinacions Ordre no importa Múltiple: $\binom{n}{k}$

2. La probabilitat de cada camí

A cada combinació concreta, l'èxit ($p$) es repeteix $k$ vegades, i el fracàs ($q$) es repeteix la resta de vegades ($n - k$). Com que són experiments independents, les probabilitats es multipliquen:

$p^k \cdot q^{n-k}$
1. Nombre de camins Combinatòria $\binom{n}{k}$
2. Probabilitat per camí Independents $p^k \cdot q^{n-k}$

3. La Fórmula Definitiva

Si multipliquem els camins possibles per la probabilitat de cadascun d'ells, obtenim la funció de probabilitat de la distribució Binomial.

$P(X = k) = \binom{n}{k} \cdot p^k \cdot q^{n-k}$

Exemple: Càlcul de probabilitats

Context del problema

El 10% dels telèfons d'un cert model s'espatllen durant el període de garantia. D'aquests, el 70% es pot reparar i, la resta, s'ha de canviar. Una empresa compra 20 telèfons mòbils d'aquest model. La probabilitat que un s'hagi de canviar és $p=0.1\cdot 0.3=0.03$ (i $q=0.97$).

$x \sim B(20 \, , \, 0.03)$

1. Cas puntual: $P(x=5)$

Quina probabilitat hi ha que s'hagin de canviar exactament 5 telèfons?

$P(x = 5) = \binom{20}{5} \cdot 0.03^5 \cdot 0.97^{15}$

$P(x = 5) \approx 0.00015$

Nota: Com que la $p$ és molt petita, és molt difícil que n'hi hagi 5 de defectuosos.

2. Escenaris múltiples: $P(x=1 \text{ o } 2)$

Quina probabilitat hi ha que n'haguem de canviar un o dos?

$P(x=1 \text{ o } 2) = P(x=1) + P(x=2)$
$\binom{20}{1} \cdot 0.03^1 \cdot 0.97^{19} + \binom{20}{2} \cdot 0.03^2 \cdot 0.97^{18}$

$\approx 0.3364 + 0.0988 = \mathbf{0.4352}$

3. Més dificultat: "Com a mínim dos" ($x \ge 2$)

Calcular $P(x=2) + P(x=3) + \dots + P(x=20)$ és massa llarg. Fem servir el contrari:

$P(x \ge 2) = 1 - P(x < 2) = 1 - [P(x=0) + P(x=1)]$
$1 - \left( \binom{20}{0} 0.97^{20} + \binom{20}{1} 0.03^1 0.97^{19} \right)$

$1 - (0.5438 + 0.3364) = \mathbf{0.1198}$

Com és la Distribució Binomial?

1. Naturalesa Discreta

La variable només pot prendre valors enters ($x = 0, 1, 2 \dots$).

Si preguntem a 20 persones, no en podem trobar "2,5" que votin un candidat. Per això el gràfic està fet de barres separades.

1. Discreta Només valors enters

2. La Mitjana o Esperança ($\mu$)

És el "centre" del gràfic. Indica quants èxits esperem aconseguir de mitjana.

$\mu = n \cdot p$
1. Discreta Només enters
2. Mitjana ($\mu$) $n \cdot p$

3. La Desviació Típica ($\sigma$)

Mesura l'amplada i la dispersió de les probabilitats respecte al valor central.

$\sigma = \sqrt{n \cdot p \cdot q}$


D'on surt la fórmula $\sigma = \sqrt{n \cdot p \cdot q}$ ?

1. La incertesa d'una sola persona ($p \cdot q$)

Matemàticament, la mesura d'incertesa (variància) d'un experiment de Bernoulli és la multiplicació de $p$ i $q$. Per què?

  • Si $p=0.99$, gairebé segur que sortirà Èxit. Hi ha poca incertesa ($0.99 \cdot 0.01 = 0.0099$).
  • Si $p=0.50$, és pura loteria. Hi ha la màxima incertesa possible ($0.5 \cdot 0.5 = 0.25$).
1. Variància individual Incertesa = $p \cdot q$

2. Variància de la suma

Com hem vist abans: si sumem experiments independents, les seves variàncies se sumen.

Si enquestem a $n$ persones, la variància total de la mostra serà sumar $p \cdot q$ un total de $n$ vegades:

Variància total $= pq + pq + pq + \dots = \mathbf{n \cdot p \cdot q}$
1. Variància 1 persona $p \cdot q$
2. Regla de la suma ($n$) Var. Total = $n \cdot p \cdot q$

3. Retorn a les unitats reals ($\sigma$)

Com que la variància està elevada al quadrat (ens donaria "vots al quadrat" o "persones al quadrat"), hem de fer l'arrel quadrada per trobar la desviació típica i parlar en unitats normals.

$\sigma = \sqrt{n \cdot p \cdot q}$

De la Binomial a la Normal

1. El límit de la calculadora

Quan $n$ és molt gran, els càlculs es fan inviables perquè els nombres es fan enormes ($\binom{1000}{300}$ no es pot fer).

1. Problema $n$ gran = ERROR

2. Aproximació per una normal

Comprova-ho al gràfic de la dreta: a mesura que puges la $n$, els rectangles es fusionen i dibuixen una campana de Gauss gairebé perfecta.

1. Problema $n$ gran = ERROR
2. Solució Forma de campana

3. La Regla d'Aproximació

Si es compleix que $n \ge 30$, podem abandonar l'estadística discreta i utilitzar directament la distribució Normal:

$B(n, p) \implies N(np \, , \sqrt{n p q})$
1. Problema$n$ gran = ERROR
2. SolucióForma de campana
3. Regla$B(n,p) \rightarrow N(\mu, \sigma)$

4. La Correcció de Yates

Per calcular $P(X \le k)$, fixa't al gràfic on acaba la barra de $k$. Com que té amplada 1, físicament s'estén fins a $k + 0.5$!

$P(X \le \mathbf{k}) \approx P(X' \le \mathbf{k + 0.5})$



🟦 Binomial 🟪 Correcció Yates 🟥 Normal

El Teorema de De Moivre-Laplace

Només importa la n?

Als llibres sovint hi diu que podem aproximar la Binomial a la Normal només si $n \ge 30$.

Però, què passa si la probabilitat d'èxit és molt petita o molt gran? Prova de baixar la $p$ a 0.05 al gràfic de la dreta i observa què li passa a la "campana".

El Mur de la Binomial

La Binomial no pot tenir valors inferiors a $0$ ni superiors a $n$.

Si $p$ és molt extrema, la mitjana $\mu = n \cdot p$ queda massa a prop del límit. La campana xoca contra el mur i queda deformada (asimètrica), fent que l'aproximació sigui invàlida per molt gran que sigui la $n$.

El Criteri de De Moivre-Laplace

Per assegurar-nos que la campana cau prou lluny dels murs i és ben simètrica, exigim que compleixi estrictament aquestes tres condicions:

$n \ge 30$
$n \cdot p \ge 5$
$n \cdot q \ge 5$


n ≥ 30 n·p = 25.0 n·q = 25.0