Pruebas de Hipotesis
Marcelo Guaman Hidalgo
13 de junio de 2017
1. Prueba de hipotesis de compracion de dos medias con σ conocida
1.1. Prueba de hipotesis comparacion de dos medias con σ conocida
Ejercicio 1
Un profesor realiza un estudio sobre el rendimiento en la clase de estadistica y sus desempeño academico.
Tomo una muestra de 32 estudiantes de primer nivel y obtuvo una media de 13.75 y una media poblacional de
2,50. (De informacion historica y otra muestra de 28 observaciones cuya media fue de 14.50 con igual desviacion
etadar poblacional que le grupo anterior. Realizar una prueba de hipotesis que verifique si el primer grupo tiene
una media menor que el segundo grupo. Considerando un nivel de significancia de 0.02.
1. Planteamiento de la hipotesis
H0 : µ1 = µ2
H1 : µ1 < µ2
2. nivel de significancia
α = 0, 02
3. Calculo del estadistico de prueba
x¯1 − x¯2
z=q 2
σ1 σ22
n1 + n2
13, 75 − 14, 50
z= q
22 2,502
32 + 28
z = −1, 27
4. Determinacion del estadistico de contraste
Para determinar el estadistico de contraste emplearemos el uso de R para lo cual utilizaremos el comando
qnorm.
> q1=round(qnorm(0.02, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE),3)
> q1
[1] -2.054
>
El valor del estadistico critico de la izquierda es zc = −2,054
5. Toma de decision
El estadistico de prueba z, se encuentra entre los valores de zc , es decir cae en la zona de no rechazo. Por
lo cual no rechazo H0 . No tengo evidencia suficiente para rechazar la hipotesis nula.
1
6. Solucion usando el software R
Para resolver el ejemplo usando R, primero debemos cargar la libreria PASWR.
> library(PASWR)
Una vez que tenemos cargada esta libreria usamos el comando z. test, de la siguiente manera, pero
trataremos de emplearlo haciendo una pequena simulacion:
> set.seed(10)
> x<- rnorm(32,13.75,2.5)
> x
[1] 13.796865 13.289369 10.321674 12.252081 14.486363 14.724486 10.729810
[8] 12.840810 9.683318 13.108804 16.504449 15.639454 13.154416 16.218612
[15] 15.603475 13.973368 11.362640 13.262124 16.063803 14.957446 12.259223
[22] 8.286783 12.062835 8.452347 10.587005 12.815846 12.031111 11.569603
[29] 13.495597 13.115549 9.115649 13.555135
> y<-rnorm(28,14.50,2.5)
> y
[1] 16.92142 14.96231 11.05014 10.91121 15.40522 10.10228 13.68864 12.87109
[9] 17.21638 12.59364 12.42834 16.58618 12.08087 14.42796 15.08131 13.74698
[17] 12.80596 16.13807 13.49841 13.66361 17.91988 19.84442 15.76455 16.46586
[25] 12.24447 15.83224 12.88526 15.22747
> z.test(x, y, alternative = "less", mu = 0, sigma.x = 2.5,
+ sigma.y = 2.5, conf.level = 0.98)
Two-sample z-Test
data: x and y
z = -2.4406, p-value = 0.007332
alternative hypothesis: true difference in means is less than 0
98 percent confidence interval:
-Inf -0.2502529
sample estimates:
mean of x mean of y
12.79125 14.37015
El ministerio de agricultura a traves del INIAP ha obtenido una nueva variedad de maiz la cual quiere
contrastar los niveles de produccion con la variedad que actualmente se usa. Para lo cual realizo un experimento,
donde la variedad anterior sembró en 16 fincas en la que se obtuvo un rendimiento promedio de 75qq por
hectarea, mientras que la variedad de maiz nueva se sembaron en 25 fincas obteniendose un rendimiento de
85qq. La desviación estandar obtenida de la muestra 1 fue de 7qq por hectárea, mientras que la nueva variedad
presenta una desviacion estandar de 8qq. Realizar una prueba de hipótesis que verifique que la nueva variedad
tiene mayor rendimiento. Considere un nivel de significancia de α = 0,01 .
1. Planteamiento de la hipotesis
H0 : µ1 = µ2
H1 : µ1 < µ2
2. nivel de significancia
α = 0, 01
3. Calculo del estadistico de prueba
2
x¯1 − x¯2
z=q 2
σ1 σ22
n1 + n2
75 − 85
z=q
52 52
16 + 25
z = −6, 25
4. Determinacion del estadistico de contraste
Para determinar el estadistico de contraste emplearemos el uso de R para lo cual utilizaremos el comando
qnorm.
> q_t1=qt(p=0.01, df=39, lower.tail = TRUE, log.p = FALSE)
> q_t1
[1] -2.425841
El valor del estadistico critico de la izquierda es zc = −2,054
5. Toma de decision
El estadistico de prueba z, cae en la zona de aceptación por tanto rechazo
6. Calculo del valor-p
z z2
e− 2
Z
valor − p = · √ dz
−∞ 2π
−2,054 z2
e− 2
Z
valor − p = · √ dz
−∞ 2π
valor − p = 0,0478
V alor − p > α por lo tanto no rechazo H0
7. Solucion usando el software R
Para resolver el ejemplo usando R, primero debemos cargar la libreria PASWR.
> library(PASWR)
Una vez que tenemos cargada esta libreria usamos el comando t. test, de la siguiente manera, pero trata-
remos de emplearlo haciendo una pequena simulacion:
> set.seed(20)
> x<- rnorm(16,75,7)
> x
[1] 83.13880 70.89853 87.49826 65.67184 71.87403 78.98724 54.77198 68.91687
[9] 71.76808 71.11121 74.85905 73.94732 70.60311 84.26255 64.35055 71.93800
> y<-rnorm(25,85,8)
> y
[1] 92.76462 85.22578 84.31374 88.11372 86.89350 83.84448 90.77784 87.95925
[9] 83.06347 73.22349 80.23072 75.82640 65.20291 80.09193 83.26951 97.72117
[17] 97.44915 93.86761 76.22127 70.11515 77.69137 94.96455 85.70284 88.38786
[25] 78.45214
> t.test(x,y,alternative ="less", mu=0, paired = FALSE, var.equal = FALSE, conf.level = 0,99 )
3
Welch Two Sample t-test
data: x and y
t = -4.3097, df = 32.906, p-value = 6.96e-05
alternative hypothesis: true difference in means is less than 0
0 percent confidence interval:
-Inf -Inf
sample estimates:
mean of x mean of y
72.78734 84.05498
