SISTEMA ESTADÍSTICO PARA LA PLANEACIÓN Y
CONTROL DE LA PRODUCCIÓN DE CHAYOTE DE
EXPORTACIÓN EN LA REGIÓN DE LAS ALTAS
MONTAÑAS DE VERACRUZ

Statistical System for Planning and Control of Export Chayote

Production in the Altas Montañas Region of Veracruz

María Bárbara Del Valle Juárez

Instituto Tecnológico Superior de Huatusco

Rebeca Pulido Saavedra

Instituto Tecnológico Superior de Huatusco

Lenin Jacobo Rosas Ortiz

Instituto Tecnológico Superior de Huatusco

Mario Ángel Páez Bulbarela

Instituto Tecnológico Superior de Huatusco
pág. 7633
DOI:
https://doi.org/10.37811/cl_rcm.v10i5.26631
Sistema estadístico para la planeación y control de la producción de chayote
de exportación en la región de las Altas Montañas de Veracruz

María Bárbara Del Valle Juárez
1
mdelvallej@huatusco.tecnm.mx

https://orcid.org/0000-0001-5504-5984

Tecnológico Nacional de México / Instituto

Tecnológico Superior de Huatusco

México.

Rebeca Pulido Saavedra

213z0193@alum.huatusco.tecnm.mx

https://orcid.org/0009-0003-0441-1261

Tecnológico Nacional de México / Instituto

Tecnológico Superior de Huatusco

México.

Lenin Jacobo Rosas Ortiz

lrosaso@huatusco.tecnm.mx

https://orcid.org/0000-0003-4974-5866

Tecnológico Nacional de México / Instituto

Tecnológico Superior de Huatusco

México.

Mario Ángel Páez Bulbarela

mpaezb@huatusco.tecnm.mx

https://orcid.org/0009-0008-3989-0158

Tecnológico Nacional de México / Instituto

Tecnológico Superior de Huatusco

México.

RESUMEN

La variabilidad de los volúmenes de producción y de la proporción de producto con calidad exportable
constituye una limitante para la planeación y el cumplimiento de los requerimientos comerciales del
chayote (Sechium edule) en la región de las Altas Montañas de Veracruz.
El objetivo de esta
investigación fue desarrollar y evaluar un sistema estadístico de apoyo para la planeación y control de
la producción de chayote de exportación, mediante la estructuración de registros históricos, el análisis
de la variabilidad productiva y la evaluación del desempeño de los pronósticos.
Se realizó un estudio
cuantitativo, aplicado, no experimental y longitudinal-retrospectivo. La metodología se desarrolló en
cuatro fases: recolección, diagnóstico y estructuración de la información; análisis exploratorio de datos;
evaluación del pronóstico; e integración de indicadores para el control y apoyo a la toma de decisiones.
La base analizada comprendió 38 registros, de los cuales 37 fueron válidos para la evaluación porcentual
del pronóstico. Los resultados evidenciaron una elevada heterogeneidad productiva, con coeficientes de
variación mensuales superiores al 70 % y distribuciones no normales. Se identificó una asociación
positiva muy alta entre la producción real y pronosticada (ρ = 0.9766; p < 0.001). En términos globales,
se registraron 22,057 unidades reales frente a 22,163 pronosticadas, con una desviación de −106
unidades, equivalente a −0.48 % respecto al pronóstico. El desempeño presentó un MAPE de 15.86 %,
una mediana del error porcentual absoluto de 10.57 %, un Bias porcentual medio de 4.88 % y un
índice
de cumplimiento del pronóstico
agregado de 99.52 %. No se encontró asociación estadísticamente
significativa entre la producción real o el error del pronóstico y la calidad exportable. Los indicadores
fueron integrados en un tablero de control en Microsoft Excel, permitiendo identificar desviaciones por
productor y periodo. Se concluye que la integración de registros históricos e indicadores estadísticos
constituye una herramienta útil para fortalecer la planeación y el seguimiento de la producción, aunque
se requiere ampliar la serie histórica para validar formalmente modelos predictivos de series de tiempo.

Palabras clave:
chayote; pronóstico estadístico; planeación de la producción; control de la producción;
agricultura de exportación.

1
Autor principal
Correspondencia:
mdelvallej@huatusco.tecnm.mx
pág. 7634
Statistical System for Planning and Control of Export Chayote Production

in the Altas Montañas Region of Veracruz

ABSTRACT

The variability in production volumes and in the proportion of export
-quality products represents a
limitation to production planning and the fulfillment of commercial requirements for chayote (
Sechium
edule
) in the High Mountains region of Veracruz, Mexico. The objective of this research was to develop
and evaluate a statistical system to support the planning and control of export chayote production
by
structuring historical records, analyzing production variability, and evaluating forecasting performance.

