# Contrastes corregidos: anual 2025 y periodos disponibles de 2026

Reanálisis del 6 de noviembre de 2026. Dirección científica: PhD. Rene Flores Clavo.

## Qué se compara

La referencia utiliza enero–diciembre de 2025 y FACTOR07 sin transformación. Los archivos analizados de 2026 cubren T1 y T2; se usa FACTOR de cada módulo y trimestre. No hay una estimación anual 2026 en este análisis.

El contraste principal compara cada indicador anual 2025 con el promedio aritmético de las dos estimaciones trimestrales de 2026. Cada trimestre aporta la mitad: es un promedio temporal de dos razones, no un factor semestral oficial ni el cociente de los totales apilados. La familia principal contiene 50 contrastes: diez indicadores por cuatro regiones y el total nacional. Como análisis secundario se comparan los dos trimestres por separado con el anual, en una familia de 100 contrastes. Se aplica Holm por familia; los IC95 % son puntuales, no simultáneos.

Se retiraron de la presentación las 200 estimaciones trimestrales 2025 construidas con peso anual. Las 100 comparaciones anteriores del “mismo trimestre entre años” se sustituyeron por anual 2025 frente a cada trimestre 2026. No se obtuvo un peso trimestral multiplicando FACTOR07 por cuatro.

## Estimación e incertidumbre

Las proporciones y medias de cada ola se calculan con su propio denominador válido y factor. Se conserva el universo de cada indicador y se construye el diseño antes de estimar los dominios. La diferencia principal es:

    D = (estimación T1 2026 + estimación T2 2026) / 2 − estimación anual 2025.

Su varianza incluye las covarianzas de las tres estimaciones. Se utilizan UPM compartidas cuando coinciden CONGLOME, UBIGEO, DOMINIO y ESTRATO. Esta correspondencia territorial es compatible con la reutilización documentada de conglomerados; no acredita seguimiento del mismo hogar ni probabilidades conjuntas de selección exactas. La sensibilidad bajo independencia omite las covarianzas entre olas.

Se prueban tres especificaciones operativas de estratificación: DOMINIO×ESTRATO (principal, para continuidad con el análisis), departamento×ESTRATO y ESTRATO. Son variantes del análisis, no una certificación de que las tres reproduzcan el estrato final del INEI. La documentación disponible no permite certificar todavía esa combinación final. La opción de una UPM por estrato es adjust; se usa Taylor y distribución t. En los contrastes de razones se adopta el menor grado de libertad de las olas implicadas.

## Cómo se considera la estacionalidad

Se añade un modelo lineal ponderado por indicador y territorio:

    indicador = intercepto + ola T1 2026 + ola T2 2026 + mes de entrevista.

Los doce meses entran como categorías; 2025 es la referencia. Se normalizan los pesos a la misma masa total por ola para el modelo, sin sustituir los factores de las estimaciones descriptivas. La inferencia usa el sandwich robusto de encuesta, estratos y UPM. En proporciones se aplica un modelo lineal de probabilidad para expresar el coeficiente en puntos porcentuales; no se extrapolan probabilidades individuales. Los dos coeficientes trimestrales y su promedio se acompañan de IC y ajuste de Holm en familias de 100 y 50, respectivamente.

Este ajuste evalúa diferencias condicionadas al mes observado, bajo el supuesto aditivo del modelo. No crea estimaciones trimestrales oficiales de 2025 ni elimina todo posible efecto estacional. La información de julio–diciembre 2026 no existe en estos archivos; no se extrapola. Un único ciclo anual permite evaluar heterogeneidad mensual, no demostrar una estacionalidad persistente ni atribuir variaciones a lluvia, clima o dengue.

También se calcula una prueba de Wald de heterogeneidad entre meses de 2025 por indicador y territorio (50 pruebas, Holm). Una diferencia mensual puede reflejar composición de muestra o cambios reales; no identifica su causa.

## Verificación y criterio de robustez

