MANUAL DE USUARIO FINAL PROGRAMA DE FUNCIONES DE DISTRIBUCION Ax+B VERSIÓN 1.0.0.0 Autor del manual: Ing. Fidelmar Merlos Villegas © Revisión del Manual: Ing. Jorge Alfonso Almanza Campos Febrero 2015 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 © Fidelmar Merlos Villegas ISBN: En tramite P á g i n a 2 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 DESCARGO DE RESPONSABLIDAD La información contenida en este documento ha sido resultado de un amplio proceso de investigación y recopilación de datos. Esta información y sus resultados obtenidos han sido sometidos a una revisión exhaustiva por parte del autor. Las menciones realizadas a marcas o productos comerciales no constituyen un reconocimiento o recomendación para su uso. Aunque se ha realizado un esfuerzo considerable para garantizar que los resultados sean correctos, el/los programas descritos en este manual son susceptibles a depuración. Por consiguiente ni el autor ni sus colaboradores se hacen responsables ni asumen ninguna relación con los resultados obtenidos con el programa, del uso de los mismos, ni de daños o litigios que resulten de la utilización de este manual o del programa referido. P á g i n a 3 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 PROLOGO El Programa de Funciones de Distribución Ax+B se ha creado con el objetivo de ser de ayuda en la modelación matemática de eventos extremos y con ello ser de colaborar en la investigación de los mismos. En esta primera revisión de la versión del software se incluyen diversas funciones de distribución, entre las cuales se incluyen la Doble Gumbel y Doble LogNornal. Teniendo para cada función de distribución diferente formas de estimar los parámetros de las mismas. P á g i n a 4 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 Contenido DESCARGO DE RESPONSABLIDAD .................................................................. 3 PROLOGO ............................................................................................. 4 1.- INTRODUCCIÓN .................................................................................... 6 2.- MARCO TEÓRICO DEL SOFTWARE................................................................ 7 2.1.- PERIODO DE RETORNO ....................................................................... 7 2.2.- FUNCIONES DE DISTRIBUCIÓN .............................................................. 7 2.3.- TEORÍA DE FUNCIONES DE DISTRIBUCIÓN DE DOS POBLACIONES ...................... 9 2.4.- ESTIMACIÓN DE PARÁMETROS .............................................................. 9 2.4.1.- METODO DE MOMENTOS ............................................................... 11 2.4.2.- METODO DE MAXIMA VEROSIMILITUD ................................................ 12 2.4.3.- METODO DE MOMENTOS L ............................................................. 12 2.4.4.- MAXIMA ENTROPIA ..................................................................... 13 2.5.- PRUEBAS DE BONDAD DE AJUSTE ......................................................... 14 2.5.1.-ERROR CUÁDRATICO .................................................................... 14 2.5.2.- KOLMOGOROV-SMIRNOF ............................................................... 15 2.5.3.- CHI CUADRADO ......................................................................... 15 2.5.4.- ANDERSON DARLING ................................................................... 15 2.5.5.- NASH-SUTOLFFE ........................................................................ 15 3.- UTILIZACIÓN DEL PROGRAMA ................................................................... 16 3.1.- REQUISITOS MINIMOS ....................................................................... 16 3.2.- INTERFAZ ..................................................................................... 16 3.2.- BARRA DE MENÚ ............................................................................. 17 3.3.- BARRA DE HERRAMIENTAS .................................................................. 17 3.4.- CONFIGURACIÓN DEL PROGRAMA ......................................................... 18 3.5.-ENTRADA DE DATOS .......................................................................... 19 3.6.- OPTIMIZACIÓN DE FUNCIONES DE DISTRIBUCIÓN DE DOS POBLACIONES ............ 20 3.7.- GRÁFICO DE FUNCIONES DE DISTRIBUCIÓN .............................................. 21 3.8.- EXTRAPOLAR DATOS ........................................................................ 22 3.9 GENERACION DE REPORTES .................................................................. 23 4.- BIBLIOGRAFIA ..................................................................................... 25 P á g i n a 5 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 1.- INTRODUCCIÓN Este manual y el software referido representan un esfuerzo por facilitar la tarea del hidrólogo y de cualquier otro profesional de las ciencias de la tierra que requiere de la aplicación de los procedimientos aquí descritos. La intención de la autor de este software es facilitar la investigación y desarrollo de las ciencias, es específico el desarrollo de metodologías para la modelación de eventos extremos de carácter hidrológico. Cabe mencionar que el programa se puede ver como parte medular del conjunto de software para hidrología encabezado por el programa SIHIMax. Si requiere información adicional a este manual puede consultar la página: http://hydrobits.com/AxmasB/. P á g i n a 6 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 2.- MARCO TEÓRICO DEL SOFTWARE Para el desarrollo de este software se ha efectuado una revisión exhaustiva de las metodologías existentes para estimar los parámetros de cada una de las funciones de distribución propuestas. Parte de los resultados de esta revisión se presentan en este apartado. 2.1.- PERIODO DE RETORNO El periodo de retorno se define como el número de años promedio en el cual el evento puede ser igualado o excedido cuando menos una vez. Existen diversos autores que han propuesto ecuaciones para estimar el periodo de retorno, el programa de funciones de distribución Ax+B emplea la ley Weibull: ππ = π+1 π Donde: Tr = Periodo de retorno n = Numero de datos de la serie m = Numero de orden de la serie 2.2.- FUNCIONES DE DISTRIBUCIÓN El programa Ax+B es capaz de realizar el ajuste funciones de distribución de probabilidad a series de datos, esto mediante la aplicación de diversos procesos para estimar los parámetros de cada una como se menciona en el apartado 2.3. Las funciones de distribución aplicadas son: Nombre de la función Ecuación π§ Normal Log-Normal de 2 parámetros Log-Normal de 3 parámetros πΉ(π‘) = β« π§2 1 ββ β2π π§ πΉ(π‘) = β« 0 π§ πΉ(π‘) = β« 0 β π β 2 β π§ 1 β2π 1 π§2 π§= β π β 2 β π§ π§ = π₯βπ π ln(π₯) β π π π§2 β π β 2 β π§π§ = ln(π₯ β π₯π) β π β2π π P á g i n a 7 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 πΉ(π₯) = π βπ Gumbel ββ(π₯βπ½) Exponencial πΉ(π₯) = π½π βπ½π₯ Exponencial de dos parámetros πΉ(π₯) = 1 β π (π₯βπ)/π Gamma de dos parámetros Gamma de tres parámetros Doble Gumbel πΉ(π₯) = πΉ(π₯) = π₯ π₯ β« ( )πβ1 π (π₯)/π β π₯ πΞ (π) πΏ πΌ 1 π₯ 1 πΞ (π) πΉ(π₯) = ππ βπ β« ( βπΌ1 (π₯βπ1 ) π§1 πΉ(π₯) = π β« Log-Normal Doble 0 (1 β π) β« Frechet + (1 β π)π βπ 1 β2π π§2 1 0 Doble Frechet πΏ π₯ β πΏ πβ1 (π₯βπΏ)/π ) π β π₯ πΌ β2π β πβ β π§ 2 β 2 π 2 πΉ(π₯) = π βπ πΉ(π₯) = ππ βπ βπΌ1 (πΏπ(π₯)βπ1 ) π§1 2 2 βπΌ2 (π₯βπ2 ) β π§ + β π§ ββ(lnβ‘(π₯)βπ½) + (1 β π)π βπ βπΌ2 (πΏπ(π₯)βπ2 ) TABLA 1.- F UNCIONES DE DISTRIBUCIÓN APLICADAS EN EL PROGRAMA AX+B P á g i n a 8 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 2.3.- TEORÍA DE FUNCIONES DE DISTRIBUCIÓN DE DOS POBLACIONES Es frecuente que las series temporales en hidrología sean de dos poblaciones, esto debido a la existencia de fenómenos ciclónicos que alteran estadísticamente la serie temporal. Por lo que se forman dos muestras ordenadas de datos, como se muestra a continuación: X1 X2 X3 X4 X5 .. .. .. .. Xn Población Ciclónica Población No Ciclónica De manera que estadísticamente tenemos: πΉ(π₯) = ππΉ1 (π₯) + (1 β π)πΉ2 (π₯) Donde: p = Cociente del número de datos de la población No Ciclónica y la muestra total. F1(x) = Probabilidad de la función No Ciclónica F2(x) = Probabilidad de la función Ciclónica Cabe mencionar que en la funcionalidad del programa el valor de p introducido esta en porcentaje y se ajusta de manera automática al número de datos reales. 2.4.- ESTIMACIÓN DE PARÁMETROS El programa Ax+B es capaz de realizar el ajuste funciones de distribución de probabilidad a series de datos, esto mediante la aplicación de diversos procesos para estimar los parámetros de cada una como se menciona a continuación: ο Distribución Normal o Ajuste por Método de Momentos/Máxima Verosimilitud o Ajuste por el Método de Momentos L ο Distribución Log-Normal de dos parámetros o Ajuste por el Método de Momentos o Ajuste por el Método de Máxima Verosimilitud o Ajuste por el Método de Stigth ο Distribución Log-Normal de tres parámetros P á g i n a 9 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 ο ο ο ο ο ο ο ο ο o Ajuste por el Método de Momentos o Ajuste por el Método de Máxima Verosimilitud Distribución Exponencial de un parámetro o Ajuste por el Método de Momentos Distribución Exponencial de dos parámetro o Ajuste por el Método de Momentos o Ajuste por el Método de Máxima Verosimilitud Distribución Gamma de dos parámetros o Ajuste por el Método de Momentos o Ajuste por el Método de Máxima Verosimilitud o Ajuste por el Método de Momentos L Distribución Gamma de tres parámetros o Ajuste por el Método de Momentos o Ajuste por el método de Máxima Verosimilitud Distribución Gumbel o Ajuste por el Método de Momentos o Ajuste por el Método de Máxima Verosimilitud o Ajuste por el Método de Momentos L o Ajuste por el Método de Máxima Entropía o Ajuste por el Método de Mínimos Cuadrados Distribución Frechet(Log-Gumbel) o Ajuste por el Método de Momentos Distribución Gumbel de dos poblaciones (Doble Gumbel) o Ajuste por el Método de Momentos o Ajuste por el Método de Máxima Verosimilitud o Ajuste por el Método de Momentos L o Ajuste por el Método de Máxima Entropía o Ajuste por el Método de Mínimos Cuadrados Distribución Log-Normal de dos poblaciones o Ajuste por el Método de Momentos o Ajuste por el Método de Máxima Verosimilitud o Ajuste por el Método de Stigth Distribución Frechet_Doble o Ajuste por el Método de Momentos P á g i n a 10 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 2.4.1.- METODO DE MOMENTOS Este es un procedimiento muy sencillo para encontrar los estimadores de los parámetros poblacionales. Consiste básicamente en plantear un sistema de ecuaciones, cuyo tamaño depende del número de parámetros a estimar. Esto se hace al igualar los momentos poblacionales con los muéstrales. Los momentos muéstrales son parámetros estadísticos comunes: ο Media: 1 π₯Μ = βππ=1 π₯π π 1 2 ο Varianza Sesgada: π π ππ π = βππ=1(π₯π β π₯Μ )2 π ο Varianza no Sesgada: 2 π πππ ππ π = 1 πβ1 βππ=1(π₯π β π₯Μ )2 ο Coeficiente de asimetría sesgado: ππ ππ π = ο Coeficiente de asimetría no sesgado: ο Coeficiente de Curtosis sesgado: ο Coeficiente de Curtosis no sesgado ο Desviación Estándar: ππππ ππ π = ππ ππ π = 1 π β (π₯ βπ₯Μ )3 π π=1 π 2 (ππ ππ π )3/2 π2 (πβ1)(πβ2) β ππ ππ π 1 π β (π₯ βπ₯Μ )4 π π=1 π 2 (ππ ππ π )2 ππ ππ π = π3 (πβ1)(πβ2)(πβ3) β ππ ππ π π = βπ 2 ο Coeficiente de variación: πΆπ = π π₯Μ P á g i n a 11 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 2.4.2.- METODO DE MAXIMA VEROSIMILITUD Seaβ‘f(x, π1 , π2 , β¦ . ππ ) una función de densidad de probabilidad de x con parámetros ππ π = 1 β¦ π. Si existe una muestra aleatoria π₯1 , π₯2 β¦ . π₯π de esta función de densidad, entonces, su función de densidad conjunta es f(π₯1 , π₯2 β¦ . π₯π : π1 , π2 β¦ β¦ ππ ) . Debido a que la muestra es aleatoria, la función de densidad conjunta se puede escribir como en la ecuación: π π(π₯1 , π₯2 β¦ . π₯π : π1 , π2 β¦ β¦ ππ ) = β π(π₯1 : π1 , π2 β¦ β¦ ππ ) π=1 La probabilidad de obtener la muestra aleatoria π₯1 , π₯2 β¦ . π₯π a partir de la población de x es proporcional al producto de sus densidades de probabilidad individual. Esta función conjunta es llamada la función de verosimilitud L π πΏ = β π(π₯1 : π1 , π2 β¦ β¦ ππ ) π=1 Donde los parametros ππ , π = 1,2 β¦ β¦ β¦ π Debido a que con ln(L) se alcanza también su máximo para valores específicos ππ π = 1 β¦ π como lo hace L, entonces, la función de verosimilitud se puede expresar como: π lnβ‘(πΏ) = ππ β π(π₯1 : π1 , π2 β¦ β¦ ππ ) π=1 2.4.3.- METODO DE MOMENTOS L Los momentos lineales surgen de combinaciones lineales de los momentos ponderados probabilísticamente y constituyen un sistema alternativo al método tradicional de los momentos convencionales para describir las funciones de distribución. Los momentos ponderados probabilísticamente, de la variable aleatoria X con una función de distribución F(X), quedan definidos por las siguientes expresiones: π1 = π½0 π2 = 2π½1 β 2π½0 π3 = 6π½2 β 6π½1 + π½0 π4 = 20π½3 β 30π½2 + 12π½1 β π½0 π 1 π½0 = β π₯π π πβ1 π½1 = π=1 1 (π β π) β π₯ π (π β 1) π π=1 P á g i n a 12 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 πβ1 1 (π β π)(π β π β 1) π½2 = β π₯ π (π β 1)(π β 2) π πβ3 π½3 = π=1 1 (π β π)(π β π β 1)(π β π β 2) β π₯π π (π β 1)(π β 2)(π β 3) π=1 Donde: π½π β‘= Momentos de Probabilidad Ponderada de Orden n πβ‘= Numero de datos de la muestra π₯π β‘= Elemento i-enesimo de la serie 2.4.4.- MAXIMA ENTROPIA Tomando en consideración que la entropía puede ser vista como una medida de la incertidumbre para una función de densidad (Shannon ,1948), se tiene que π = β‘ β β« π(π¦) ln(π(π¦)) β π¦ En 1957 Jaynes propuso que el mejor ajuste es una distribución de probabilidad que minimice la información previa necesaria, que será aquella que maximice la entropía de Shannon. Entonces, la función de densidad de máxima entropía f(y) definida en el intervalo [a; b] es aquella que cumpla: π πππ₯ = β β« π(π¦) ln(π(π¦)) β π¦ π π β« π(π¦)β π¦ = 1 π π β« π¦ π π(π¦)β π¦ = ππ π Donde: ππ β‘ Son los M primeros momentos con respecto al origen π(π¦) Función de ajuste P á g i n a 13 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 2.5.- PRUEBAS DE BONDAD DE AJUSTE Para encontrar la función de distribución que cuenta con el mejor ajuste se han incluido en el software las siguientes pruebas de bondad: Error Cuadrático Kolmogorov- smirnov Chi Cuadrado Anderson Darling Nash-Sutclffe ο ο ο ο ο Cabe mencionar que el programa efectúa las pruebas de Kolmogorov-Smirnov y ChiCuadrado en base a la probabilidad acumulada teórica y real del dato observado. 2.5.1.-ERROR CUÁDRATICO Esta técnica permite seleccionar el mejor ajuste de los distintos modelos, fue propuesto por Kite en 1988. Este estadístico se obtiene con la ecuación: 1β 2 βππ=1(πΜππ β πππ )2 πΈπΈπ΄ = [ ] ππ β ππ Donde: π π ππ Son los eventos ππ ordenados de mayor a menor con un periodo de retorno asignado: π= ππ β1 π y una probabilidad de excedencia π = 1 β 1 π ππ : Longitud del registro analizado π: Numero de orden π πΜπ : Eventos estimados con alguna función de probabilidad para cada periodo de retorno ππ : Numero de parámetros de la distribución ajustada La distribución de mejor ajuste será aquella que proporcione el mínimo valor del estadístico E.E.A. Si una o más distribuciones tienen valores similares del E. E.A, entonces se deberá optar por aquella distribución que tenga el menor número de parámetros. P á g i n a 14 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 2.5.2.- KOLMOGOROV-SMIRNOF La prueba Kolmogorov Smirnov considera la desviación de la función de distribución de probabilidades de la muestra P(x) de la función de probabilidades teórica mediante la ecuación: Dn ο½ max( P( x) ο Po( x)) 2.5.3.- CHI CUADRADO Esta prueba mide las diferencias entre las frecuencias observadas (fo) y las frecuencias calculadas (fc) por medio de una distribución teórica, para ello emplea el estadístico Ο² k ο£2 ο½ ο₯ i ο½1 ( fo ο fc )2 fc Si el estadístico Ο²=0 significa que las distribuciones teórica y empírica ajustan exactamente. 2.5.4.- ANDERSON DARLING La prueba de Anderson-Darling es una prueba no paramétrica sobre si los datos de una muestra provienen de una distribución específica. La fórmula para el estadístico A determina si los datos vienen de una distribución con función acumulativa. De manera tal que π΄2 = βπ β π π π=β π=1 2π β 1 [lnβ‘(πΉ(π₯π ) + lnβ‘(1 β πΉ(π₯π+1βπ )) π 2.5.5.- NASH-SUTOLFFE A pesar de que en México el criterio de Nash-Sutcliffe es poco conocido es uno de los más usados en Hidrología fuera del país. Se define como: βπ1(ln(πππππ’ππβ π) β ln(πππ πππ£πβ π))2 πΈ =1β π β (ln(πππ πππ£πβ π) β Μ Μ Μ Μ Μ Μ Μ Μ Μ Μ Μ Μ Μ Μ Μ Μ Μ Μ Μ ln(πππ πππ£πβ π))2 1 P á g i n a 15 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 3.- UTILIZACIÓN DEL PROGRAMA 3.1.- REQUISITOS MINIMOS Para la utilización del programa la computadora del usuario final debe cumplir con ciertos requisitos mínimos de software y hardware: ο ο ο ο ο ο Sistema Operativo Windows 7 o superior Memoria Ram de 1GB Espacio libre en disco duro 2 Mb Velocidad mínima del procesador 3.2 GHz Resolución mínima de pantalla 1366x768 Net Framework 4.0 3.2.- INTERFAZ La interfaz del programa cuenta con: 7 1 6 2 3 5 4 Interfaz del Programa ο ο ο ο ο ο ο 1) 2) 3) 4) 5) 6) 7) Barra de Menú Tabla de Resultados Barra de Herramientas Grafico Tabla Resumen de la ejecución del Método Tabla de Parámetros de la distribución seleccionada Estadísticos de la Serie P á g i n a 16 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 3.2.- BARRA DE MENÚ La barra de Menú del programa cuenta con las siguientes opciones: ο ο ο ο Barra de Menú Menú Archivo: Desde este menú se tiene la posibilidad de abrir un nuevo archivo de entrada de datos y finalizar la ejecución del programa. Opciones: En este apartado se pueden configurar las opciones del programa y se puede optimizar de forma automática el parámetro p de las funciones de dos poblaciones en base al error cuadrático. Créditos: Desde este apartado puede visualizar los créditos de autoría del programa y la forma de ponerse en contacto con el/los autores. Ayuda: Desde este Menú puede acceder a la ayuda en línea, manuales y foros; de forma adicional es posible visualizar el contrato de licencia de este software. 3.3.- BARRA DE HERRAMIENTAS El programa cuenta con una barra de herramientas la cual nos permite manipular las opciones graficas del programa, ordenar los resultados en base a alguna de la pruebas de bondad de ajuste. Barra de Herramientas del Programa Ax+B Opciones de la barra de herramientas: ο Recalcular: Permite recalcular el ajuste en base a los cambios hechos en la configuración del programa (decimales, parámetro p, funciones de distribución objetivo). ο Ordenar: Permite ordenar los ajustes a partir de las pruebas de bondad de ajuste efectuadas. ο Extrapolar Datos: Obtiene el resultado del ajuste para cada uno de los periodos de retorno incluidos en la configuración. ο : Muestra gráficamente el ajuste de la función seleccionada en la tabla resumen, al cambiar de fila seleccionada, si esta opción esta deshabilitada se mostrara el en el grafico el ajuste de todas las funciones aplicadas. ο : Al seleccionar esta opción el programa mostrara al usuario en el eje βxβ del grafico la Probabilidad Acumulada, si esta opción esta deshabilitada se mostrara el periodo de retorno. : El grafico se mostrara en escala logarítmica si el usuario selecciona esta opción. De no ser así la escala será aritmética. ο P (%): Parámetro p de las funciones de distribución de dos poblaciones en porcentaje. ο P á g i n a 17 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 3.4.- CONFIGURACIÓN DEL PROGRAMA Para configurar el programa debe dirigirse a Opciones/Configuración en la barra de menú. Se desplegara la siguiente ventana desde la cual podrá configurar el programa: Ventana de Configuración En la primera pestaña se puede agregar o remover los periodos de retorno a extrapolar y modificar la cantidad de decimales a mostrar. Para remover un periodo de retorno de la lista basta con seleccionar la fila de dicho periodo y presionar la tecla Supr. Si se desea agregar un periodo adicional se tiene que dar click en el botón el cual desplegara un espacio de texto en el cual el usuario tecleara el periodo de retorno y finalizara con la tecla ENTER. En la segunda pestaña se muestran las distribuciones que el usuario puede aplicar, en esta misma es posible desactivar las que el usuario desee. P á g i n a 18 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 Configuración del Programa, Segunda Pestaña 3.5.-ENTRADA DE DATOS Para que el programa lea una nueva serie de datos el usuario deberá ir a Archivo/Abrir Archivo o utilizar la combinación de teclas Ctrl+O, donde se desplegara la siguiente ventana: Ventana para leer archivo de entrada P á g i n a 19 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 El archivo de entrada tendrá el formato β.txtβ que se refiere a un archivo de texto y deberá contener la siguiente estructura: Archivo de entrada de datos .txt Una vez que el archivo sea leído por el programa se ejecutara el ajuste de las funciones de distribución, este proceso podrá tomar algún tiempo dependiendo de la cantidad de datos que contenga la serie y la cantidad de funciones ajustadas. 3.6.- OPTIMIZACIÓN DE FUNCIONES DE DISTRIBUCIÓN DE DOS POBLACIONES Una vez que el usuario el programa ha realizado el primer ajuste a la serie el usuario tiene la posibilidad de optimizar el parámetro p en las funciones de dos poblaciones. Para realizar el ajuste el programa cuenta con Opciones/Optimizar Funciones de Dos Poblaciones. El programa desplegara el siguiente formulario: Ventana de Optimización de Funciones de dos Poblaciones En este el usuario colocara los rangos (en %) de iteración de la variable, cabe mencionar que el número mínimo de datos de una serie que toma en consideración el programa es 3. P á g i n a 20 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 3.7.- GRÁFICO DE FUNCIONES DE DISTRIBUCIÓN Posterior al ajuste el programa muestra al usuario un gráfico con las funciones de distribución y la serie original como el siguiente: Gráfico de Funciones de Probabilidad General El usuario tiene la posibilidad de copiar este grafico al portapapeles desde el botón . De forma adicional el usuario puede modificar la gráfica des de la barra de tareas, haciendo que esta se muestre en escala logarítmica, que grafique probabilidad acumulada o periodo de retorno y tiene la posibilidad de elegir entre mostrar el ajuste seleccionado o todos. P á g i n a 21 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 Gráfico de Funciones de Probabilidad Individual 3.8.- EXTRAPOLAR DATOS Una vez que el usuario ha elegido un ajuste, el proceso para extrapolar datos es simple. Basta con seleccionar el renglón de la distribución elegida y dar clic en el botón βExtrapolarβ de la barra de herramientas, acto seguido aparecerá en la pantalla la siguiente ventana: Ventana de Extrapolación de Datos P á g i n a 22 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 3.9 GENERACION DE REPORTES El software ofrece la posibilidad de generar reportes en formato βExcelβ (.xls) para su posterior impresión, para ello es necesario dar click en la opción Archivo\Generar Reporte Menú de Generación de Reportes Se elegirá la ruta donde se almacenara el archivo y el programa generara el reporte de forma automática como se muestra a continuación: PROGRAMA DE FUNCIONES DE DISTRIBUCION AX+B RESUMEN DE RESULTADOS Funcion Gumbel_Doble_MC Gumbel_Doble_ML LogNormal_2p_Stigth_Doble LogNormal_2p_Doble Frechet Gumbel_Doble_Momentos Frechet_Doble LogNormal_2p_MV_Doble Exponencial_2p Gumbel_Doble_MV Gamma_3p LogNormal_3p LogNormal_3p_MV Gumbel Gumbel_ML LogNormal_2p_Stigth LogNormal_2p Gumbel_ME Gamma_2p Gumbel_MV LogNormal_2p_MV Gamma_2p_ML Gamma_2p_MV Exponencial_2p_MV Normal Normal_ML Gumbel_MC Gamma_3p_MV Exponencial_1p Error_Cuad Orden 2.3894 2.4877 2.4955 2.5276 2.6834 2.7374 2.7637 2.7896 3.1196 3.1336 3.1472 3.3405 3.458 3.4874 3.624 3.6879 3.7012 4.0147 4.2299 4.3226 4.349 4.4878 4.7453 4.9664 5.5671 5.672 7.0579 9.3654 26.634 Kolmogorov 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 0.1092 0.0826 0.0928 0.0917 0.0855 0.0785 0.1126 0.0925 0.4814 0.0806 0.1063 0.1019 0.0996 0.114 0.1171 0.1215 0.1197 0.1248 0.1334 0.1181 0.1327 0.1421 0.1491 0.2431 0.1715 0.179 0.2088 0.3299 0.4133 Orden 2 11 3 7 5 4 1 12 6 29 2 10 9 8 13 14 17 16 18 20 15 19 21 22 26 23 24 25 27 28 Anderson Darling 1.7455 0.2368 0.7709 0.7101 0.2896 0.4287 0.1124 0.7955 1.2152 0.6616 7.4248 0.4419 0.4539 0.7859 0.9103 1.0508 0.9632 1.2423 1.457 0.9057 1.4055 1.7502 1.9985 2.9606 2.931 3.1266 6.5023 3.5593 3.0104 Orden 3 20 2 9 8 3 4 1 11 16 7 29 5 6 10 13 15 14 17 19 12 18 21 22 24 23 26 28 27 25 Ji Cuadrado 0.0896 0.1083 0.1205 0.1222 0.1433 0.1158 0.163 0.1596 0.7171 0.1819 0.1893 0.2002 0.1994 0.2398 0.2567 0.2851 0.2842 0.3202 0.4269 0.3698 0.4395 0.5597 0.7305 0.8138 1.6638 3.0965 2.2051 2.1219 3.5891 Orden 4 1 2 4 5 6 3 8 7 22 9 10 12 11 13 14 16 15 17 19 18 20 21 23 24 25 28 27 26 29 NashSutcliffe 0.9575 0.9736 0.9775 0.9775 0.9691 0.9692 0.9627 0.9755 0.948 0.9694 0.9561 0.9688 0.9671 0.964 0.9638 0.9573 0.9571 0.9587 0.9348 0.9543 0.9525 0.9405 0.9399 0.8377 0.8518 0.8828 0.8208 0.572 -11.4561 Orden 6 14 4 1 2 7 6 12 3 20 5 17 8 9 10 11 15 16 13 23 18 19 21 22 26 25 24 27 28 29 Generación ce Reportes en Excel P á g i n a 23 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 PROGRAMA DE FUNCIONES DE DISTRIBUCION AX+B Exponencial_1p PARAMETROS DE CALCULO Alfa 0.02068045 Tr 34 17 11.3333 8.5 6.8 5.6667 4.8571 4.25 3.7778 3.4 3.0909 2.8333 2.6154 2.4286 2.2667 2.125 2 1.8889 1.7895 1.7 1.619 1.5455 1.4783 1.4167 1.36 1.3077 1.2593 1.2143 1.1724 1.1333 1.0968 1.0625 1.0303 Dato 96.46 83.23 72.69 68.32 62.07 61.08 61.07 58.69 53.32 51.6 49.14 48.38 47.5 46.18 44 43.83 43.74 43.54 42.85 42.03 41.98 41.62 41.5 41.05 40.14 38.63 38.32 37.73 33.42 33.12 31.03 30.89 26.56 Calculado 170.5166 136.9996 117.3934 103.4826 92.6925 83.8764 76.4224 69.9655 64.2702 59.1755 54.5668 50.3593 46.4889 42.9054 39.5693 36.4485 33.517 30.7531 28.1387 25.6584 23.2992 21.0497 18.9003 16.8423 14.8684 12.9719 11.1469 9.3884 7.6915 6.0522 4.4667 2.9315 1.4435 Error Prob. Real Cuadratico 5484.3843 0.0294 2891.1701 0.0588 1998.3938 0.0882 1236.4067 0.1176 937.7377 0.1471 519.6745 0.1765 235.6974 0.2059 127.138 0.2353 119.9062 0.2647 57.3878 0.2941 29.4498 0.3235 3.9178 0.3529 1.0223 0.3824 10.7229 0.4118 19.6314 0.4412 54.4862 0.4706 104.5092 0.5 163.5038 0.5294 216.4215 0.5588 268.0277 0.5882 348.972 0.6176 423.1356 0.6471 510.7473 0.6765 586.0119 0.7059 638.655 0.7353 658.3399 0.7647 738.3754 0.7941 803.2472 0.8235 661.9532 0.8529 732.6634 0.8824 705.609 0.9118 781.678 0.9412 630.8368 0.9706 Prob. Teorica 0.136 0.1788 0.2224 0.2434 0.277 0.2828 0.2828 0.2971 0.332 0.344 0.362 0.3677 0.3744 0.3848 0.4025 0.404 0.4047 0.4064 0.4122 0.4193 0.4197 0.4229 0.4239 0.4279 0.436 0.4498 0.4527 0.4583 0.501 0.5041 0.5264 0.5279 0.5774 Generación de Reportes en Excel P á g i n a 24 | 25 MANUAL DE USUARIO FINAL AX+B VERSION 1.0.0.0 4.- BIBLIOGRAFIA ο Abramowitz, M. y Stegun, I.(1972), "Handbook of mathematical functions", Dover Publications, N. Y. ο Campos Aranda (2007), βIntroducción a los métodos numéricos: Software en Basic y aplicaciones en Hidrología Superficialβ ο Campos Aranda (2007), βEstimación y Aprovechamiento del Escurrimientoβ, San Luís de Potosí, México: Editorial Campos. ο Escalante Sandoval (2005), βTécnicas Estadísticas en Hidrologíaβ, UNAM, D.F. México. P á g i n a 25 | 25
© Copyright 2024