A quantitative, applied, non
-experimental, longitudinal-retrospective study was conducted. The
methodology was developed in four phases: data col
lection, diagnosis, and structuring; exploratory data
analysis; forecast evaluation; and integration of indicators for production control and decision
-making
support. The analyzed database comprised 38 records, of which 37 were valid for percentage
-based
f
orecast evaluation. The results revealed high production heterogeneity, with monthly coefficients of
variation exceeding 70% and non
-normal distributions. A very strong positive association was identified
between actual and forecasted production (
ρ = 0.9766; p < 0.001). Overall, 22,057 actual units were
recorded compared with 22,163 forecasted units, resulting in a deviation of −106 units, equivalent to

−0.48% relative to the forecast. Forecast performance showed a MAPE of 15.86%, a median absolute

percenta
ge error of 10.57%, a mean percentage bias of 4.88%, and an aggregate forecast fulfillment
index of
99.52%. No statistically significant association was found between actual production or forecast
error and export
-quality product. The indicators were integrated into a Microsoft Excel control
dashboard, enabling the identification of deviations by produc
er and period. It is concluded that the
integration of historical records and statistical indicators provides a useful tool for strengthening

production plan
ning and monitoring; however, a longer historical time series is required to formally
validate predictive time
-series models.
Keywords:
chayote; statistical forecasting; production planning; production control; export agriculture.
Artículo recibido 20 agosto 2026

Aceptado para publicación: 20 setiembre 2026
pág. 7635
INTRODUCCIÓN

La planeación de la producción agrícola requiere disponer de información confiable que permita
anticipar los volúmenes disponibles y coordinar la oferta con los requerimientos del mercado. Esta
necesidad adquiere mayor importancia en productos perecederos destinados a exportación, debido a que
las variaciones en los volúmenes cosechados y en la calidad pueden afectar el cumplimiento de los
compromisos comerciales. En la región de las Altas Montañas de Veracruz, la producción de chayote
(Sechium edule) constituye una actividad de relevancia económica en la que participan pequeños y
medianos productores vinculados con empresas comercializadoras. Las características propias del
cultivo, las diferencias entre unidades productivas y la variabilidad temporal hacen necesario fortalecer
los mecanismos utilizados para planear y controlar la producción.

En Agro Cuautlapan Produce S.P.R. de R.L., unidad de estudio de esta investigación, se identificaron
limitaciones relacionadas con la generación, registro y aprovechamiento de la información productiva.
Los registros presentaban falta de estandarización, uso de formatos físicos, información incompleta,
ausencia de una base centralizada y dependencia del conocimiento empírico. Estas condiciones
dificultaban la integración y consulta sistemática de los datos y, principalmente, su transformación en
información cuantitativa para comparar los volúmenes estimados con los realmente obtenidos. Por tanto,
el problema de investigación no radica únicamente en la disponibilidad de registros, sino en la limitada
capacidad para estructurarlos, analizarlos y utilizarlos como apoyo para anticipar la producción e
identificar desviaciones.

Esta problemática resulta relevante porque las diferencias entre la producción pronosticada y la real
pueden repercutir directamente en la planeación comercial. Una sobreestimación puede generar
compromisos superiores a la disponibilidad efectiva del producto, mientras que una subestimación
dificulta anticipar excedentes y organizar su comercialización. En consecuencia, además de generar
estimaciones, resulta necesario evaluar sistemáticamente su precisión y dirección mediante indicadores
que permitan reconocer desviaciones por productor y periodo. Este enfoque adquiere particular
importancia ante la heterogeneidad observada en los registros productivos, por lo que una evaluación
basada exclusivamente en valores agregados puede ocultar comportamientos individuales relevantes.
pág. 7636
Desde el enfoque estadístico, los métodos de pronóstico permiten utilizar información histórica para
generar estimaciones que apoyen la planeación. Los modelos de series de tiempo, constituyen
alternativas para analizar tendencia y estacionalidad (Box et al., 2015; Hyndman & Athanasopoulos,
2018). Su aplicación en el sector agrícola ha sido documentada en estudios sobre diferentes cultivos,
mostrando la utilidad del análisis de series temporales para apoyar la estimación de la producción
(Sernaqué et al., 2022; Krishna Priya & Kausalya, 2024; Kusuma, 2024; Maheshnath et al., 2024). De
manera complementaria, el control estadístico proporciona herramientas para analizar la variabilidad y
establecer mecanismos sistemáticos de seguimiento del desempeño de los procesos (Montgomery,
2019).