R 4.3.3, survey 4.2.1 y data.table 1.14.10 producen 450 contrastes de razones y 450 coeficientes/contrastes ajustados: tres variantes de estratos por 150 comparaciones. Una implementación separada en NumPy reproduce las razones ponderadas, la linealización por conglomerado y el sandwich de mínimos cuadrados ponderados directamente desde los archivos suministrados. Verifica los 900 puntos y errores estándar: diferencia máxima en puntos de 8,58×10⁻¹³ y en EE de 8,09×10⁻¹⁴, con tolerancia de 10⁻⁹. Coincidir numéricamente no certifica el diseño oficial, el supuesto longitudinal ni causalidad.

Se denomina “robusta en las sensibilidades aplicadas” a una diferencia que conserva dirección y Holm<0,05 en las tres variantes de estratos, tanto en los contrastes brutos con UPM compartidas y bajo independencia como en el ajuste por mes con UPM compartidas. No implica que se hayan probado todos los diseños o mecanismos de dependencia posibles. No se ha supuesto independencia para afirmar significación por defecto.

## Resultados principales

En la especificación principal, 11 de los 50 contrastes brutos del promedio T1/T2 frente al anual superan Holm; 8 de los 50 contrastes ajustados por mes lo hacen. Cuatro contrastes cumplen el criterio completo de robustez definido arriba. Estos recuentos incluyen los mismos indicadores para ámbitos relacionados y no representan cuatro hallazgos independientes.

| Resto de la Costa | Anual 2025 | Promedio T1/T2 2026 | Diferencia bruta, IC95 % | Ajustada por mes, IC95 % | Holm ajustado por mes |
|---|---:|---:|---:|---:|---:|
| Agua 24 horas todos los días, entre hogares con red y respuesta válida | 22,1 % | 30,0 % | +7,9 pp [5,3; 10,6] | +7,1 pp [4,7; 9,4] | 2,39×10⁻⁷ |
| Horas diarias equivalentes, entre hogares con red y datos válidos | 8,7 h | 10,6 h | +1,90 h [1,33; 2,47] | +1,87 h [1,38; 2,36] | 5,04×10⁻¹² |

Ambas diferencias cumplen el criterio completo de robustez. Hay mayor continuidad en los periodos observados de 2026 frente a la referencia anual de 2025, pero persiste una brecha importante: el promedio de continuidad completa en el Resto de la Costa sigue alrededor del 30 %. No se concluye una mejora de todo el año 2026 ni un efecto de políticas o intervenciones.

En Lima Metropolitana, el contraste bruto del promedio muestra +6,2 pp en continuidad, pero tras el ajuste por mes no supera Holm (p ajustado≈0,087). Las horas tampoco superan Holm tras ese ajuste (≈0,190). No se presentan como cambios robustos.

Los otros dos contrastes robustos corresponden al indicador P110E=2–4 (origen de la muestra de agua en depósito), en el Resto de la Costa y en el total nacional. Este indicador no equivale a presencia del vector ni al porcentaje de todos los hogares que almacenan agua.

## Cómo repetir

    Rscript contrastes_temporales.R /ruta/enaho /ruta/resultados
    python3 verificar_temporales.py --originales /ruta/enaho --resultados /ruta/resultados
    python3 integrar_temporales.py --resultados /ruta/resultados --output /ruta/datos

El directorio ENAHO debe contener 1031-Modulo01, 1031-Modulo04, 1039-Modulo01, 1039-Modulo04, 1048-Modulo01 y 1048-Modulo04, con sus CSV originales completos o los extractos cotejados. La búsqueda admite subdirectorios. No se publican filas individuales, identificadores de UPM ni hogares/personas en las salidas.

Documentación del cálculo: https://r-survey.r-forge.r-project.org/pkgdown/docs/reference/svyglm.html y https://r-survey.r-forge.r-project.org/pkgdown/docs/reference/regTermTest.html. Factores: diccionarios INEI suministrados, anual 2025 y trimestral 2026 (módulo 100, p. 26; módulo 400, p. 134). El cotejo de originales se conserva en COTEJO_DOCUMENTAL.md y cotejo_originales.json.