Método

Ejercicios: Test de bondad de ajuste

χ2=(OE)2E,gl=k1p\chi^{2} = \sum \frac{(O-E)^{2}}{E}, \quad gl = k - 1 - p

Si no te acordás cómo se usa, repasá la teoría.

Las 12 actividades recorren los distintos tipos de ejercicio del tema, con tres variantes de cada uno. Si te trabás en una, mirá la solución y seguí con la siguiente.

1. Se tiró un dado 180180 veces y las caras del 11 al 66 salieron 2626; 2525; 3838; 3535; 3232; 2424 veces. De tabla, el valor crítico al 5%5\% con 55 grados de libertad es 11,070511{,}0705. Decidí si el dado es equilibrado.

Ver solución

1. Planteamos las hipótesis. H0H_0: las seis caras son igualmente probables. H1H_1: al menos una no lo es.

2. Calculamos las frecuencias esperadas. Bajo la nula cada cara tiene probabilidad 16\dfrac{1}{6}, así que E=1806=30E = \dfrac{180}{6} = 30 para las seis.

3. Verificamos la condición. Todas las esperadas valen 3030, muy por encima de 55: el test se puede aplicar.

4. Calculamos el estadístico. χ2=(OE)2E=(2630)230+(2530)230+(3830)230+(3530)230+(3230)230+(2430)230=5,6667\chi^{2} = \sum \dfrac{(O - E)^{2}}{E} = \dfrac{(26 - 30)^{2}}{30} + \dfrac{(25 - 30)^{2}}{30} + \dfrac{(38 - 30)^{2}}{30} + \dfrac{(35 - 30)^{2}}{30} + \dfrac{(32 - 30)^{2}}{30} + \dfrac{(24 - 30)^{2}}{30} = 5{,}6667

5. Damos los grados de libertad. k1p=610=5k - 1 - p = 6 - 1 - 0 = 5. No se descuenta nada más porque no se estimó ningún parámetro de los datos: la distribución propuesta estaba completamente especificada.

6. Decidimos. 5,66675{,}6667 no supera a 11,070511{,}0705, así que no se rechaza: los datos son compatibles con un dado equilibrado.

Respuesta: χ2=5,6667\chi^{2} = 5{,}6667 contra 11,070511{,}0705 con 55 grados: no se rechaza el ajuste

2. Se probó el ajuste de una distribución a una muestra agrupada en 66 clases, estimando 11 parámetro a partir de los propios datos. El estadístico dio 10,67810{,}678. De tabla, los valores críticos al 5%5\% son 9,48779{,}4877 con 44 grados y 11,070511{,}0705 con 55. Decidí el test con los grados correctos.

Ver solución

1. Escribimos la regla de los grados de libertad. gl=k1pgl = k - 1 - p, con kk la cantidad de clases y pp la cantidad de parámetros estimados de los datos.

2. Calculamos. gl=611=4gl = 6 - 1 - 1 = 4.

3. Explicamos por qué se descuenta. Cada parámetro estimado de la muestra hace que el modelo ya venga ajustado a esos mismos datos. El ajuste va a parecer mejor de lo que es, y descontar grados de libertad compensa ese efecto.

4. Decidimos con los grados correctos. 10,67810{,}678 supera a 9,48779{,}4877: se rechaza el ajuste.

5. Vemos qué pasaría con los grados sin descontar. Con 55 grados el crítico sería 11,070511{,}0705, y la decisión sería no rechazar: la conclusión cambia.

6. En qué dirección se equivoca el error. Usar más grados de los que corresponde da un crítico más alto y hace más difícil rechazar: el modelo pasa el test con más facilidad de la que merece.

Respuesta: Con gl=4gl = 4 y crítico 9,48779{,}4877: se rechaza el ajuste

3. Un modelo teórico predice que los casos se reparten en las categorías Bajo, Medio, Alto con probabilidades 0,250{,}25, 0,50{,}5 y 0,250{,}25. En 8080 casos observados se contaron 1717; 4040; 2323. De tabla, el crítico al 5%5\% con 22 grados es 5,99155{,}9915. Realizá el test de bondad de ajuste.

Ver solución

1. Planteamos las hipótesis. H0H_0: los datos siguen el reparto propuesto. H1H_1: no lo siguen.

