Pronósticos

Column 1

Modelo Aplicado: ARIMA(13, 1, 13)

Column 2

Función de Auto-Correlación: Lag Observado a 7 Días

Interpretación del Análisis de Series Temporales: Protección Ambiental

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

Validación

Column 1

Pruebas de Diagnóstico: Validación Estadística del Modelo

Column 2

Prueba de Cuantiles-Cuantiles

Pruebas Analíticas: Shapiro y Box


    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

Conclusiones

Column 1

Column 2

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.

Column 3

Column 4

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.

Column 5

Resultados

Column 1

Column 2

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.)

Column 3

Column 4




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

Column 5

Datos

Column 1

Column 2

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

Column 3