Estimación Puntual

El uso de una estadística para estimar, o hallar un valor cercano al valor real, de un parámetro. Ejemplo: Uso del promedio muestal $\overline{X}$, para acercarnos al promedio poblacional $\mu$.

Intervalo de Confianza

Intervalo en el cual esperamos que esté el valor del parámetro con una confiabilidad $\alpha$ pre-establecida.

Distribución de Muestreo

Distribución de probabilidad de una estadística basada en muestras aleatorias.

Estimación por intervalo

Para una variable aleatoria $X$. Queremos estimar un parámetro. ¿cómo sabemos en que rango de valores estará el parámetro? Queremos tener:

• Alta confiabilidad $\alpha$ : que el parámetro muy seguramente esté en el intervalo.
• Intervalos angostos
• No conocemos los parámetros ($\mu$, $\sigma$)
• Podemos estimar en muestras ($\overline{x}$, $S$)
• Sabemos el tamaño de la muestra, $n$.

Intervalo de Confianza para $\mu$

(La fórmula se desarrolla en Contento 2019, pg. 274 )

Dado un valor de confianza $(1 - \alpha) \cdot 100 \%$, un tamaño de muestra $n$, una desviación estándar de la muestra $s$, y un promedio de la muestra de $\overline{x}$, el intervalo de confianza es:

$\overline{x} \pm t_{\alpha/2;n-1}\frac{s}{\sqrt{n}}$

• $t_{\alpha/2;n-1}$ es una distribución t-student con $n-1$ grados de libertad.

• $\alpha / 2$ es la probabilidad a cola derecha.

• Para calcularla en python:

  • Se usa la librería scipy.stats y de allí se calcula con la función t.ppf.

  • Los argumentos de la función son $1-\alpha$ y los grados de libertad

  • Por ejemplo, se calcula $t_{\alpha/2;n-1}$, para $\alpha=0.1$ y $n=10$, así:

  from scipy.stats import t
  t.ppf(1-0.05,10)
  

  y si lo queremos usar después se puede asignar el valor a una variable.

Ejemplo, Contento 2019, pg 275.

Una agencia ambiental mide el promedio de masa ($\mu g$) en 1 $m^3$ de aire. Si hay $n=5$ mediciones, con $\overline{x}=61$, y $s=5.244$, ¿cuál es el intervalo de confianza de 95% ?

Respuesta

Como queremos que la confianza sea del 95%, entonces tenemos que: $(1 - \alpha) 100 = 95\%$

Esto quiere decir que $\alpha = 0.05$ y por lo tanto $\alpha/2 = 0.025$

Como son 5 mediciones, tenemos que $n-1 = 5-1 = 4$

Con estos datos calculamos la distribución $t$:

$t_{0.025;4}$ se calcula así:

from scipy.stats import t
t.ppf(1-0.025,4)

$\overline{x}-t_{0.025;4}(5.244/2.36)=54.49$

$\overline{x}+t_{0.025;4}(5.244/2.36)=67.51$

El intervalo sería, entonces: (54.49,67.51)

Recursos

• Estimación Puntual y por Intervalo https://youtu.be/fcf4MZq0WiE?si=WSx-usuTEq8rZcvq
• Ejemplo https://youtu.be/a7xO1dtJcTE?si=L2QTHugAVSAVSefs
• Contento, capítulo 5