2. Calculamos las frecuencias esperadas. Bajo: 800,25=2080 \cdot 0{,}25 = 20 · Medio: 800,5=4080 \cdot 0{,}5 = 40 · Alto: 800,25=2080 \cdot 0{,}25 = 20, que suman 8080.

3. Verificamos la condición. La esperada más chica es 2020, mayor que 55: el test se puede aplicar.

4. Calculamos el estadístico. (1720)220+(4040)240+(2320)220=0,9\dfrac{(17 - 20)^{2}}{20} + \dfrac{(40 - 40)^{2}}{40} + \dfrac{(23 - 20)^{2}}{20} = 0{,}9

5. Damos los grados de libertad. k1p=310=2k - 1 - p = 3 - 1 - 0 = 2: no se estimó ningún parámetro, las probabilidades venían del modelo.

6. Decidimos. 0,90{,}9 no supera a 5,99155{,}9915, así que no se rechaza: los datos son compatibles con el reparto propuesto.

Respuesta: χ2=0,9\chi^{2} = 0{,}9 contra 5,99155{,}9915: no se rechaza el ajuste

4. Se quiere probar si la cantidad de fallas por turno sigue una Poisson de parámetro 1,21{,}2, con 100100 turnos observados agrupados en las clases 00, 11, 22, 33, 44 o más. Calculá las frecuencias esperadas y verificá si se cumple la condición del test.

Ver solución

1. Calculamos las probabilidades de cada clase bajo el modelo. P(X=0)=0,3012P(X = 0) = 0{,}3012 · P(X=1)=0,3614P(X = 1) = 0{,}3614 · P(X=2)=0,2169P(X = 2) = 0{,}2169 · P(X=3)=0,0867P(X = 3) = 0{,}0867 · P(X4)=10,9662=0,0338P(X \geq 4) = 1 - 0{,}9662 = 0{,}0338.

2. Verificamos que sumen uno. 11, como corresponde: las clases cubren todos los valores posibles.

3. Calculamos las esperadas. 00: 1000,301230,1194100 \cdot 0{,}3012 \approx 30{,}1194 · 11: 1000,361436,1433100 \cdot 0{,}3614 \approx 36{,}1433 · 22: 1000,216921,686100 \cdot 0{,}2169 \approx 21{,}686 · 33: 1000,08678,6744100 \cdot 0{,}0867 \approx 8{,}6744 · 44 o más: 1000,03383,3769100 \cdot 0{,}0338 \approx 3{,}3769

4. Revisamos la condición. 44 o más: 3,37693{,}3769, que no llegan a 55.

5. Decidimos cómo seguir. Hay que agrupar las clases con esperada chica con la vecina antes de calcular el estadístico. Al hacerlo baja la cantidad de clases, y con ella los grados de libertad.

6. Anotamos los grados de libertad que corresponderían. k1pk - 1 - p, con p=0p = 0 si el parámetro 1,21{,}2 viene de afuera, o p=1p = 1 si se lo estimó con el promedio de los propios datos.

Respuesta: No se cumple en la clase 44 o más, con 3,37693{,}3769: hay que agrupar clases

5. Se tiró un dado 300300 veces y las caras del 11 al 66 salieron 4242; 4848; 4343; 5050; 4646; 7171 veces. De tabla, el valor crítico al 5%5\% con 55 grados de libertad es 11,070511{,}0705. Decidí si el dado es equilibrado.

Ver solución

1. Planteamos las hipótesis. H0H_0: las seis caras son igualmente probables. H1H_1: al menos una no lo es.

2. Calculamos las frecuencias esperadas. Bajo la nula cada cara tiene probabilidad 16\dfrac{1}{6}, así que E=3006=50E = \dfrac{300}{6} = 50 para las seis.

3. Verificamos la condición. Todas las esperadas valen 5050, muy por encima de 55: el test se puede aplicar.

4. Calculamos el estadístico. χ2=(OE)2E=(4250)250+(4850)250+(4350)250+(5050)250+(4650)250+(7150)250=11,48\chi^{2} = \sum \dfrac{(O - E)^{2}}{E} = \dfrac{(42 - 50)^{2}}{50} + \dfrac{(48 - 50)^{2}}{50} + \dfrac{(43 - 50)^{2}}{50} + \dfrac{(50 - 50)^{2}}{50} + \dfrac{(46 - 50)^{2}}{50} + \dfrac{(71 - 50)^{2}}{50} = 11{,}48