En relación con el chayote, investigaciones recientes han abordado diferentes aspectos del cultivo en
Veracruz. Castillo-Martínez y Romero-Jaen (2024) destacan su importancia en la región de las Altas
Montañas; Ruiz-García et al. (2025) analizan factores ambientales y productivos asociados con su
cultivo en Coscomatepec; mientras que Herrera-Ubaldo et al. (2026) estudian la producción,
comercialización y los modelos de negocio relacionados con el sector chayotero veracruzano. Estos
antecedentes contribuyen a comprender las características productivas y comerciales del cultivo; sin
embargo, persiste la necesidad de integrar los registros históricos de producción con indicadores
cuantitativos que permitan evaluar el comportamiento de los pronósticos y facilitar el seguimiento
operativo de la producción destinada a exportación.

En este sentido, el aporte de la investigación consiste en estructurar y analizar la información histórica
disponible y evaluar el desempeño de los pronósticos mediante indicadores como el
Error Porcentual
Absoluto Medio (MAPE), Bias e índice de cumplimiento del pronóstico (ICP).
La incorporación
conjunta de estos indicadores permite analizar la magnitud y dirección de los errores, así como el nivel
de cumplimiento respecto a los volúmenes estimados. Los resultados se integran en un tablero de control
desarrollado en Microsoft Excel para facilitar la visualización de las desviaciones por productor y
periodo y proporcionar información objetiva que apoye la toma de decisiones. Este planteamiento
permite pasar de registros operativos dispersos a un esquema sistemático de análisis, evaluación y
seguimiento de la producción.
pág. 7637
Por lo anterior, el objetivo de esta investigación fue desarrollar y evaluar un sistema estadístico de
apoyo
para la planeación y control de la producción de chayote de exportación en la región de las Altas
Montañas de Veracruz, mediante la estructuración de registros históricos, el análisis de la variabilidad
productiva, la evaluación del desempeño de los pronósticos y la integración de indicadores para apoyar
la toma de decisiones.

METODOLOGÍA

La metodología se desarrolló en cuatro fases secuenciales: recolección, diagnóstico y estructuración de
la información; análisis exploratorio de datos; evaluación del pronóstico; e integración de indicadores
para el control y apoyo a la toma de decisiones.

Fase 1. Recolección, diagnóstico y estructuración de la información

Se recopilaron e integraron los registros históricos y operativos relacionados con la producción de
chayote, considerando variables como productor, periodo de corte, superficie cultivada, volumen de
producción y clasificación del producto. Como etapa inicial se realizó un diagnóstico de la calidad de la
información mediante un diagrama de causa–efecto (Ishikawa), con el propósito de identificar los
factores que podían generar inconsistencias, variabilidad, omisiones o pérdida de datos durante su
registro.

Posteriormente, la información fue sometida a un proceso de depuración, codificación, validación y
estructuración, eliminando registros inconsistentes y homogeneizando criterios de captura. Como
resultado de esta fase se obtuvo una base de datos organizada y confiable para sustentar los análisis
estadísticos posteriores.

Fase 2. Análisis exploratorio de datos

A partir de la base de datos estructurada se realizó el análisis exploratorio de la producción, con el
propósito de caracterizar su comportamiento e identificar patrones relevantes entre productores y
periodos. Se empleó estadística descriptiva mediante medidas de tendencia central y dispersión,
complementada con el análisis de la variabilidad de los volúmenes registrados.

Asimismo, se analizaron las diferencias en los niveles de producción entre productores, la distribución
temporal de los volúmenes y la relación entre la superficie cultivada y la producción obtenida. Este
pág. 7638
análisis permitió identificar productores con mayor participación, periodos con concentraciones
elevadas de producción y diferencias importantes en la escala productiva.

Los resultados de esta fase permitieron establecer la estructura y comportamiento de los datos que
posteriormente fueron utilizados para evaluar los pronósticos de producción.

Fase 3. Evaluación del pronóstico y análisis de su desempeño

Con base en la información histórica previamente estructurada y analizada, se evaluó el componente de
pronóstico mediante la comparación de los volúmenes estimados con los valores reales registrados por
productor y periodo
.
El análisis se concentró en evaluar el desempeño de los pronósticos operativos disponibles mediante la
comparación de los valores estimados y observados.

Para evaluar el desempeño del sistema se utilizaron cuatro indicadores principales:

•
MAPE (Error Porcentual Absoluto Medio): cuantificó la magnitud promedio del error de
pronóstico y permitió comparar la precisión obtenida entre productores.

•
Bias (sesgo): permitió determinar la tendencia sistemática del pronóstico hacia la
sobreestimación o subestimación de la producción.

•
índice de cumplimiento del pronóstico: permitió evaluar el nivel de cumplimiento de la
producción respecto a los requerimientos establecidos.

