{"id":892,"date":"2010-04-16T09:05:12","date_gmt":"2010-04-16T14:05:12","guid":{"rendered":"http:\/\/blog.espol.edu.ec\/duval\/?p=892"},"modified":"2010-04-16T15:59:27","modified_gmt":"2010-04-16T20:59:27","slug":"2010-04-16","status":"publish","type":"post","link":"https:\/\/blog.espol.edu.ec\/duval\/2010\/04\/16\/2010-04-16\/","title":{"rendered":"2010-04-16 Verificaci\u00f3n satisfactoria de algoritmos"},"content":{"rendered":"<p>Observando los \u00faltimos dos resultados en cuanto al % de violaciones, se lleg\u00f3 a la conclusi\u00f3n de que exist\u00eda un error recurrente. Este afectaba siempre a la primera Muestra y a cualquier otra que se iniciara con las mismas condiciones.<\/p>\n<p>Al parecer el error era generado al momento de leer los archivos, porque el procedimiento es el siguiente:<\/p>\n<ul>\n<li>Hay 3 m\u00f3dulos: Generador aleatorio de las muestras; Generador de Pol\u00edtica Optima usando la media de las muestras anteriores; y Buscador de Pol\u00edtica Robusta usando la \u00f3ptima y las muestras de caudales.<\/li>\n<li>El 3ro espera los resultados del 2do y el 2do espera los resultados del 1ro.<\/li>\n<li>En la primera iteraci\u00f3n el 2do comenzaba a leer datos que todav\u00eda no eran completamente fijados en memoria y lo mismo hac\u00eda el 3ro.<\/li>\n<\/ul>\n<p>Lo solucionamos utilizando la funci\u00f3n <a href=\"http:\/\/msdn.microsoft.com\/en-us\/library\/ms686298\">Sleep<\/a> mientras se esperaba para leer el archivo del m\u00f3dulo anterior.<\/p>\n<p>Realizamos una prueba con M=5 y N=100 y este fue el resultado.<br \/>\n<img decoding=\"async\" src=\"http:\/\/blog.espol.edu.ec\/duval\/files\/2010\/04\/03_resultados-preliminares-05-muestras.png\" alt=\"\" \/><\/p>\n<p>Como podemos notar, ya no tiene el error anterior, adem\u00e1s los resultados son como los esper\u00e1bamos (del 30 al 50% de violaciones).<\/p>\n<p>Con esto resuelto, procedimos a probar uno de las pol\u00edticas cuasi robustas obtenidas (la 1ra, con 42% de violaciones)en Matlab, y como deber\u00eda ser, dieron los mismos porcentajes de fallos en su muestra correspondiente.<\/p>\n<p>Luego la probamos con todas las 10.000 series y obtuvimos 44% de violaciones, lo cual es muy satisfactorio porque quiere decir que con esta pol\u00edtica satisfacemos el 56% de las 10.000 series de caudales. \u00a1Lo cual nos trae mucha tranquilidad y alegr\u00eda!<\/p>\n<p>Con esto en mente probamos la mejor de las 5 pol\u00edticas encontradas (35 % de violaciones en su muestra) y nos dio 42% de violaciones en las 10.000 series, lo que representa 58% de satisfacci\u00f3n. Otro dato significativo es que 23 fue el n\u00famero m\u00e1ximo de violaciones que provoc\u00f3 en las no pudo satisfacer y en general fueron menores de 10.<\/p>\n<p>Y para hacer una \u00faltima verificaci\u00f3n, vamos a aumentar el % de variaci\u00f3n del 30 al 50% y realizar una ejecuci\u00f3n con M=5 y N=100, esto tomara aproximadamente 1h 20m.<\/p>\n<p>Lo obtenido de esta ejecuci\u00f3n es:<br \/>\n<img decoding=\"async\" src=\"http:\/\/blog.espol.edu.ec\/duval\/files\/2010\/04\/03_resultados-05-muestras-50porc-variacion.png\" alt=\"\" \/><br \/>\nLa mejor tuvo un 19% de violaciones y al probarla en la 10.000 muestras obtuvo 31% de violaciones, es decir, 69% de satisfacci\u00f3n y con 17 violaciones como m\u00e1ximo. Puede mejorar!<\/p>\n<p>Lo siguiente en hacer es una variaci\u00f3n en el m\u00f3dulo de b\u00fasqueda del robusto, la base actual para la b\u00fasqueda es la Pol\u00edtica \u00d3ptima, lo que deseamos hacer ir almacenando la mejor Pol\u00edtica Robusta y si esta es mejor que la \u00d3ptima de la Muestra actual, usarla como base para la b\u00fasqueda.<\/p>\n<p>Bueno, los resultados obtenidos con este cambio no cambiaron casi nada, no porque no fuera efectivo si no que hay dos condiciones que considerar.<\/p>\n<ol>\n<li>Los \u00d3ptimos son todos muy parecidos y siempre (hasta el momento) son menores que la Robusta. ver las gr\u00e1ficas.<\/li>\n<li>Solo est\u00e1bamos guardando como Pol\u00edtica Robusta aquella que estuviera dentro del rango aceptable (10%), y hasta el momento no hemos obtenido ninguna.<\/li>\n<\/ol>\n<p>Por estas razones estoy proponiendo lo siguiente:<\/p>\n<ul>\n<li>Seguir generando un n\u00famero M de muestras, cada una de N series.<\/li>\n<li>Calcular una y solo una pol\u00edtica \u00f3ptima, usando la media de una de las muestras. (pr\u00e1cticamente todas las \u00f3ptimas son iguales en cuanto a su costo)<\/li>\n<li>Usar esa pol\u00edtica como base inicial para buscar una \u00f3ptima y para las siguientes muestras usar la mejor Robusta que vayamos encontrando.<\/li>\n<\/ul>\n<p>Probando lo \u00faltimo esperamos tener mejores resultados. En un principio vamos a considerar como mejor a aquel que tenga menor % de violaciones, en un siguiente experimento lo haremos tomando en cuenta el costo de operaci\u00f3n que generan la pol\u00edticas, luego definiremos como usar las dos al mismo tiempo.<\/p>\n<p>No vamos a colocar ning\u00fan gr\u00e1fico como los anteriores pero, describir\u00e9 el comportamiento.<\/p>\n<ul>\n<li>Con la \u00d3ptima de la 1ra Muestra la b\u00fasqueda comenz\u00f3 con 100% de violaciones y finaliz\u00f3 en 17%.<\/li>\n<li>En la 2da muestra, tomando como base la Robusta anterior, inici\u00f3 con 29% de violaciones, en la generaci\u00f3n 1.700 estaba en el 13% -parece no va a mejorar- y finaliz\u00f3 con 13%.<\/li>\n<li>En la 3ra muestra, tomando como base la Robusta anterior, inici\u00f3 con 18% de violaciones, en la generaci\u00f3n 500 estaba en el 12%, en la generaci\u00f3n 800 estaba en el 8% -lo m\u00e1s probable es que no cambie- y termin\u00f3 en 8%.<\/li>\n<li>En la 4ta muestra, tomando como base la Robusta anterior, inici\u00f3 con 4%  de violaciones -lo cual es excelente-, aunque puede significar simplemente que la misma pol\u00edtica satisface mejor a la 4ta que a la 3ra, es muy probable que no mejore pero con esto hemos superado las expectativas que ten\u00edamos. Termin\u00f3 con 4%.<\/li>\n<li>En la 5ta muestra, tomando como base la Robusta anterior, inici\u00f3 con 2%   de violaciones lo cual es wahu. Bueno hasta el lunes<\/li>\n<\/ul>\n<p>Todo esto fue trabajado con M = 5, N = 100,\u00a0 un 50% de variaci\u00f3n y con <strong>5.000 generaciones<\/strong> para la b\u00fasqueda del Robusto. Lo m\u00e1s probable es que no podamos presentar la Pol\u00edtica \u00d3ptima casi robusta hoy por el tiempo (terminar\u00e1 la ejecuci\u00f3n a las 16h30 y sola hay carros hasta las 16h00), pero el lunes estar\u00e1 a primera hora.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Observando los \u00faltimos dos resultados en cuanto al % de violaciones, se lleg\u00f3 a la conclusi\u00f3n de que exist\u00eda un error recurrente. Este afectaba siempre a la primera Muestra y a cualquier otra que se iniciara con las mismas condiciones. Al parecer el error era generado al momento de leer los archivos, porque el procedimiento [&hellip;]<\/p>\n","protected":false},"author":2965,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[17155],"tags":[],"class_list":["post-892","post","type-post","status-publish","format-standard","hentry","category-optimo-robusto"],"_links":{"self":[{"href":"https:\/\/blog.espol.edu.ec\/duval\/wp-json\/wp\/v2\/posts\/892","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/blog.espol.edu.ec\/duval\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/blog.espol.edu.ec\/duval\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/duval\/wp-json\/wp\/v2\/users\/2965"}],"replies":[{"embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/duval\/wp-json\/wp\/v2\/comments?post=892"}],"version-history":[{"count":37,"href":"https:\/\/blog.espol.edu.ec\/duval\/wp-json\/wp\/v2\/posts\/892\/revisions"}],"predecessor-version":[{"id":896,"href":"https:\/\/blog.espol.edu.ec\/duval\/wp-json\/wp\/v2\/posts\/892\/revisions\/896"}],"wp:attachment":[{"href":"https:\/\/blog.espol.edu.ec\/duval\/wp-json\/wp\/v2\/media?parent=892"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/duval\/wp-json\/wp\/v2\/categories?post=892"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/duval\/wp-json\/wp\/v2\/tags?post=892"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}