108
La cuenca del río Usumacinta
desde la perspectiva del cambio climático
esta figura, a diferencia del gráfico de la derecha de la Fig. 4, se
muestran discontinuidades en cada punto donde se registra una
precipitación diferente de los datos, lo que les da a las funciones
una apariencia escalonada. El procedimiento, a partir de esta
familia de funciones, para mapear las precipitaciones a una
distribución normal estandarizada es semejante al descrito en la
sección anterior. La serie de SPIs resultante de utilizar funciones
ECDF en vez de funciones de distribución Gamma, se muestra en
el gráfico inferior de la Fig. 6.
Es importante destacar que la implementación en el lenguaje R,
desarrollada por Wheatley (2010), deriva de la implementación
alternativa de Edwards y McKee (1977) que toma como base
las funciones ECDF, y que la implementación desarrollada
para el presente trabajo, a su vez, se fundamenta en esas dos.
Lo anterior obedece a la ventaja que tienen estas funciones al
trabajar directamente con los datos, en vez de hacerlo con una
aproximación de ellos. Una desventaja del uso de las funciones
ECDF es que, en ellas, todas las precipitaciones máximas arrojan
una probabilidad de 1 y, por consiguiente, éstos valores al ser
mapeados mediante la inversa de la función normal estandarizada
a sus correspondientes valores de SPI, dan un valor infinito.
Para lidiar con este problema, se recurre a un procedimiento
de optimización que reajusta, para esos valores de precipitación,
valores de SPI que mantienen la distribución normal estandarizada,
esto es, con media 0 y desviación estándar 1 (Wheatley, 2010).
