TEST DE BONDAD DE

AJUSTE: EJERCICIOS
HIDROLOGÍA CIV 231
DOCENTE: MS. C. ING. CÉSAR LUIS VISCARRA PINTO
EJERCICIO Nº 1:
CHI CUADRADO Y SMIRNOV KOLMOGOROV
Dada la serie histórica de caudales medios anuales en m3/s, que
corresponde a un registro de 38 años. Realizar el test de bondad de
ajuste Chi Cuadrado y Smirnov Kolmogorov, para ver si se ajustan a
una distribución normal.

121,3 26,7 110,1 63,4 122,4 64,2 59,6

144,9 92,8 95,6 76,3 162,1 110,2 40,3
142,4 58,8 48,8 52,3 97,2 144,7 112,2
205,8 57,4 148,3 36,3 52,5 109,2 137,1
114,5 79,0 67,5 88,0 165,6 48,5 32,9
72,5 76,9 70,0

C.L.V.P.
SOLUCIÓN: Chi cuadrado
1. La hipótesis será:
Ho: frecuencia observada = frecuencia esperada
Ha: frecuencia observada = frecuencia esperada
2. Ordenando los datos de menor a mayor, se tiene:
Ordenando los datos de menor a mayor, se tiene:
26,7 32,9 36,3 40,3 48,5 48,8 52,3
52,5 57,4 58,8 59,6 63,4 64,2 67,5
70,0 72,5 76,3 76,9 79,0 88,0 92,8
95,6 97,2 109,2 110,1 110,2 112,2 114,5
121,3 122,4 137,1 142,4 144,7 144,9 148,3
162,1 165,5 205,8
C.L.V.P.
SOLUCIÓN: Chi cuadrado
3. Cálculo de la distribución de frecuencias
3.1. Cálculo del número de intervalos, según Yevjevich:
NC  1 1,33 ln(N) NC  1 1,33 ln(38) NC  5,84  6
3.2. Cálculo de la amplitud de cada intervalo:

Xmax  Xmin 205,8  26,5 x

x  x  x  35,82  36  18
NC  1 6 1 2

3.3. Cálculo de los intervalos de clase, marcas de clase, frecuencia absoluta

observada, frecuencia relativa.

C.L.V.P.
SOLUCIÓN: Chi cuadrado
3.3. Tabla de distribución de frecuencias
Marca Frecuencia Frecuencia Frecuencia
Intervalos
de clase absoluta relativa acumulada
(1)
(2) (3) (4) (5)
8,7 – 44,7 26,7 4 0,1053 0,1053
44,7 – 80,7 62,7 15 0,3947 0,5000
80,7 – 116,7 98,7 9 0,2368 0,7368
116,7 – 152,7 134,7 7 0,1842 0,9211
152,7 – 188,7 170,7 2 0,0526 0,9737
188,7 – 224,7 206,7 1 0,0263 1,0000

C.L.V.P.
SOLUCIÓN: Chi cuadrado
3.4. Cálculo de la media y la desviación estándar para datos
agrupados.
k k 2
 _

_  x ifi   x i  x  fi
x i 1
 90,17 i 1    43,03
S
N N 1
4. Cálculo de la frecuencia esperada, utilizando la distribución
teórica normal. Para el calculo se realizaron los siguientes cálculos:
_

Columna (3): xx variable estandarizada de la distribución

Z
normal para S x = limites de cada intervalo, de la
columna (2)

C.L.V.P.
SOLUCIÓN: Chi cuadrado
4. Cálculo de la frecuencia esperada, utilizando la distribución
teórica normal. Para el calculo se realizaron los siguientes cálculos:
Columna (4): área bajo la curva normal
Columna (5): área para cada intervalo de clase, se obtiene restando
los valores de la columna (4), si los signos de Z de la columna (3) son
iguales, y sumando los valores de la columna (4), si los signos de Z
son diferentes.
Columna (6): Columna (5) x N= 38, se redondea en forma adecuada,
de tal manera que la suma de las frecuencias absolutas sea igual a N
= 38.
Columna (7): son los mismos valores de la columna (3) de la tabla
anterior.

C.L.V.P.
 SOLUCIÓN: Chi cuadrado
4. Cálculo de la frecuencia esperada, utilizando la distribución
teórica normal.
_ Área bajo la
Limite Z  x  x curva Frecuencia
Intervalos S Frecuencia Frecuencia
de clase normal de relativa
(1) Absoluta Observada
(2) (3) 0aZ (5)
(4)
8,7 -1,89 0,4706 -
8,7 – 44,7 44,7 -1,057 0,3554 0,1152 4,38 = 5 4
44,7 – 80,7 80,7 -0,220 0,0871 0,2687 10,21 = 10 15
80,7 – 116,7 116,7 0.617 0,2324 0,3195 12,14 = 12 9
116,7 – 152,7 152,7 1,453 0,4265 0,1931 7,34 =7 7
152,7 – 188,7 188,7 2,290 0,4890 0,0625 2,38 = 3 2
188,7 – 224,7 224,7 3,126 0,4991 0,0101 0,38 = 1 1
C.L.V.P.
SOLUCIÓN: Chi cuadrado
5. Cálculo de Xc2
Se conoce:  
2
k
i  ei 2
c
i 1 ei
Sustituyendo valores:

2 
4  5 2


15  10 2


9  12 2


7  7 2


2  3 2


1 12
 c2  3,78
