Generar script que integre el procedimiento y los análisis estadísticos necesarios para obtener el modelo de lógica difusa de encharcamientos e inundaciones

[iaga]

A partir de los resultados obtenidos en los scripts:

• Índice de vulnerabilidad. modelo_fuzzy_inunda.Rmd
• Análisis estadísticos. script_exploratorios.Rmd
• Cuantiles. descriptivo_fuzzy\categorias_fuzzy_inunda.Rmd

Generar un script en R que integre todos los resultados del modelo de lógica difusa para que puedan ser utilizados en la obtención de las reglas de transición entre categorías.

Insumos:

Los insumos se encuentran ubicados en la ruta: \SIG\desarrollo\sig_megadapt\procesamiento\fuzzy_inunda

Producto:

Funcion descriptiva fuzzy.Rmd \SIG\desarrollo\sig_megadapt\procesamiento\fuzzy_inunda\descriptivo_fuzzy

[iaga]

Características del producto final:

El producto final debe ser un script general que permite, a partir de cualquier clasificación (Weber-Fechner, Cuantiles), obtener los valores defusificados para las variables indicadas (por pares) y la estadística respectiva. Además, se generarán los data frames y las tablas correspondientes.

• @JosueMA generará la primera versión de integración del script
• @serranoycandela, @iaga y @aleestrada depurarán la versión final

Primera versión de integración del script

Responsable: @JosueMA

<Documentar en esta sección>

Fecha de entrega: 23/04/2019

Actividades a realizar para depurar la primera versión del script :

Responsable: @JosueMA

• [ ] Cumplir con el protocolo de nombramiento de archivos (quitar mayúsculas, quitar espacios). Es necesario cambiar el título del script.
• [x] La tabla generada para el análisis de la capacidad debe ser de color blanco. La tabla no debe incluir el gradiente de colores.
• [x] Generar las tablas de contingencia correctamente. Las tablas de contingencia son idénticas para encharcamientos/inundaciones y contienen frecuencias incorrectas.

Responsable: @serranoycandela, @iaga y @aleestrada

• [x] Las tablas de contingencia deben contener el nombre de la variable que se está mostrando. Es necesario especificar a qué se refiere "cat_1", "cat_2".
• [x] La variable capacidad debe ser la variable normalizada. El script está utilizando la variable capacidad sin normalizar, al cambiar a la variable capacidad_n la estadística no compila.
• [x] La capacidad promedio normalizada debe incluirse dentro de la tabla. El script actual está calculando únicamente la capacidad promedio sin normalizar.

Encharcamientos.

Responsables: @iaga, @serranoycandela y @aleestrada Archivo generado: modelo_general_encharcamientos.Rmd

• [x] Datos promedio (mínimo). Generar un data frame que contenga:

cat_1	cat_2	num_agebs	mean_freq_ench_n	mean_capacidad_n	mean_escorrentia_n

• [x] Datos normalizados. Generar una tabla que contenga por categoría:

• Número de AGEBS

• Capacidad mínima y máxima

• Promedio normalizado

• Diferencia entre la capacidad máxima y mínima

• [x] Datos normalizados. Generar una data frame que contenga:

• cat_1

• cat_2

• Número de AGEBS

• Capacidad mínima y máxima

• Promedio normalizado

• Diferencia entre la capacidad máxima y mínima

• Fuzzy value

Inundaciones.

Responsables: @iaga, @serranoycandela y @aleestrada Archivo generado: modelo_general_inundaciones.Rmd

• [x] Datos promedio (mínimo). Generar un data frame que contenga:

cat_1	cat_2	num_agebs	mean_freq_inund_n	mean_capacidad_n	mean_escorrentia_n

• [x] Datos normalizados. Generar una tabla que contenga por categoría:

• Número de AGEBS

• Capacidad mínima y máxima

• Promedio normalizado

• Diferencia entre la capacidad máxima y mínima

• [x] Datos normalizados. Generar una data frame que contenga:

• cat_1

• cat_2

• Número de AGEBS

• Capacidad mínima y máxima

• Promedio normalizado (capacidad)

• Diferencia entre la capacidad máxima y mínima

[JosueMA]

Los scripts y los reportes compilados, se encuentran en:

• Inundaciones: (LANCIS)\SIG\desarrollo\sig_megadapt\procesamiento\fuzzy_inunda\descriptivo_fuzzy\modelo_general_inundaciones_v2

• Encharcamientos: (LANCIS)\SIG\desarrollo\sig_megadapt\procesamiento\fuzzy_inunda\descriptivo_fuzzy\modelo_general_encharcamientos_v2

[iaga]

Tabla de contingencia para la capacidad (por cuantiles)

Observaciones:

• A mayor cantidad de encharcamientos, mayor capacidad.
• Al comparar con los resultados obtenidos en el modelo de panel, el coeficiente de la capacidad es un parámetro positivo o cero (dependiente del cluster).
• Por lo tanto, es necesario definir cuál será la acción del SAC en el modelo general @serranoycandela

Por ejemplo, para el cluster 2 con el modelo de panel:

Coefficients:	Estimate	Std. Error	t-value	Pr
(Intercept)	-21.7634447	1.6350629	-13.3105	< 2.2e-16 ***
log(preci_v)	1.0582101	0.1180459	8.9644	< 2.2e-16 ***
log(capacidad + 1e-05)	0.1143820	0.0259775	4.4031	1.114e-05 ***
log(escorren + 1e-05)	-0.1480207	0.0187083	-7.9120	3.831e-15 ***
log(rejillas + 1e-05)	-0.0382930	0.0239645	-1.5979	0.1101959
log(bombeo + 1e-05)	-0.1602258	0.0544871	-2.9406	0.0033069 **
log(centra + 1e-05)	0.0068618	0.0186141	0.3686	0.7124326
log(alcant + 1e-05)	0.0798467	0.0228800	3.4898	0.0004921 ***

[edithvilla]

@serranoycandela La coordinadora de MEGADAPT está solicitando que la apoyemos con la siguiente información:

Con qué puntos pueden apoyar?, Ale e Ileana van a apoyar con el punto 1 y 3

Mensaje: "Estoy escribiendo un articulo sobre los modelos mentales de MEGADAPT. Necesito unos estadísticas de nuestro base de datos, y espero que me pueden ayudar:

1. % de colonias / AGEBS que definimos como “wáter scarce” para propósito de Megadapt (y como definimos escasez de agua)

2. % de AGHEBS expuestas a echarcamientos (o alguna cifra que representa que tanto la ciudad sufre inundaciones)

3. Y luego, para Magdalena Contreras, Xochimilco y Izta, como su situación compara a la ciudad en general (% AGEB Izta comparado a % AGEB ciudad en general).

4. También me acuerda que hiciste un análisis de las estadísticas socio-económicas de las AGEBs expuestas. Me encantaría revisar estos datos de nuevo. Seria bueno tener esta análisis a la mano de todos modos.

Podemos reunirnos esta semana para este propósito? Estoy copiando Paola y Luis para no hacer lata con la organización de sus tiempos!

Muchas gracias"