•
Proporción no exportable: permitió cuantificar la proporción del producto que no fue
clasificado con calidad exportable.

El cálculo de estos indicadores por productor permitió identificar diferencias en la precisión de las
estimaciones y reconocer aquellos casos que presentaron mayores desviaciones respecto al
comportamiento esperado.

Fase 4. Control, visualización y apoyo a la toma de decisiones

Finalmente, los resultados obtenidos fueron integrados en un esquema de control y seguimiento
operativo. Como complemento
se planteó un esquema de seguimiento basado en principios de Control
Estadístico de Procesos (SPC)
para analizar la estabilidad y variabilidad asociada con las etapas de
clasificación y empaque del chayote, particularmente aquellas relacionadas con la obtención de producto
con calidad de exportación.
pág. 7639
La información histórica, la producción pronosticada, la producción real, las desviaciones y los
indicadores MAPE, Bias, índice de cumplimiento del pronóstico y proporción no exportable fueron
integrados en un tablero de control desarrollado en Microsoft Excel.
Esta herramienta permitió
visualizar el comportamiento de los indicadores por productor y periodo, facilitando la identificación de
desviaciones y áreas críticas.

RESULTADOS Y DISCUSIÓN

Fase 1. Recolección, diagnóstico y estructuración de la información

Se integraron los registros históricos y operativos de producción, pronóstico y clasificación del chayote.
Antes del análisis se revisaron las condiciones de generación y registro de la información mediante un
diagrama de causa–efecto (Ishikawa), identificando riesgos asociados con ausencia de criterios formales
de medición, falta de estandarización, registros incompletos, uso de formatos físicos, rotación de
personal, dependencia del conocimiento empírico, ausencia de una base centralizada y limitaciones
tecnológicas. Posteriormente se efectuó la depuración, codificación, validación y estructuración de la
base de datos.

Tabla 1. Criterios de validación de la base utilizada en el análisis

Criterio
Resultado verificado Tratamiento metodológico
Registros totales
38 Se conservaron para totales de producción.
Registros válidos para
MAPE/Bias

37
Se excluyó el caso con producción real = 0.
Registros válidos para

ICP individual

37
Se excluyó el caso con pronóstico = 0.
Caso Real = 0 y

Pronóstico = 0

Pánfilo González, octubre

No se calcula porcentaje; el registro no altera
los totales.

Variable de calidad
Porcentaje exportable Se conserva como indicador de clasificación.
Proporción no

exportable derivada

100
− calidad exportable
Debe interpretarse como proporción no
exportable, salvo evidencia de descarte real.
pág. 7640
Figura 1. Diagrama de Ishikawa utilizado para identificar limitaciones en la confiabilidad y
estructuración de los datos.

Fase 2. Análisis exploratorio de datos

Una vez validada la base, se caracterizó el comportamiento de la producción mediante estadística
descriptiva, medidas de dispersión y análisis de distribución. Se calcularon media, mediana, desviación
estándar y coeficiente de variación (CV). La normalidad se evaluó con Shapiro–Wilk. Debido a que las
distribuciones mensuales resultaron no normales (p < 0.05), las asociaciones se complementaron con el
coeficiente de correlación de Spearman, más robusto ante asimetría y valores extremos.

Tabla 2. Estadística descriptiva de la producción real por mes

Periodo
n Total Media Mediana
Desv.

estándar

CV (%)

Agosto
8 5,282 660.25 514.00 465.75 70.54
Septiembre
8 5,590 698.75 437.50 546.32 78.19
Octubre
21* 11,185 532.62** 405.00** 426.97** 80.16**
Global
37* 22,057 596.14** 445.00** 455.08** 76.34**
Fuente: elaboración propia con base en los registros validados del estudio.

* Para el análisis descriptivo inferencial se excluyó el registro 0/0 por corresponder a una ausencia de producción y
pronóstico. ** Estadísticos calculados sobre las 21 observaciones positivas de octubre y 37 observaciones positivas globales.
pág. 7641
Tabla 3. Pruebas de normalidad y correlación

Análisis
Estadístico p-valor Interpretación
Shapiro
–Wilk, agosto W = 0.8083 0.0351 No normal (α = 0.05)
Shapiro
– Wilk,
septiembre

W = 0.8114
0.0378 No normal (α = 0.05)
Shapiro
–Wilk, octubre W = 0.8716 0.0102 No normal (α = 0.05)
Shapiro
–Wilk, global W = 0.8716 0.0005 No normal (α = 0.05)
Spearman Real
–
Pronóstico

ρ
= 0.9766 < 0.001
Asociación positiva

muy alta

Spearman Real
–
Calidad