c
5 10 12 7 3 1
6. Calculo del Xt2
Grados de libertad: v=k–1–h
v=6–1–2=3
Nivel de significancia: a = 0,05 = 5%

C.L.V.P.
SOLUCIÓN: Chi cuadrado
Con los valores de: v = 3 y a = 0,05; se tiene:  t2  7,81

7. Los criterios de decisión son:

 c2   t2 se acepta la hipótesis
 c2   t2 se rechaza la hipótesis
en nuestro caso se tiene:  c2  3,78   t2  7,81

Los datos se ajustan a la distribución normal, con un nivel

de significación del 5% ó 95% de probabilidad.

C.L.V.P.
SOLUCIÓN: Smirnov Kolmogorov
1. Cálculo de P(x): Ordenando los datos de caudales en forma creciente y
calculando la probabilidad empírica P(x), usando la formula de Weibull:
m
P(x) 
n 1
Se obtienen las columnas (2) y (3) de la tabla siguiente.
2. Calculo de x y S, de datos no agrupados.

2
_
1 1  _

x   x i  92,32 S  i   42,80
 x  x
_
n n 1 
xx
3. Calculo de la variable estandarizada Z; usando la ecuación: Z
S
Con dicha ecuación se obtiene la columna (4)

C.L.V.P.
SOLUCIÓN: Smirnov Kolmogorov
4. Calculo de F(Z) = F(X), usando la tabla de distribución normal se
obtiene la columna (5). Para los valores positivos de Z, los valores se
obtienen en forma directa. Para valores negativos de Z, los resultados se
obtienen como 1 – valor de la tabla.
5. Calculo de  = |F(Z) - P(X)|
A partir de las columnas (3) y (5) de la tabla de resultados se obtiene los  = |F(Z)
- P(X)|, la misma que se muestra en la columna (6)
6. Cálculo del  = max:
De la tabla de cálculos se observa que:  = max |F(z) – P(x)| = 0,1089

C.L.V.P.
SOLUCIÓN: Smirnov Kolmogorov
7. Calculo de o critico; de la tabla del test de Smirnov Kolmogorov se tiene:

1,36 1,36
0    0  0,22
n 38
8. Criterio de decisión

  0,1089   0  0,22

Se concluye que los datos de caudales se ajustan a la distribución

normal, con un nivel de significancia del 5% o una probabilidad del
95%.

C.L.V.P.
SOLUCIÓN: Smirnov Kolmogorov
_
Q=X P(x) = xx F(z)
m Z =|F(z)-=(x)|
m3/s m/(n+1) S
(1) (6)
(2) (3) (4) (5)
1 26,7 0,0256 -1,53 0,0630 0,0374
2 32,9 0,0513 -1,39 0,0823 0,0310
3 36,3 0,0769 -1,31 0,0951 0,0182
4 40,3 0,1026 -1,22 0,1112 0,0086
5 48,5 0,1282 -1,02 0,1539 0,0257
6 48,8 0,1538 -1,02 0,1539 0,0001
7 52,3 0,1795 -0,94 0,1736 0,0059
8 52,5 0,2051 -0,93 0,1762 0,0289
9 57,4 0,2308 -0,82 0,2061 0,0247
10 58,8 0,2564 -0,78 0,2177 0,0387
11 59,6 0,2821 -0,76 0,2236 0,0585
12 63,4 0,3077 -0,68 0,2486 0,0594
13 64,2 0,3333 -0,66 0,2546 0,0787
C.L.V.P.
SOLUCIÓN: Smirnov Kolmogorov
_
Q=X P(x) = xx F(z)
m Z =|F(z)-=(x)|
m3/s m/(n+1) S
(1) (6)
(2) (3) (4) (5)
14 67,5 0,3590 -0,58 0,2810 0,0780
15 70,0 0,3846 -0,52 0,3015 0,0831
16 72,5 0,4103 -0,46 0,3228 0,0875
17 76,3 0,4359 -0,37 0,3557 0,0802
18 76,9 0,4617 -0,36 0,3594 0,1021
19 79,0 0,4872 -0,31 0,3783 0,1089
20 88,0 0,5128 -0,10 0,4602 0,0526
21 92,8 0,5385 0,01 0,5040 0,0345
22 95,6 0,5641 0,08 0,5319 0,022
23 97,2 0,5897 0,11 0,5438 0,0459
24 109,2 0,6154 0,39 0,6517 0,0363
25 110,1 0,6410 0,42 0,6628 0,0218
26 110,2 0,6667 0,42 0,6628 0,0039
C.L.V.P.
SOLUCIÓN: Smirnov Kolmogorov
_
Q=X P(x) = xx F(z)
m Z =|F(z)-=(x)|
m3/s m/(n+1) S
(1) (6)
(2) (3) (4) (5)
27 112,2 0,6923 0,46 0,6772 0,0151
28 114,5 0,7179 0,52 0,6985 0,0194
29 121,3 0,7436 0,60 0,7517 0,0081
30 122,4 0,7692 0,70 0,7580 0,0112
31 137,1 0,7949 1,05 0,8531 0,0582
32 142,4 0,8205 1,17 0,8790 0,0585
33 144,7 0,8462 1,22 0,8888 0,0424
34 144,9 0,8718 1,23 0,8907 0,0189
35 148,3 0,8974 1,31 0,9084 0,0075
36 162,1 0,9231 1,63 0,9484 0,0253
37 165,5 0,9487 1,71 0,9564 0,0077
38 205,8 0,9744 2,65 0,9960 0,0216

C.L.V.P.