La Ciencia está obligada a informar al público lo que ocurre en el entorno. Pero no basta con opiniones ni comentarios: eso no es ciencia. Es charlatanería. Los únicos informes VÁLIDOS provienen de datos de mediciones con instrumentos. ¿El propósito? Predecir lo que ocurrirá y pronosticar CUÁNDO ocurrirá.
Las ciencias químicas, ciencias ambientales y ciencia de datos, en conjunto, producen este tipo de información para tomar DECISIONES acertadas, en tiempo y forma. Los modelos ARIMA son el mecanismo legítimo para pronosticar. Este es un informe con validez estadística para el pronóstico de calidad de agua residual tratada a 10 días. ¿Preguntas? Q.I. Víctor Soto Del Toro, victor[a]apestoff.com, Tijuana, BC, México
Shapiro-Wilk normality test
data: df.residuales.bod.out$Residuales
W = 0.98874, p-value = 0.3166
Box-Ljung test
data: df.residuales.bod.out$Residuales
X-squared = 0.050704, df = 1, p-value = 0.8218
ANTECEDENTES
En la ciencia ambiental, el control de la contaminación del agua se ejerce de diversas maneras.
Una de ellas es la ESTIMACIÓN (tan exacta y precisa como sea posible), de la CANTIDAD de contaminantes de interés presentes.
Pero no sólo eso: es aún más crítico pronosticar CUÁNDO ocurrirán los cambios de dichas cantidades.
Lo que a final de cuentas más nos interesa, es conocer de qué tamaño es el problema y el momento en que ocurrirá.
Así, el desafío es doble:
1. ¿Cuánta es la magnitud de la contaminación del agua?
2. ¿Cuándo se presentarán diferentes niveles de contaminación?
¿Y cómo se responde eso?
Con dos métodos de la Ciencia de Datos: predicción y pronóstico.
¿Por qué es esto relevante? ¿Para qué sirve?
Para el desarrollo y ejecución de Estrategias de Control:
1. Estimamos los niveles de contaminación para verificar el cumplimiento de la normatividad y a la vez, cuantificar el volumen de recursos necesarios para mitigarla.
2. Pronosticamos incrementos repentinos de la contaminación; con ello nos prevenimos con bastante tiempo de anticipación (DÍAS o inclusive SEMANAS) porque al monitorear los datos, estaremos siempre listos y alertas en todo momento para entrar cómodamente en acción, cuando los datos nos lo indiquen. Así se evitan sorpresas.
Es así como los métodos de Regresión Múltiple y Análisis de Series Temporales de la Ciencia de Datos nos permiten vigilar procesos que cambian con el transcurso del tiempo, para así evitar o al menos mitigar, costosos eventos de pérdida del control.
Es a partir de esta sofisticada tecnología de análisis de datos, como se toman las decisiones más inteligentes, responsables, serias, válidas, informadas y acertadas.
HIPÓTESIS
1. La demanda bioquímica de oxígeno en la descarga de la planta tratadora de aguas residuales, no exhibe tendecia.
2. La variación de la demanda bioquímica de oxígeno en la descarga, es constante.
DISCUSIÓN DE RESULTADOS
1. En el gráfico de monitoreo diario de DBO en la sección de Pronósticos podemos observar que sí existe una tendencia a la baja, desde el inicio hacia el final de la serie de datos. Una de las posibles causas es la temperatura del clima: a mayores temperaturas, menores niveles de DBO. Para propósitos de pronósticos futuros, las tendencias en los datos provocan interferencias. Por ello, durante la limpieza de los datos se removió todo componente de tendencia.
2. En la sección de Validación, estudiamos los RESIDUALES, que son los valores que resultan de la diferencia entre las mediciones realizadas de DBO en el laboratorio y el ajuste de la regresión del modelo ARIMA. El análisis estadístico muestra con detalle que los residuales provienen solamente de una única distribución y además, la conducta de la dispersión sugiere una curva Normal. Esto se confirma con la prueba de Shapiro-Wilk.
3. En la Función de Auto-Correlación de la misma sección de Validación, vemos que los residuales están prácticamente libres de autocorrelación y sólo queda “ruido blanco”. Esto significa que toda la información de fluctuación de DBO en el tiempo contenida en los datos, ha sido extraída. La prueba de Ljung-Box lo confirma.
HALLAZGOS Y CONCLUSIONES
1. Los datos de mediciones en serie de tiempo para demanda bioquímica de oxígeno, en la descarga de una planta tratadora de aguas residuales, exhiben 1) tendencia a la baja y 2) variación constante.
2. De lo anterior, fue posible construir un modelo autorregresivo integrado de promedio móvil, mejor conocido como ARIMA y el orden resultó de (13, 1, 13).
3. CONCLUIMOS que el modelo ARIMA obtenido es confiable para hacer pronósticos del nivel de contaminación del agua por DBO para esta planta. El límite máximo del pronóstico es hasta 10 días después de la última medición.
En la sección de Resultados se ofrecen los 10 pronósticos de DBO con sus respectivos intervalos de confianza al 80% y 95% y en Pronósticos se representan gráficamente en color púrpura tenue y obscuro, respectivamente.
Pronósticos de DBO en la descarga de aguas residuales.
| Día | Lim. Inf. 95% | Lim. Inf. 80% | Promedios | Lim. Sup. 80% | Lim. Sup. 95% |
|---|---|---|---|---|---|
| 141 | 8.48 | 11.54 | 17.31 | 23.07 | 26.13 |
| 142 | 7.48 | 10.88 | 17.28 | 23.69 | 27.09 |
| 143 | 8.47 | 12.06 | 18.83 | 25.61 | 29.20 |
| 144 | 7.99 | 11.63 | 18.51 | 25.39 | 29.04 |
| 145 | 8.19 | 11.97 | 19.12 | 26.26 | 30.04 |
| 146 | 8.17 | 12.09 | 19.49 | 26.89 | 30.80 |
| 147 | 9.82 | 13.84 | 21.42 | 29.00 | 33.02 |
| 148 | 5.78 | 10.38 | 19.08 | 27.77 | 32.38 |
| 149 | 4.33 | 9.14 | 18.23 | 27.31 | 32.12 |
| 150 | 6.10 | 11.05 | 20.41 | 29.78 | 34.73 |
Interpretación:
1. La primera columna, Día, indica el calendario del pronóstico, desde el día 141 al 150.
2. La columna, Promedios, indica las concentraciones pronosticadas de Demanda Bioquímica de Oxígeno (en mg/L) esperadas para cada uno de los 10 días.
3. Las demás columnas indican los límites estimados, para niveles de confianza del 80% y 95%.
(Estos valores se representan gráficamente con colores púrpura claro y obscuro en la sección de Pronósticos.)
Todos los cálculos estadísticos, elementos gráficos y el cuerpo de este documento se escribieron completamente en código del lenguaje R.
Realizado por: Q. I. Víctor Soto Del Toro.
victor[a]apestoff.com
Fuente:
Statistics For Experimenters, An Introduction to Statistics, Data Analysis and Model Building
Box, Hunter, Hunter, 1a. Ed 1978, pp 626-628
Este es el texto clásico que los profesionales de las ciencias aplicadas (física, química, biología, materiales, medio ambiente, etc), consultan para especializarse en Ciencia de Datos.
| Día | DBO Afluente, (mg/L) | DBO Efluente (mg/L) | Temperatura (Celsius) | SST (mg/L) |
|---|---|---|---|---|
| 1 | 202 | 37 | 53 | 1100 |
| 2 | 202 | 37 | 53 | 1100 |
| 3 | 141 | 16 | 51 | 1100 |
| 4 | 173 | 27 | 48 | 1100 |
| 5 | 233 | 29 | 44 | 1120 |
| 6 | 214 | 23 | 48 | 1090 |
| 7 | 208 | 26 | 49 | 1130 |
| 8 | 227 | 26 | 48 | 1180 |
| 9 | 227 | 26 | 53 | 1540 |
| 10 | 201 | 12 | 51 | 1540 |
| 11 | 182 | 27 | 50 | 1240 |
| 12 | 199 | 28 | 50 | 1200 |
| 13 | 182 | 22 | 48 | 1300 |
| 14 | 222 | 24 | 49 | 1210 |
| 15 | 206 | 25 | 48 | 1250 |
| 16 | 206 | 25 | 52 | 1320 |
| 17 | 147 | 10 | 50 | 1320 |
| 18 | 173 | 20 | 44 | 1160 |
| 19 | 204 | 18 | 44 | 1330 |
| 20 | 194 | 14 | 46 | 1340 |
| 21 | 287 | 28 | 50 | 1200 |
| 22 | 185 | 26 | 53 | 1880 |
| 23 | 120 | 12 | 53 | 1330 |
| 24 | 120 | 12 | 45 | 1330 |
| 25 | 189 | 22 | 48 | 1280 |
| 26 | 193 | 25 | 47 | 1250 |
| 27 | 172 | 18 | 43 | 1300 |
| 28 | 197 | 27 | 40 | 1260 |
| 29 | 193 | 24 | 44 | 1290 |
| 30 | 193 | 24 | 44 | 1380 |
| 31 | 287 | 16 | 43 | 1380 |
| 32 | 185 | 36 | 46 | 1250 |
| 33 | 120 | 30 | 44 | 1110 |
| 34 | 120 | 17 | 45 | 1130 |
| 35 | 189 | 31 | 52 | 1160 |
| 36 | 193 | 23 | 52 | 1380 |
| 37 | 172 | 23 | 47 | 1410 |
| 38 | 197 | 8 | 48 | 1410 |
| 39 | 193 | 26 | 45 | 1380 |
| 40 | 193 | 20 | 45 | 1340 |
| 41 | 192 | 19 | 46 | 1380 |
| 42 | 239 | 24 | 50 | 1320 |
| 43 | 177 | 18 | 47 | 1330 |
| 44 | 177 | 18 | 50 | 1350 |
| 45 | 166 | 6 | 48 | 1350 |
| 46 | 190 | 16 | 50 | 1310 |
| 47 | 171 | 19 | 48 | 1360 |
| 48 | 201 | 13 | 50 | 1310 |
| 49 | 188 | 19 | 53 | 1300 |
| 50 | 174 | 19 | 60 | 1230 |
| 51 | 174 | 19 | 50 | 1330 |
| 52 | 136 | 7 | 50 | 1330 |
| 53 | 159 | 14 | 42 | 1470 |
| 54 | 216 | 18 | 44 | 1400 |
| 55 | 204 | 13 | 51 | 1440 |
| 56 | 250 | 21 | 47 | 1310 |
| 57 | 233 | 22 | 54 | 1370 |
| 58 | 233 | 22 | 53 | 1330 |
| 59 | 199 | 8 | 50 | 1330 |
| 60 | 209 | 16 | 46 | 1360 |
| 61 | 185 | 15 | 46 | 1310 |
| 62 | 203 | 17 | 47 | 1170 |
| 63 | 225 | 28 | 50 | 1360 |
| 64 | 225 | 23 | 52 | 1190 |
| 65 | 225 | 23 | 51 | 1290 |
| 66 | 119 | 10 | 50 | 1290 |
| 67 | 174 | 16 | 46 | 1220 |
| 68 | 190 | 18 | 52 | 1280 |
| 69 | 238 | 22 | 50 | 1290 |
| 70 | 225 | 24 | 50 | 1240 |
| 71 | 257 | 24 | 54 | 1250 |
| 72 | 257 | 24 | 50 | 1030 |
| 73 | 121 | 10 | 51 | 1030 |
| 74 | 163 | 14 | 50 | 1170 |
| 75 | 144 | 15 | 50 | 1300 |
| 76 | 172 | 10 | 49 | 1250 |
| 77 | 226 | 18 | 48 | 1190 |
| 78 | 184 | 11 | 48 | 1260 |
| 79 | 184 | 11 | 51 | 1280 |
| 80 | 184 | 8 | 49 | 1280 |
| 81 | 177 | 17 | 48 | 1170 |
| 82 | 160 | 18 | 47 | 1190 |
| 83 | 183 | 17 | 50 | 1230 |
| 84 | 144 | 25 | 53 | 1070 |
| 85 | 159 | 18 | 53 | 1090 |
| 86 | 159 | 18 | 50 | 1130 |
| 87 | 118 | 8 | 52 | 1130 |
| 88 | 159 | 19 | 47 | 1110 |
| 89 | 147 | 14 | 52 | 1080 |
| 90 | 154 | 16 | 54 | 1110 |
| 91 | 179 | 14 | 56 | 1020 |
| 92 | 222 | 15 | 47 | 1040 |
| 93 | 222 | 15 | 48 | 1180 |
| 94 | 162 | 10 | 50 | 1180 |
| 95 | 196 | 22 | 49 | 1150 |
| 96 | 212 | 25 | 53 | 1220 |
| 97 | 205 | 18 | 56 | 1210 |
| 98 | 192 | 15 | 55 | 1210 |
| 99 | 170 | 9 | 53 | 1310 |
| 100 | 170 | 9 | 52 | 1330 |
| 101 | 106 | 9 | 54 | 1330 |
| 102 | 140 | 27 | 53 | 1020 |
| 103 | 143 | 25 | 50 | 1310 |
| 104 | 161 | 19 | 52 | 1240 |
| 105 | 181 | 11 | 54 | 1110 |
| 106 | 156 | 15 | 54 | 1120 |
| 107 | 156 | 15 | 57 | 1180 |
| 108 | 139 | 6 | 56 | 1180 |
| 109 | 167 | 17 | 54 | 1160 |
| 110 | 155 | 13 | 56 | 1220 |
| 111 | 172 | 11 | 54 | 1210 |
| 112 | 212 | 15 | 55 | 1120 |
| 113 | 253 | 13 | 58 | 1060 |
| 114 | 253 | 13 | 53 | 1140 |
| 115 | 134 | 6 | 55 | 1140 |
| 116 | 177 | 14 | 53 | 1100 |
| 117 | 192 | 13 | 52 | 1180 |
| 118 | 179 | 22 | 55 | 1250 |
| 119 | 193 | 13 | 56 | 1170 |
| 120 | 216 | 10 | 58 | 1120 |
| 121 | 216 | 10 | 56 | 1140 |
| 122 | 137 | 8 | 56 | 1140 |
| 123 | 185 | 18 | 55 | 1140 |
| 124 | 156 | 13 | 58 | 1210 |
| 125 | 160 | 8 | 58 | 1250 |
| 126 | 186 | 14 | 60 | 1120 |
| 127 | 200 | 14 | 59 | 1210 |
| 128 | 200 | 14 | 58 | 1310 |
| 129 | 166 | 8 | 58 | 1310 |
| 130 | 166 | 13 | 57 | 1200 |
| 131 | 158 | 15 | 61 | 1250 |
| 132 | 213 | 15 | 56 | 1270 |
| 133 | 205 | 14 | 58 | 1300 |
| 134 | 171 | 11 | 58 | 1300 |
| 135 | 171 | 11 | 62 | 1430 |
| 136 | 194 | 17 | 62 | 1430 |
| 137 | 191 | 16 | 59 | 1660 |
| 138 | 169 | 19 | 64 | 1610 |
| 139 | 173 | 19 | 64 | 2010 |
| 140 | 183 | 21 | 59 | 1480 |