ρ
= −0.0045 0.9791
Sin asociación

monotónica

Spearman Error
–
Calidad

ρ
= −0.0297 0.8615
Sin asociación

monotónica

Figura 2. Variabilidad de la producción real por periodo.

Los coeficientes de variación elevados muestran heterogeneidad productiva. Además, la prueba de
Shapiro–Wilk indica que no es apropiado asumir normalidad. La correlación de Spearman confirma que
los pronósticos siguen estrechamente el orden relativo de la producción real, pero la base analizada no
demuestra una relación estadísticamente significativa entre volumen real o error de pronóstico y el
porcentaje de calidad exportable.
pág. 7642
Fase 3. Evaluación del pronóstico y análisis de su desempeño

A partir de la información estructurada se evaluó el desempeño de los pronósticos operativos mediante
la comparación entre los valores de producción real y producción pronosticada por productor y periodo.
Debido a que la base disponible contiene los valores finales reales y pronosticados, el análisis se
concentró en validar las métricas de desempeño de dichos pronósticos y no en identificar o estimar
formalmente el modelo estadístico que los generó. En consecuencia, los resultados obtenidos
corresponden a una evaluación del desempeño predictivo de los registros disponibles.

Fórmulas empleadas: Error = Real − Pronóstico; APE = |Real − Pronóstico| / Real × 100; MAPE =
promedio de los APE válidos; Bias porcentual = (Real − Pronóstico) / Real × 100;
índice de
cumplimiento del pronóstico
= Real / Pronóstico × 100.
Tabla 4. Desempeño agregado del pronóstico por periodo

Periodo

n

válid

o

Real

total

Pronóstic

o total

Desviació

n

Desv

. %

vs

pron

.

MAP

E %

Median

a APE

%

Bias

medi

o %

ICP

(%)

Agosto
8 5,282 4,620 662
14.3

3

12.05
10.12 11.07
114.3

3

Septiembr

e

8
5,590 5,973 −383
−

6.41

12.74
9.38 −4.21 93.59
Octubre
21
11,18

5

11,570
−385
−

3.33

18.50
14.11 5.99 96.67
Total
37
22,05

7

22,163
−106
−

0.48

15.86
10.57 4.88 99.52
pág. 7643
Figura 3. Comparación de producción real y pronosticada por mes.

Figura 4. Error porcentual absoluto por productor y periodo.

Fuente: elaboración propia con base en los datos analizados, salvo indicación en la figura.

La revisión aritmética confirma que el MAPE global correcto es 15.86 % y que la mediana del error
porcentual absoluto es 10.57 %. El
índice de cumplimiento del pronóstico es 99.52 %, resultado de
22,057 unidades reales frente a 22,163 pronosticadas. La desviación agregada es −106 unidades,
equivalente a −0.48 % respecto al pronóstico. Agosto presentó subestimación agregada (14.33 % por
encima del pronóstico), mientras que septiembre y octubre presentaron producción agregada inferior a
lo pronosticado.
pág. 7644
Valor extremo validado. Miguel Bautista en octubre registra producción real de 80 y pronóstico de 7.
El APE es 91.25 % y el
índice de cumplimiento del pronóstico 1,142.86 %. Estos valores son
matemáticamente correctos y deben tratarse como un caso extremo de fuerte subestimación, no como
un error de fórmula.

Fase 4. Control, visualización y apoyo a la toma de decisiones

Los indicadores obtenidos se integraron en un tablero de control en Microsoft Excel para facilitar el
seguimiento por productor y periodo. El control se orienta a detectar desviaciones de pronóstico y
cambios en la proporción exportable. Para el componente de Control Estadístico de Procesos (SPC), la
base disponible permite plantear el seguimiento mediante gráficas por atributos para la fracción
exportable y mediante gráficas por variables para volúmenes, pero no contiene subgrupos temporales
suficientes para recalcular límites de control X̄–R o p de manera estadísticamente defendible en esta
revisión. Por tanto, estos límites deben incorporarse únicamente cuando existan observaciones
secuenciales por corte y subgrupo.

Tabla 5. Indicadores clave para el tablero de control

Indicador
Resultado global Criterio de interpretación
MAPE
15.86 %
Error porcentual absoluto
promedio de 37 observaciones
válidas.

Mediana APE
10.57 % La mitad de los errores
individuales es
≤ 10.57 %.
Bias porcentual medio
4.88 %
Sesgo medio hacia
subestimación; sensible a
valores extremos.

índice de cumplimiento del

pronóstico
99.52 %
Cumplimiento agregado
respecto al volumen
pronosticado.

Calidad exportable promedio
30.84 % Proporción media clasificada
como exportable.