5. Damos los grados de libertad. k1p=610=5k - 1 - p = 6 - 1 - 0 = 5. No se descuenta nada más porque no se estimó ningún parámetro de los datos: la distribución propuesta estaba completamente especificada.

6. Decidimos. 11,4811{,}48 supera a 11,070511{,}0705, así que se rechaza: hay evidencia de que el dado no es equilibrado.

Respuesta: χ2=11,48\chi^{2} = 11{,}48 contra 11,070511{,}0705 con 55 grados: se rechaza el ajuste

6. Se probó el ajuste de una distribución a una muestra agrupada en 55 clases, estimando 11 parámetro a partir de los propios datos. El estadístico dio 13,48313{,}483. De tabla, los valores críticos al 5%5\% son 7,81477{,}8147 con 33 grados y 9,48779{,}4877 con 44. Decidí el test con los grados correctos.

Ver solución

1. Escribimos la regla de los grados de libertad. gl=k1pgl = k - 1 - p, con kk la cantidad de clases y pp la cantidad de parámetros estimados de los datos.

2. Calculamos. gl=511=3gl = 5 - 1 - 1 = 3.

3. Explicamos por qué se descuenta. Cada parámetro estimado de la muestra hace que el modelo ya venga ajustado a esos mismos datos. El ajuste va a parecer mejor de lo que es, y descontar grados de libertad compensa ese efecto.

4. Decidimos con los grados correctos. 13,48313{,}483 supera a 7,81477{,}8147: se rechaza el ajuste.

5. Vemos qué pasaría con los grados sin descontar. Con 44 grados el crítico sería 9,48779{,}4877, y la decisión sería rechazar: la misma, pero por casualidad.

6. En qué dirección se equivoca el error. Usar más grados de los que corresponde da un crítico más alto y hace más difícil rechazar: el modelo pasa el test con más facilidad de la que merece.

Respuesta: Con gl=3gl = 3 y crítico 7,81477{,}8147: se rechaza el ajuste

7. Un modelo teórico predice que los casos se reparten en las categorías Bajo, Medio, Alto con probabilidades 0,250{,}25, 0,50{,}5 y 0,250{,}25. En 160160 casos observados se contaron 3434; 7474; 5252. De tabla, el crítico al 5%5\% con 22 grados es 5,99155{,}9915. Realizá el test de bondad de ajuste.

Ver solución

1. Planteamos las hipótesis. H0H_0: los datos siguen el reparto propuesto. H1H_1: no lo siguen.

2. Calculamos las frecuencias esperadas. Bajo: 1600,25=40160 \cdot 0{,}25 = 40 · Medio: 1600,5=80160 \cdot 0{,}5 = 80 · Alto: 1600,25=40160 \cdot 0{,}25 = 40, que suman 160160.

3. Verificamos la condición. La esperada más chica es 4040, mayor que 55: el test se puede aplicar.

4. Calculamos el estadístico. (3440)240+(7480)280+(5240)240=4,95\dfrac{(34 - 40)^{2}}{40} + \dfrac{(74 - 80)^{2}}{80} + \dfrac{(52 - 40)^{2}}{40} = 4{,}95

5. Damos los grados de libertad. k1p=310=2k - 1 - p = 3 - 1 - 0 = 2: no se estimó ningún parámetro, las probabilidades venían del modelo.

6. Decidimos. 4,954{,}95 no supera a 5,99155{,}9915, así que no se rechaza: los datos son compatibles con el reparto propuesto.

Respuesta: χ2=4,95\chi^{2} = 4{,}95 contra 5,99155{,}9915: no se rechaza el ajuste

8. Se quiere probar si la cantidad de fallas por turno sigue una Poisson de parámetro 1,21{,}2, con 150150 turnos observados agrupados en las clases 00, 11, 22, 33, 44 o más. Calculá las frecuencias esperadas y verificá si se cumple la condición del test.

Ver solución

