marzo-08

Distribuciones Estadísticas

Valores extremos. Distribución de Gumbel

Esta ley de distribución de frecuencias se utiliza para el estudio de los valores extremos. Por
ejemplo, si hemos elegido el día mas caudaloso o de mayor precipitación de cada año de una serie
de años.
La probabilidad de que se presente un valor inferior a x es:

b= (x – u)
b
F (x)  ee (1) siendo: (2)

 =  y / sx (3)
y 
ux (4)

e = base de los logaritmos neperianos
x = media aritmética de la muestra
sx = desviación típica de la muestra
y , y = consultar en la tabla adjunta, según
el número de datos de la muestra2] nº datos y y
10 0,4952 0,9496
(Ver nota3 ) 15 0,5128 1,0206
Mediante las expresión anteriores podremos calcular 20 0,5236 1,0628
la frecuencia a partir del valor x, es decir: calcularcon 25 0,5309 1,0914
qué frecuencia (o periodo de retorno) se presentará un 30 0,5362 1,1124
cierto caudal o precipitación. 35 0,5403 1,1285
40 0,5436 1,1413
Para solucionar el caso inverso (qué caudal o 45 0,5463 1,1518
precipitación se producirán cada n años) debemos 50 0,5485 1,1607
despejar b en la expresión (1), obteniendo: 55 0,5504 1,1682
60 0,5521 1,1747
b = –ln (–ln (F(x))) (5) 65 0,5535 1,1803
70 0,5548 1,1854
Y, finalmente, despejando x en (2): 75 0,5559 1,1898
x =(b/  + u (6) 80
85
0,5569
0,5578
1,1938
1,1974
90 0,5586 1,2007
Ejemplo.- De una serie de 55 caudales extremos (el 95 0,5593 1,2037
caudal diario máximo de cada año)4, hemos calculado: 100 0,5600 1,2065

2 y , yson, respectivamente, la media y la desviación típica de una serie de


valores yi (i = 1 a N ; N = nº de datos de la muestra) que dependen solamente
del número de datos, y que corresponden a la siguiente expresión:
(Cálculo de estos parámetros en un documento Excel, en la web,
sección “Complementos”)
Wanielista (1997) utiliza y =0,5772; y = 1,2825 sin considerar el número de
datos

Pág. 1
 Media= 21,97 m3/seg
Desv típica=13,22 m3/seg
Calcular: a) Probabilidad de que se supere un caudal de 60 m3/s. b) Qué caudal se superará el
1% de los años

a) ¿Cual será la probabilidad de que (el día más caudaloso del año) el caudal supere el valor
de 60 m3/seg?
1º) De acuerdo con la tabla adjunta, para 55 datos, tomamos los valores:
y = 0,5504 y = 1,1682

2º) Calculamos  y u:
   y / sx = 1,1682 / 13,22 = 0,0884
u = x - y /    – 0,5504 / 0,0884 = 15,741

3º) Calculamos el exponente b:

b=  (x–u) = 0,0884 · (60 – 15,741) =3,917

4º) Aplicamos la ecuación de Gumbel (1) para el caudal del problema (60 m3/s). La probabilidad
de que se presente un caso menor que x será:
3, 917
F (x) e e  0,9803  98,03%

Por tanto, la probabilidad de que se presente un caso mayor que x será:

1- F (x) = 1 -0,9803 = 0,0197 (= 1,97%)
Finalmente, el periodo de retorno es el inverso de la probabilidad:
Periodo de retorno= 1/0,0197 = 50,8 años

b) Caso inverso: Calcular qué caudal se superará un 1% de los casos.

Si un caudal se superará el 1% de los años, será inferior el 99%, es decir: F(x) =0,99.
Calculamos  y u mediante las exprresiones (3) y (4). En este caso, ya las hemos calculado en
el apartado anterior.

Aplicando la fórmula (5):

b = –ln (–ln (F(x))) = –ln (–ln (0,99)) = 4,600

Aplicando la fórmula (6):