Proporción no exportable
69.16 %
100
− calidad exportable; no
equivale necesariamente a
merma física.
pág. 7645
Figura 5.
Índice de cumplimiento del pronóstico (ICP) agregado por mes respecto a la referencia
operativa de 98 %.

Figura 6. Bias porcentual por productor: valores positivos indican subestimación y negativos
sobreestimación.

Con una referencia operativa del
índice de cumplimiento del pronóstico ≥ 98 %, agosto supera el
objetivo y el resultado agregado del periodo también lo cumple; septiembre (93.59 %) y octubre (96.67
%) quedan por debajo. El Bias por productor permite identificar la dirección del error, pero debe
interpretarse junto con la desviación agregada, ya que valores extremos pueden cambiar el promedio del
sesgo porcentual. En particular, octubre muestra desviación agregada negativa (−3.33 % respecto al
pronóstico) aunque su Bias porcentual medio es positivo (5.99 %) debido al efecto de observaciones
extremas.
pág. 7646
D
ISCUSIÓN DE RESULTADOS
Los resultados evidenciaron una elevada heterogeneidad en la producción de chayote entre productores
y periodos, con coeficientes de variación superiores al 70 % y distribuciones no normales. Esto indica
que la planeación no debe sustentarse únicamente en valores promedio, sino considerar las diferencias
entre productores y la variabilidad temporal. Este comportamiento es congruente con lo reportado por
Ruiz-García et al. (2025), quienes destacan la influencia de factores ambientales en el cultivo de Sechium
edule, y por Herrera-Ubaldo et al. (2026), quienes describen diferencias en los sistemas de producción
y comercialización del chayote en Veracruz. Por ello, la variabilidad observada puede estar asociada con
factores productivos, ambientales y de manejo que deberán incorporarse en futuros análisis.

La asociación positiva muy alta entre la producción real y la pronosticada (ρ = 0.9766;
p < 0.001)
muestra que las estimaciones reproducen adecuadamente el comportamiento relativo de los volúmenes

registrados; sin embargo, esta relación no implica precisión absoluta. El MAPE global de 15.86 % y la

mediana del error porcentual absoluto de 10.57
% evidencian que existen diferencias individuales
relevantes, aun cuando la desviación acumulada fue de solo −0.48 %, resultado de 22,057 unidades

reales frente a 22,
163 pronosticadas. Además, agosto presentó producción superior al pronóstico,
mientras que septiembre y octubre registraron valores inferiores, lo que confirma que los resultados

agregados pueden ocultar comportamientos particulares.

La evaluación mediante métricas de error es consistente con estudios de pronóstico agrícola. Sernaqué

et al. (2022) compararon ARIMA y Holt
-Winters para series de producción de cereales; Kusuma (2024)
evaluó Holt
-Winters y Seasonal ARIMA en producción de arroz; Krishna Priya y Kausalya (2024)
utilizaron ARIMAX para variables relacionadas con producción de trigo; y Maheshnath et al. (2024)

aplicaron ARIMA para pronosticar producción de maíz. Estos antecedentes respaldan la necesidad de

evaluar cuantitativame
nte el desempeño de los pronósticos antes de utilizarlos como apoyo para la
planeación. En el presente estudio, la identificación de valores extremos, como el caso registrado en

octubre con un APE de 91.25 %, confirma también la importancia de analizar los
errores por productor
y no únicamente mediante indicadores globales.

Por otra parte, no se encontró asociación estadísticamente significativa entre la producción real o el error

del pronóstico y la calidad exportable. Este resultado sugiere que el volumen producido no explica por
pág. 7647
sí mismo la proporción de producto que cumple con las condiciones de exportación, por lo que será

necesario incorporar variables agronómicas, ambientales, de clasificación y manejo poscosecha. Esta

interpretación coincide con los antecedentes de Ruiz
-García et al. (2025) y Herrera-Ubaldo et al. (2026),
quienes muestran la importancia de considerar factores adicionales en el análisis de la producción

chayotera.

Finalmente, la integración del MAPE, Bias, índice de cumplimiento del pronóstico y proporción
no
exportable en un tablero de control permite transformar los registros históricos en información útil para
el seguimiento y la toma de decisiones. Este enfoque es consistente con Montgomery (2019), quien
señala que el análisis de la variabilidad constituye un elemento fundamental del control estadístico de
procesos. No obstante, la base disponible aún no contiene suficientes observaciones secuenciales para
establecer límites de control formalmente validados, por lo que el sistema debe considerarse actualmente
como una herramienta de monitoreo y apoyo a la toma de decisiones. La ampliación de la serie histórica
permitirá, en investigaciones posteriores, comparar formalmente modelos como ARIMA, Holt-Winters
o ARIMAX mediante procedimientos de validación fuera de muestra y determinar cuál presenta mayor
capacidad predictiva para la producción de chayote de exportación en la región de las Altas Montañas
de Veracruz.