1. Calculamos las probabilidades de cada clase bajo el modelo. P(X=0)=0,3012P(X = 0) = 0{,}3012 · P(X=1)=0,3614P(X = 1) = 0{,}3614 · P(X=2)=0,2169P(X = 2) = 0{,}2169 · P(X=3)=0,0867P(X = 3) = 0{,}0867 · P(X4)=10,9662=0,0338P(X \geq 4) = 1 - 0{,}9662 = 0{,}0338.

2. Verificamos que sumen uno. 11, como corresponde: las clases cubren todos los valores posibles.

3. Calculamos las esperadas. 00: 1500,301245,1791150 \cdot 0{,}3012 \approx 45{,}1791 · 11: 1500,361454,215150 \cdot 0{,}3614 \approx 54{,}215 · 22: 1500,216932,529150 \cdot 0{,}2169 \approx 32{,}529 · 33: 1500,086713,0116150 \cdot 0{,}0867 \approx 13{,}0116 · 44 o más: 1500,03385,0653150 \cdot 0{,}0338 \approx 5{,}0653

4. Revisamos la condición. Todas las esperadas superan 55, así que la condición se cumple.

5. Decidimos cómo seguir. Se puede calcular el estadístico directamente sobre estas cinco clases.

6. Anotamos los grados de libertad que corresponderían. k1pk - 1 - p, con p=0p = 0 si el parámetro 1,21{,}2 viene de afuera, o p=1p = 1 si se lo estimó con el promedio de los propios datos.

Respuesta: Se cumple: todas las esperadas superan 55

9. Se tiró un dado 180180 veces y las caras del 11 al 66 salieron 3636; 2727; 2929; 2525; 2929; 3434 veces. De tabla, el valor crítico al 5%5\% con 55 grados de libertad es 11,070511{,}0705. Decidí si el dado es equilibrado.

Ver solución

1. Planteamos las hipótesis. H0H_0: las seis caras son igualmente probables. H1H_1: al menos una no lo es.

2. Calculamos las frecuencias esperadas. Bajo la nula cada cara tiene probabilidad 16\dfrac{1}{6}, así que E=1806=30E = \dfrac{180}{6} = 30 para las seis.

3. Verificamos la condición. Todas las esperadas valen 3030, muy por encima de 55: el test se puede aplicar.

4. Calculamos el estadístico. χ2=(OE)2E=(3630)230+(2730)230+(2930)230+(2530)230+(2930)230+(3430)230=2,9333\chi^{2} = \sum \dfrac{(O - E)^{2}}{E} = \dfrac{(36 - 30)^{2}}{30} + \dfrac{(27 - 30)^{2}}{30} + \dfrac{(29 - 30)^{2}}{30} + \dfrac{(25 - 30)^{2}}{30} + \dfrac{(29 - 30)^{2}}{30} + \dfrac{(34 - 30)^{2}}{30} = 2{,}9333

5. Damos los grados de libertad. k1p=610=5k - 1 - p = 6 - 1 - 0 = 5. No se descuenta nada más porque no se estimó ningún parámetro de los datos: la distribución propuesta estaba completamente especificada.

6. Decidimos. 2,93332{,}9333 no supera a 11,070511{,}0705, así que no se rechaza: los datos son compatibles con un dado equilibrado.

Respuesta: χ2=2,9333\chi^{2} = 2{,}9333 contra 11,070511{,}0705 con 55 grados: no se rechaza el ajuste

10. Se probó el ajuste de una distribución a una muestra agrupada en 66 clases, estimando 22 parámetros a partir de los propios datos. El estadístico dio 5,0865{,}086. De tabla, los valores críticos al 5%5\% son 7,81477{,}8147 con 33 grados y 11,070511{,}0705 con 55. Decidí el test con los grados correctos.

Ver solución

1. Escribimos la regla de los grados de libertad. gl=k1pgl = k - 1 - p, con kk la cantidad de clases y pp la cantidad de parámetros estimados de los datos.

2. Calculamos. gl=612=3gl = 6 - 1 - 2 = 3.

3. Explicamos por qué se descuenta. Cada parámetro estimado de la muestra hace que el modelo ya venga ajustado a esos mismos datos. El ajuste va a parecer mejor de lo que es, y descontar grados de libertad compensa ese efecto.