x =(b /  + u = 4,600 / 0,0884 + 15,741 = 67,8 m3/s

La simplificación mostrada a lo largo de todas estas páginas (indicando que los valores medios se ajustan
a Gauss, y los valores extremos se ajustan a Gumbel) es solamente válida con fines didácticos, para una
primera aproximación al tema.
Existen muchas otras distribuciones, entre las que destacan, como más utilizadas, la lognormal (los
logaritmos de los valores son los que se ajustan a la ley de Gauss) o la ley Pearson III, adoptada por las
agencias federales en USA. Ver, por ejemplo en Viessman, 2003, capítulo 3.
En España los organismos oficiales para precipitaciones máximas aplican la distribución SQRT-max5

Pág. 2
Probabilidad, periodo de retorno y riesgo de fallo
A lo largo de los apartados anteriores se ha estado utilizando indistintamente probabilidad (por
ejemplo: un 2% de los años) y expresiones como “cada 50 años”.
Es evidente que si un suceso se presenta (por término medio) cada 10 años, su probabilidad es
de 0,10 (10%).
Análoga e inversamente, si la probabilidad de que algo suceda es de 0,04 (4%), ello quiere
decir que, en promedio, sucederá 4 veces en 100 años, o sea: una vez cada 25 años.
Estos conceptos se relacionan mediante la expresión:
1
Periodo de retorno 
Probabilidad
Chow et al. (1994, p.393) ofrecen una elaborada demostración de la fórmula anterior (¡ !).
En Hidrología se utiliza más el periodo de retorno que la probabilidad. Así, se habla de la
crecida de 50 años en lugar de referirse a la crecida con probabilidad 0,02 o de la precipitación
con retorno de 100 años en vez de la precipitación con probabilidad 0,01.
Supongamos que hemos calculado un cierto caudal que corresponde al retorno de 50 años. La
probabilidad de que se produzca el año próximo será de 0,02 (=1/50); y la probabilidad de que se
produzca el siguiente año será de 0,02 y así cada año. Necesitamos conocer la probabilidad de que
se alcance ese caudal en los próximos n años:

Probabilidad de que un suceso de retorno T se produza el próximo año ...........................1/T

“ “ “ NO se produzca el próximo año (*) ................1-(1/T)
“ “ “ NO se produzca los próximos dos años(**) ....[1-(1/T)]. [1-
(1/T)]
“ “ “ NO se produzca los próximos n años..............[1-(1/T)]n
“ “ “ SI se produzca los próximos n años (*)........1-[1-(1/T)]n

Vamos a denominar a la última expresión obtenida arriba es el riesgo de fallo (R), es decir: la
probabilidad de que sí se produzca alguna vez un suceso de periodo de retorno T a lo largo de un
periodo de n años:
n
 1
R  1 1 
 T 

Ejemplo: Se va a construir un canal cuya vida útil es de 75 años. Si el caudal supera el valor
correspondiente al periodo de retorno de 100 años, se desbordará. Calcular la probabilidad de que
se produzca un desbordamiento en alguno de los próximos 75 años
 1 
R  1  1  75  0,529  52,9%


 100 
Por tanto, existe un 52,9% de probabilidad de que el caudal de retorno 100 años se alcance en
alguno de los próximos 75 años.
Se produce la siguiente paradoja: si consideramos un caudal con retorno de 100 años, parece
seguro que se presente en alguno de los próximos 100 años. Pero si aplicamos la fórmula anterior,
haciendo T= 100 y n=100, y obtenemos 0,633 , es decir solamente un 63,3 %

probabilidad de obtener un 3 en un dado= 1/6. Probabilidad de obtener dos 3 seguidos = 1/6.1/6 =1/36
Pág. 3
 Cálculo inverso: evaluación del periodo de retorno a partir del riesgo
Muchas veces, esta sencilla fórmula debe aplicarse a la inversa: despejar T a partir del riesgo de
fallo (R) y del número de años (n). Depejando T en la última fórmula obtenemos:
1
T 
ln(1  R) 
 1  exp
 
 n 
donde: R =riesgo de que se produzca el suceso de probabilidad 1/T durante los próximos n años
T = periodo de retorno en años
exp( x) = e x
Ejemplo: se está diseñando una obra cuya vida útil se calcula en 50 años y se admite que en
ese periodo el riesgo sea de un 10% (probabilidad de que en esos 50 años se produzca un caudal
superior a un valor determinado). Calcular dicho caudal.
En la fórmula anterior basta con hacer: R = 0,10; n = 50 años; y despejar T. Con estos datos
obtenemos un periodo de retorno T = 475 años.
En este ejemplo, el paso siguiente sería estudiar estadísticamente las series históricas de
caudales de ese cauce para evaluar el caudal correspondiente a un retorno de 475 años
(Probabilidad = 1/475 = 0,0021).

Pág. 4
Pág. 5