CONCLUSIONES

Los resultados obtenidos permiten concluir que la integración de registros históricos, herramientas

estadísticas e indicadores de desempeño constituye una alternativa viable para fortalecer la planeación

y el control de la producción de chayote de exportac
ión en la región de las Altas Montañas de Veracruz.
El análisis evidenció una elevada heterogeneidad productiva entre productores y periodos, con

coeficientes de variación superiores al 70 % y distribuciones que no presentaron normalidad. Esta

condición de
muestra que la planeación no debe sustentarse únicamente en valores promedio, sino
considerar el comportamiento particular de los productores y la variabilidad temporal de la producción.

Asimismo, la alta asociación entre la producción real y pronosticada
(ρ = 0.9766; p < 0.001) indica que
las estimaciones reproducen adecuadamente el comportamiento relativo de los volúmenes registrados,

aunque la presencia de errores individuales confirma la necesidad de mantener un seguimiento

desagregado.
pág. 7648
La evaluación del desempeño mostró un MAPE global de 15.86 %, una mediana del error porcentual

absoluto de 10.57 % y un índice agregado de cumplimiento del pronóstico de 99.52 %. Estos resultados

permiten considerar al pronóstico como un apoyo cuantitativo
para la planeación; sin embargo, su
interpretación debe realizarse conjuntamente con las desviaciones por periodo y productor. La diferencia

agregada de −106 unidades entre la producción real y la pronosticada representa apenas −0.48 %

respecto al volumen
estimado, pero este equilibrio global oculta comportamientos distintos: agosto
presentó una producción superior a la pronosticada, mientras que septiembre y octubre registraron

volúmenes inferiores a los previstos. Por ello, el resultado agregado no debe
asumirse como evidencia
suficiente de una precisión uniforme entre productores y periodos.

Por otra parte, no se encontró asociación estadísticamente significativa entre la producción real o el error

de pronóstico y el porcentaje de calidad exportable, lo que indica que el volumen producido no explica

por sí mismo la proporción de producto que c
umple con las condiciones requeridas para exportación.
En consecuencia, la calidad debe analizarse mediante variables adicionales relacionadas con las

condiciones productivas, agronómicas, ambientales, de clasificación y de manejo poscosecha. Esta

necesida
d adquiere relevancia considerando la diversidad de líneas de investigación existentes alrededor
de
Sechium edule y la necesidad de continuar ampliando el conocimiento sobre los factores que
intervienen en su producción y aprovechamiento (Barrera
-Guzmán et al., 2024). De igual forma, la
proporción no exportable no debe interpretarse directamente como merma física mie
ntras no exista
evidencia que permita determinar el destino final del producto que no alcanza la clasificación de

exportación.

En función de lo anterior, el sistema estadístico propuesto cumple su propósito como herramienta de

apoyo
para la toma de decisiones, al transformar registros históricos en información estructurada,
evaluar el desempeño de los pronósticos e integrar indicadores para su seguimiento mediante un tablero

de control. No obstante, los resultados deben circunscribirse a la información analiza
da y no constituyen
todavía una validación formal de un modelo específico de series de tiempo. Como continuidad de la

investigación, se recomienda ampliar la serie histórica y disponer de observaciones secuenciales por

corte que permitan comparar formalmen
te distintas alternativas de pronóstico. En este sentido, existen
antecedentes metodológicos basados en Holt
-Winters y sus variantes (Liu & Wu, 2022; Pleños, 2022;
pág. 7649
Priyanka et al., 2022), así como aproximaciones que comparan o complementan modelos ARIMA con

técnicas de inteligencia artificial (Neog, 2024; Suthar et al., 2023), los cuales pueden constituir

referentes para futuras etapas de validación predictiva. Estas
comparaciones deberán realizarse mediante
procedimientos de entrenamiento y evaluación fuera de muestra para determinar cuál alternativa

presenta mayor capacidad de generalización bajo las características particulares de la producción de

chayote.

Finalmente, futuras investigaciones deberán incorporar variables climáticas, agronómicas y de manejo

poscosecha, evaluar el desempeño del esquema en otros ciclos productivos y unidades de producción y

analizar la estabilidad temporal de los indicadores. La
ampliación de la información permitiría avanzar
de un sistema de evaluación y seguimiento del pronóstico hacia un modelo predictivo formalmente

validado y, posteriormente, establecer límites estadísticos de control sustentados en observaciones