4. Decidimos con los grados correctos. 5,0865{,}086 no supera a 7,81477{,}8147: no se rechaza el ajuste.

5. Vemos qué pasaría con los grados sin descontar. Con 55 grados el crítico sería 11,070511{,}0705, y la decisión sería no rechazar: la misma, pero por casualidad.

6. En qué dirección se equivoca el error. Usar más grados de los que corresponde da un crítico más alto y hace más difícil rechazar: el modelo pasa el test con más facilidad de la que merece.

Respuesta: Con gl=3gl = 3 y crítico 7,81477{,}8147: no se rechaza el ajuste

11. Un modelo teórico predice que los casos se reparten en las categorías Bajo, Medio, Alto con probabilidades 0,250{,}25, 0,50{,}5 y 0,250{,}25. En 8080 casos observados se contaron 2020; 3535; 2525. De tabla, el crítico al 5%5\% con 22 grados es 5,99155{,}9915. Realizá el test de bondad de ajuste.

Ver solución

1. Planteamos las hipótesis. H0H_0: los datos siguen el reparto propuesto. H1H_1: no lo siguen.

2. Calculamos las frecuencias esperadas. Bajo: 800,25=2080 \cdot 0{,}25 = 20 · Medio: 800,5=4080 \cdot 0{,}5 = 40 · Alto: 800,25=2080 \cdot 0{,}25 = 20, que suman 8080.

3. Verificamos la condición. La esperada más chica es 2020, mayor que 55: el test se puede aplicar.

4. Calculamos el estadístico. (2020)220+(3540)240+(2520)220=1,875\dfrac{(20 - 20)^{2}}{20} + \dfrac{(35 - 40)^{2}}{40} + \dfrac{(25 - 20)^{2}}{20} = 1{,}875

5. Damos los grados de libertad. k1p=310=2k - 1 - p = 3 - 1 - 0 = 2: no se estimó ningún parámetro, las probabilidades venían del modelo.

6. Decidimos. 1,8751{,}875 no supera a 5,99155{,}9915, así que no se rechaza: los datos son compatibles con el reparto propuesto.

Respuesta: χ2=1,875\chi^{2} = 1{,}875 contra 5,99155{,}9915: no se rechaza el ajuste

12. Se quiere probar si la cantidad de fallas por turno sigue una Poisson de parámetro 1,21{,}2, con 200200 turnos observados agrupados en las clases 00, 11, 22, 33, 44 o más. Calculá las frecuencias esperadas y verificá si se cumple la condición del test.

Ver solución

1. Calculamos las probabilidades de cada clase bajo el modelo. P(X=0)=0,3012P(X = 0) = 0{,}3012 · P(X=1)=0,3614P(X = 1) = 0{,}3614 · P(X=2)=0,2169P(X = 2) = 0{,}2169 · P(X=3)=0,0867P(X = 3) = 0{,}0867 · P(X4)=10,9662=0,0338P(X \geq 4) = 1 - 0{,}9662 = 0{,}0338.

2. Verificamos que sumen uno. 11, como corresponde: las clases cubren todos los valores posibles.

3. Calculamos las esperadas. 00: 2000,301260,2388200 \cdot 0{,}3012 \approx 60{,}2388 · 11: 2000,361472,2866200 \cdot 0{,}3614 \approx 72{,}2866 · 22: 2000,216943,372200 \cdot 0{,}2169 \approx 43{,}372 · 33: 2000,086717,3488200 \cdot 0{,}0867 \approx 17{,}3488 · 44 o más: 2000,03386,7538200 \cdot 0{,}0338 \approx 6{,}7538

4. Revisamos la condición. Todas las esperadas superan 55, así que la condición se cumple.

5. Decidimos cómo seguir. Se puede calcular el estadístico directamente sobre estas cinco clases.

6. Anotamos los grados de libertad que corresponderían. k1pk - 1 - p, con p=0p = 0 si el parámetro 1,21{,}2 viene de afuera, o p=1p = 1 si se lo estimó con el promedio de los propios datos.

Respuesta: Se cumple: todas las esperadas superan 55

Las cuentas de este material están verificadas una por una. Si aun así encontrás algo raro, .