Test Hipotesis 1
Test Hipotesis 1
Test Hipotesis 1
Contrastes de
hipótesis
Profesor: Daniel Isidoro Ramírez
Email: [email protected]
Departamento de Métodos Estadísticos
Universidad de Zaragoza
Inferencia estadística
𝑯𝟎
𝜷
𝛼/2
𝛼/2
𝑃 =1 −𝛼
Test de hipótesis
Test de hipótesis:
Población Normal(1) Población NO-normal
Prueba
𝑫𝒊𝒔𝒕𝒓𝒊𝒃𝒖𝒄𝒊𝒐𝒏𝒆𝒔 𝒕, 𝝌𝟐 𝒚 𝑭 Comparaciones sobre las medianas
Sobre una media:
Test bilateral: 𝑥ഥ = 𝜇0 Test de Wilcoxon para una muestra
Test t para una muestra
Test unilateral: 𝑥ഥ ≤ 𝜇0 𝑜 𝑥ഥ ≥ 𝜇0 (también llamado Mann-Whitney)
t9; 0.025=-2.26
Normalmente se hace la comparación
tc=-3.16 con los valores positivos:
𝛼/2 ഥ − 𝝁𝟎
𝒙
¿ 𝑬𝒔 ≷ 𝒕𝒏−𝟏;𝟎.𝟗𝟕𝟓 ?
𝛼/2 0 𝑺Τ 𝒏
Poblaciones Normales
𝒕𝒄 = −𝟑. 𝟏𝟔
+2.26
𝛼
0
Poblaciones Normales
0.2 0.2
0.18 0.18
0.16 0.16
0.14 0.14
0.12 0.12
0.1 0.1
0.08 0.08
0.06 0.06
2 2
0.04 𝜒𝑛−1;𝛼 Τ2 2
𝜒𝑛−1;
0.04 𝜒𝑛−1;1−𝛼
0.02 1−𝛼Τ2 0.02
0 0
0 2 4 6 8 10 12 14 16 18 20 0 2 4 6 8 10 12 14 16 18 20
Poblaciones Normales
0.2
0.18
0.16 Frecuentemente las Tablas de la distribución
0.14 𝜒 2 dan el valor que deja una probabilidad 𝛼 a
0.12 su derecha, en lugar de a su izquierda:
2
𝑃 𝜒 2 > 𝜒𝑛−1; 𝛼 =𝛼
0.1
0.08
0.06
0.04 2 2
𝜒𝑛−1; 1−𝛼 ՜ 𝜒𝑛−1; 𝛼
0.02
0
0 2 4 6 8 10 12 14 16 18 20
Poblaciones Normales
𝑛−1 · 𝑆 2
2º Calculamos 𝜒𝑐2 =
𝜎02
2 𝑛−1 · 𝑆 2 2
𝑆𝑖 𝜒𝑛−1; 𝛼 Τ2 ≤ 𝜒𝑐2 = ≤ 𝜒𝑛−1;1− 𝛼 Τ2 ֜ 𝐴𝑐𝑒𝑝𝑡𝑎𝑚𝑜𝑠 𝐻𝑜
𝜎02
3º Decisión
𝑛−1 · 𝑆 2 2 𝑛−1 · 𝑆 2 2
𝑆𝑖 𝜒𝑐2 = > 𝜒𝑛−1;1− 𝛼 Τ2 𝑜 𝜒𝑐2 = < 𝜒𝑛−1; 𝛼 Τ2 ֜ 𝑅𝑒𝑐ℎ𝑎𝑧𝑎𝑚𝑜𝑠 𝐻𝑜
𝜎02 𝜎02
𝑛−1 · 𝑆 2 2
𝑆𝑖 𝜒𝑐2 = ≤ 𝜒𝑛−1;1−𝛼 ֜ 𝐴𝑐𝑒𝑝𝑡𝑎𝑚𝑜𝑠 𝐻𝑜
𝜎02
3º Decisión
𝑛−1 · 𝑆 2 2
𝑆𝑖 𝜒𝑐2 = > 𝜒𝑛−1;1−𝛼 ֜ 𝑅𝑒𝑐ℎ𝑎𝑧𝑎𝑚𝑜𝑠 𝐻𝑜
𝜎02
𝑛−1 · 𝑆 2 2
𝑆𝑖 𝜒𝑐2 = ≥ 𝜒𝑛−1; 𝛼 ֜ 𝐴𝑐𝑒𝑝𝑡𝑎𝑚𝑜𝑠 𝐻𝑜
𝜎02
3º Decisión
𝑛−1 · 𝑆 2 2
𝑆𝑖 𝜒𝑐2 = < 𝜒𝑛−1; 𝛼 ֜ 𝑅𝑒𝑐ℎ𝑎𝑧𝑎𝑚𝑜𝑠 𝐻𝑜
𝜎02
Comparamos el cociente de las varianzas muestrales (𝑆12 Τ𝑆22 𝑐𝑜𝑛 𝑆12 > 𝑆22 ) con 1
𝑆12
Por tanto ~ 𝐹𝑛1−1; 𝑛2−1 . Fijamos un N.S. α (α = 0.05 e.g.)
𝑆22
1º H0: 𝜎12 = 𝜎22 Test bilateral o de dos colas (con N.S. α):
HA: 𝜎12 ≠ 𝜎22 Región de aceptación: 𝐹𝑛1 −1; 𝑛2 −1; 𝛼Τ2 ; 𝐹𝑛1 −1; 𝑛2 −1; 1−𝛼Τ2
𝑆12
2º Calculamos 𝐹𝑐 =
𝑆22
𝑆𝑖 𝛼 Τ2 ≤ 𝑃 𝐹 < 𝐹𝑐 ≤ 1 − 𝛼 Τ2 ֜ 𝐴𝑐𝑒𝑝𝑡𝑎𝑚𝑜𝑠 𝐻𝑜
3º O bien ቐ
𝑆𝑖 𝑃 𝐹 < 𝐹𝑐 < 𝛼 Τ2 𝑜 𝑃 𝐹 < 𝐹𝑐 > 1 − 𝛼 Τ2 ֜ 𝑅𝑒𝑐ℎ𝑎𝑧𝑎𝑚𝑜𝑠 𝐻𝑜
Como en la distribución 𝜒 2 , las Tablas de la distribución F suelen dar el valor que deja una probabilidad
𝛼 a su derecha, en lugar de a su izquierda: P 𝐹 > 𝐹𝑛1 −1; 𝑛2 −1; 𝛼 = 𝛼
Poblaciones Normales
Comparamos el cociente de las varianzas muestrales (𝑆12 Τ𝑆22 𝑐𝑜𝑛 𝑆12 > 𝑆22 ) con 1
𝑆12
Por tanto ~ 𝐹𝑛1−1; 𝑛2−1 . Fijamos un N.S. α (α = 0.05 e.g.)
𝑆22
3º O bien ቐ 0,1
𝛼
𝑆𝑖 𝑃 𝐹 < 𝐹𝑐 > 1 − 𝛼 ֜ 𝑅𝑒𝑐ℎ𝑎𝑧𝑎𝑚𝑜𝑠 𝐻𝑜 0
0 1 2 3 4 5
Poblaciones Normales
Comparamos el cociente de las varianzas muestrales (𝑆12 Τ𝑆22 𝑐𝑜𝑛 𝑆12 > 𝑆22 ) con 1
𝑆12
Por tanto ~ 𝐹𝑛1−1; 𝑛2−1 . Fijamos un N.S. α (α = 0.05 e.g.)
𝑆22
3º O bien ቐ 0,1
𝛼
𝑆𝑖 𝑃 𝐹 < 𝐹𝑐 < 𝛼 ֜ 𝑅𝑒𝑐ℎ𝑎𝑧𝑎𝑚𝑜𝑠 𝐻𝑜 0
0 1 2 3 4 5
Poblaciones Normales
Se utilizan test distintos según que las varianzas de las 2 muestras puedan considerarse
iguales o no (𝑆12 Τ𝑆22 = 1 ՜ 𝑃𝑟𝑢𝑒𝑏𝑎 𝐹)
𝑥1 −𝑥2 𝑆12 𝑆22
Si 𝜎12 = 𝜎22 : ~ 𝑡𝑛1 +𝑛2 −2 𝐸𝑟𝑟𝑜𝑟 𝑒𝑠𝑡á𝑛𝑑𝑎𝑟 𝑑𝑒 𝑙𝑎 𝑑𝑖𝑓𝑒𝑟𝑒𝑛𝑐𝑖𝑎: 𝑠𝑥1 −𝑥2 = 𝑛1
+ 𝑛2
𝑆2 2
1 + 𝑆2
𝑛1 𝑛2
𝑆2 2
1 + 𝑆2
𝑥1 −𝑥2 𝑛1 𝑛2
Si 𝜎12 ≠ 𝜎22 : ~ 𝑡𝜈 𝑑𝑜𝑛𝑑𝑒 𝜈 = 2 2 2 [solo lo utilizaremos en Rcmdr]
𝑆2 2 𝑠2
1 ൗ𝑛1 + 𝑠2 ൗ𝑛2
1 + 𝑆2
𝑛1 𝑛2 𝑛1 −1 𝑛2 −1
Comparamos la media muestral (𝑑)ҧ con 0 donde la muestra procede de una población
normal de varianza σ2 desconocida que estimamos como S2(d)
ത
𝑑−0
Por tanto ~𝑡𝑛−1 .
𝑆𝑑 Τ 𝑛
𝐻0 : 𝐸 𝑑ҧ = 0 𝐻0 : 𝐸 𝑑ҧ ≤ 0 𝐻0 : 𝐸 𝑑ҧ ≥ 0
1º 𝑇𝑒𝑠𝑡 𝑑𝑒 𝑑𝑜𝑠 𝑐𝑜𝑙𝑎𝑠 ቐ ; 𝑇𝑒𝑠𝑡 𝑑𝑒 𝑢𝑛𝑎 𝑐𝑜𝑙𝑎 ቐ 𝑜ቐ
𝐻𝐴 : 𝐸 𝑑ҧ = 0 𝐻𝐴 : 𝐸 𝑑ҧ > 0 𝐻𝐴 : 𝐸 𝑑ҧ < 0
𝑑ത
2º Calculamos 𝑡𝑐 =
𝑆𝑑 Τ 𝑛
𝐶𝑜𝑚𝑝𝑎𝑟𝑎𝑚𝑜𝑠 𝑡𝑐 ≶ 𝑡𝑛−1
3º Decisión 𝑆𝑖 𝑡𝑐 ∈ 𝑅𝑒𝑔𝑖ó𝑛 𝑑𝑒 𝑎𝑐𝑒𝑝𝑡𝑎𝑐𝑖ó𝑛 ֜ 𝐴𝑐𝑒𝑝𝑡𝑎𝑚𝑜𝑠 𝐻𝑜
𝑆𝑖 𝑡𝑐 ∉ 𝑅𝑒𝑔𝑖ó𝑛 𝑑𝑒 𝑎𝑐𝑒𝑝𝑡𝑎𝑐𝑖ó𝑛 ֜ 𝑅𝑒𝑐ℎ𝑎𝑧𝑎𝑚𝑜𝑠 𝐻𝑜
Poblaciones Normales
Norte Sur
100 95
107 90 a) Poblaciones normales ; 𝜎12 = 𝜎22 :
115 102
𝑥ҧ 1 −𝑥ҧ2 𝑆12 𝑆22
96 100 ~𝑡𝑛1 +𝑛2 −2;1−𝛼Τ2 ; 𝐼𝐶1−𝛼 = 𝑥1ҧ − 𝑥ҧ2 ± 𝑡𝑛1 +𝑛2−2;1−𝛼Τ2 +
𝑛1 𝑛2
102 98 𝑆2 2
1 + 𝑆2
110 90 𝑛1 𝑛2
109 100
103 95 𝑥1ҧ − 𝑥ҧ2 = 8.4
99 88 1 − 𝛼 = 0.95֜𝑡𝑛1 +𝑛2 −2;1−𝛼Τ2 = 𝑡18;0.975 = 2.10
106 105 𝑰𝑪𝟎.𝟗𝟓 = 𝟑. 𝟎𝟐 , 𝟏𝟑. 𝟕𝟖
𝒙 104.7 96.3 𝑆12 𝑆22
+ = 2.56
𝑺 5.8 5.7 𝑛1 𝑛2
𝒏 10 10
S.E. 1.8 1.8
Poblaciones Normales
Norte Sur
100 95 b) H0: 𝑥ҧ1 = 𝑥ҧ 2 ¿ 𝑜 𝑝𝑜𝑑𝑟í𝑎 𝑠𝑒𝑟 𝑥ҧ1 > 𝑥ҧ 2 ?
107 90
115 102 Poblaciones normales ; 𝜎12 ≠ 𝜎22 :
96 100
102 98 𝑥ҧ1 − 𝑥ҧ 2
110 90 𝑡𝑐 = ~𝑡𝑛1 +𝑛2 −2 ;
109 100 𝑛1 − 1 ∙ 𝑆12 + 𝑛2 − 1 ∙ 𝑆22 1 1
103 95 ∙ +
99 88
𝑛1 + 𝑛2 − 2 𝑛1 𝑛2
106 105
𝒙 104.7 96.3 8.4 𝑃 𝑡18 > 𝑡𝑐 = 0.00015
𝑺 5.8 5.7
𝑡𝑐 = = 4.454 ֜ ቊ
1.886 𝑃 𝑡18 > 𝑡𝑐 = 0.00031
𝒏 10 10
S.E. 1.8 1.8
Media: La media de una muestra procedente de una población NO normal ¿Es igual a
un valor prefijado 𝜇0 ?. Muestra con n GRANDE (n >>30)
1 1
Muestra: 𝑥1 , 𝑥2 … 𝑥𝑛 ՜ 𝑥ҧ = σ𝑛𝑖=1 𝑥𝑖 ; 𝑆2 = σ𝑛 𝑥𝑖 − 𝑥ҧ 2
𝑛 𝑛−1 𝑖=1
Pruebas No Paramétricas
EJ.6. Una encuesta de tamaño 100 ha resultado con 56 individuos a favor de una
propuesta y 44 en contra. a) Con esos resultados ¿se puede aceptar al N.S. 𝛼 = 0,05
que la mitad de la población está de acuerdo con la propuesta (𝑝 = 0,5)?. b) Si otra
encuesta con 𝑛 = 56 da como resultado 21 individuos a favor de la propuesta, se puede
aceptar al N.S. 𝛼 = 0,05 que los resultados de ambas encuestas son similares?
Poblaciones no normales
Pruebas No Paramétricas
Prueba 𝝌𝟐 de bondad del ajuste
Comprobar si el número de observaciones en una serie de grupos (o intervalos)
se corresponde con uno dado.
Si tenemos una serie de “p” clases con un número de observaciones “oi”
(i=1,2…p) en cada una y queremos comparar esos valores con los conteos
esperados “ei”. Se tiene:
𝑝 2
𝑜𝑖 − 𝑒𝑖 2
~𝜒𝑝−1
𝑒𝑖
𝑖=1
Lo que permite establecer con qué P son los oi = ei en la población de la que
proceden las muestras
𝑝 𝑜𝑖 −𝑒𝑖 2 2
En general, el valor calculado 𝜒𝑐2 = σ𝑖=1 se compara con 𝜒𝑝−1−𝜁 donde
𝑒𝑖
“p” es el número de clases y “𝜁” el número de parámetros calculados para
calcular los valores esperados “ei”
En el caso de tablas de doble entrada con “p” filas y “q” columnas” es
𝜁 = 𝑝−1 · 𝑞−1
Poblaciones no normales
Pruebas No Paramétricas
Prueba 𝝌𝟐 de bondad del ajuste
Si tenemos una serie de “p” clases con un número de observaciones “oi”
(i=1,2…p) en cada una y queremos comparar esos valores con los conteos
esperados “ei”.
Fijamos un N.S. 𝛼. Puesto que:
𝑝 2
𝑜𝑖 − 𝑒𝑖 2
𝜒𝑐2 = ~ 𝜒𝑝−1−𝜁
𝑒𝑖
𝑖=1
1º H0: 𝑜𝑖 = 𝑒𝑖
HA: 𝑜𝑖 ≠ 𝑒𝑖
𝑝 𝑜𝑖 −𝑒𝑖 2
2º Calculamos 𝜒𝑐2 = σ𝑖=1
𝑒𝑖
2
𝑆𝑖 𝜒𝑐2 ≤ 𝜒𝑛−1;1−𝛼 ֜ 𝐴𝑐𝑒𝑝𝑡𝑎𝑚𝑜𝑠 𝐻𝑜
3º Decisión ቐ 2
𝑆𝑖 𝜒𝑐2 > 𝜒𝑛−1;1−𝛼 ֜ 𝑅𝑒𝑐ℎ𝑎𝑧𝑎𝑚𝑜𝑠 𝐻𝑜
Recordar que las Tablas de la distribución 𝜒 2 suelen dar el valor que deja una probabilidad α a su
2
derecha, en lugar de a su izquierda: 𝑃 𝜒 2 > 𝜒𝑛−1; 𝛼 =𝛼
Poblaciones no normales
Construimos la Tabla:
Grupo 1 Grupo 2 Grupo 3 Total 3
2
1º. 𝜒𝑐2 ~ 𝜒3−1 ; 𝑝𝑜𝑟𝑞𝑢𝑒 𝜁 = 0 [Como las proporciones esperadas “𝑝𝑖 ” vienen dadas, no ha sido
necesario calcular ningún parámetro para obtenerlas]
𝐻0 : 𝑜𝑖 = 𝑒𝑖 𝑦 𝐻𝐴 : 𝑜𝑖 ≠ 𝑒𝑖
𝑝 𝑜𝑖 −𝑒𝑖 2 83−80 2 42−40 2 35−40 2 67
2º. 𝜒𝑐2 = σ𝑖=1 = + + = = 0.8375
𝑒𝑖 80 40 40 80
2 2 2
3º. a) 𝜒3−1−0;1−0.05 = 𝜒2; 0.95 = 5.99 ; 𝜒𝑐2 = 0.8375 ≪ 5.99 = 𝜒2; 0.95
֜ 𝑁𝑜 𝑟𝑒𝑐ℎ𝑎𝑧𝑎𝑚𝑜𝑠 𝐻0 𝑎𝑙 𝑁. 𝑆. 𝛼 = 0.05
b) 𝑃 𝜒22 ≤ 𝜒𝑐2 = 𝑑𝑖𝑠𝑡𝑟. 𝑐ℎ𝑖𝑐𝑢𝑎𝑑 0.8375; 2; 1 = 0.3421. La P de obtener un valor de 𝜒22 >
0.8375 es de 0.6579. [También se puede usar directamente la función distr.chicuad.cd(0.8375; 2)].