secuenciale
s. De esta manera, la investigación podría contribuir no sólo al seguimiento de la producción,
sino también a una planeación anticipada de los volúmenes destinados a exportación y a una toma de

decisiones sustentada en evidencia cuantitativa

REFERENCIAS BIBLIOGRÁFICAS

Barrera-Guzmán, L. A., Cadena-Iñiguez, J., Mojica-Zárate, H. T., Cisneros-Solano, V. M., & Alejo-
Bello, S. (2024).
Trends and evolution in the scientific research of Sechium edule (Jacq.) Swartz.
Agro Productividad, 17(8), 181
–187. https://doi.org/10.32854/agrop.v17i7.2956
Box, G. E. P., Jenkins, G. M., Reinsel, G. C., & Ljung, G. M. (2015). Time series analysis: Forecasting

and control (5th ed.). John Wiley & Sons.

Castillo-Martínez, S. I., & Romero-Jaen, M. (2024). Chayote: El tesoro verde de las Altas Montañas de
Veracruz–raíces de tradición y sabor. Revista Mexicana de Agroecosistemas, 11(2), 179–185.
https://doi.org/10.60158/rma.v11i2.444

Herrera-Ubaldo, M., Santoyo-Cortés, V. H., Muñoz-Rodríguez, M., & Figueroa-Rodríguez, K. A.
(2026). Diferenciación de la producción y comercialización de chayote en Veracruz y modelos
de negocio asociados. Revista de Geografía Agrícola, (75), e25004.
https://doi.org/10.5154/r.rga.2025.01.004

Hyndman, R. J., & Athanasopoulos, G. (2018). Forecasting: Principles and practice (2nd ed.). OTexts.
pág. 7650
Krishna Priya, S. R., & Kausalya, N. (2024). Forecasting wheat crop area, production & productivity in

India using ARIMAX model. Journal of Agriculture, Biology and Applied Statistics, 3(2), 57
–
68. https://doi.org/10.47509/JABAS.2024.v03i02.02

Kusuma, M. (2024). Perbandingan metode Holt
-Winters dan Seasonal ARIMA dalam memprediksi
produksi beras di Indonesia. Jurnal Agribisnis, 26(2). https://doi.org/10.31849/agr.v26i02.17672

Liu, L., & Wu, L. (2022). Holt
-Winters model with grey generating operator and its application.
Communications in Statistics
—Theory and Methods, 51(11), 3542–3555.
https://doi.org/10.1080/03610926.2020.1797804

Maheshnath, M., Vijaya Kumari, R., Suhasini, K., Srinivasa Reddy, D., & Meena, A. (2024). Forecasting

maize production in Telangana State using ARIMA model. Archives of Current Research

International, 24(6), 223
–229. https://doi.org/10.9734/acri/2024/v24i6780
Montgomery, D. C. (2019). Introduction to statistical quality control (8th ed.). John Wiley & Sons.

Neog, B. (2024). An improved ARIMA model using ANN and SVM for forecasting rapeseed & mustard

production. International Journal of Intelligent Systems and Applications in Engineering,

12(21s), 3152
–3162.
Pleños, M. (2022). Time series forecasting using Holt
-Winters exponential smoothing: Application to
abaca fiber data. Problems of World Agriculture, 22(2), 17
–29.
https://doi.org/10.22630/PRS.2022.22.2.6

Priyanka, P. A., Nandakumar, E., & Sangeetha, R. (2022). Holt
-Winters forecast on seasonal time series
crop insurance data with structural outliers (2000
–2018) using MS Excel. Asian Journal of
Agricultural Extension, Economics & Sociology, 40(9), 136
–145.
https://doi.org/10.9734/AJAEES/2022/v40i930986

Ruiz-García, V. H., Ruiz-García, P., Pérez-Nieto, J., & Monterroso-Rivas, A. I. (2025).
Chayote
(Sechium edule var. virens levis) cultivation and environmental factors.
Agro Productividad,
18(9), 91–99. https://doi.org/10.32854/8xb2vh64

Sernaqué, H., Meca, M., Zapata, E., Marchan, B., Medina, J., Nole, D., Aldana, C., Saavedra, Y., Trelles,
L., Chuquihuanca, N., & Mendoza, G. (2022).
Comparison of ARIMA and Holt-Winters
pág. 7651
forecasting models for time series of cereal production in Peru. Applied Human Factors and

Ergonomics International, 22. https://doi.org/10.54941/ahfe1001007

Suthar, B., Pundir, R. S., & Popat, R. (2023). Groundnut price forecasting in Gondal market of Gujarat:

A comparison of ARIMA and ANN models. Indian Journal of Agricultural Marketing, 37(1),

51
–60. https://doi.org/10.22004/ag.econ.399877