Abstract
Population growth curves in two predator species in the region of Irapuato, Guanajuato, Mexico, were estimated using multivariate multiple regression techniques through third-degree polynomials functions defined by: Y1(T, PP)=1.058e+02-1.300e+01* T-7.817e-01^* PP+4.058e-01* T2+1.037e-02* PP2+3.922e-02*T PP-5.360e-04*T PP2 Y2(T, PP) = 1.246e + 02-1.761e+01*T - 8.349e-01* PP+6.359e - 01*T2 + 6.767e-02* PP2 - 1.166e-02*T PP-2.865e-04*TPP2 where T is the average monthly standard temperature, PP is the average monthly accumulated rainfall and Y1 and Y2 are the numbers of Hippodamia convergens and Olla v-nigrum adults, respectively. The Roy criterion was used in the multivariate analysis of variance of the multivariate multiple regression, which statistical test value was 80.379 with a probability for its corresponding approximation F of 0.0001295, consequently, the estimated polynomials for the population growth curves were adequate. Multivariate association measures Wilks (ALH = 0.987978) and Roy (θ= 0.9877118) were calculated, allowing estimation of growth curves as prediction tools. These models of polynomial growth curves were compared to each other through the statistic W=1978.3, which yielded a probability value P = 0**, concluding that the growth curves were different. The critical point for temperature and rainfall in both species coincided with the point (13.7° C, 2.0 mm), in which the maximum numbers of individuals were obtained, resulting 3.4 to H. convergens and 2.2 for O. v-nigrum.
Introducción
Los coccinélidos se usan ampliamente en el manejo integra- do de plagas, así como en programas aumentativos de control biológico, tal es el caso de Olla v-nigrum (Mulsant, 1866) que depreda al psílido de los cítricos Diaphorina citri (Kuwa- yama, 1908) (Hemiptera: Psyllidae) y debido a su marcada respuesta numérica es un elemento potencial de control bio- lógico (Michaud 2001). En México, en la región del Bajío Guanajuatense, las poblaciones de la catarinita roja Hippo- damia convergens (Guérin-Méneville, 1842) (Coleoptera: Coccinellidae) emigran de las zonas bajas hacia las sierras aledañas donde forman grandes agregados en la base de los árboles, posteriormente retornan a los valles en febrero (Sa- lazar y Salas 1994), igual ocurre en el estado de California, EEUU (Lee 1980); ambas especies de coccinélidos ajustan su ciclo de vida a la disponibilidad de áfidos para su alimen- tación (Nicholls 2008). Algunas otras especies como Adalia bipunctata (Linneo, 1758) (Coleoptera: Coccinellidae) se adaptan a diferentes temperaturas, ya que incrementan el tamaño del cuerpo y reducen la sobrevivencia de las pupas cuando las temperaturas son bajas, siendo más susceptibles a temperaturas mayores (Sorensen et al. 2013).
La fluctuación de las poblaciones insectiles, sus inte- racciones tróficas y la ecología de su comportamiento, se han descrito a través de la creación de modelos ecológicos (Hilker y Westerwoff 2007). Aunque un modelo ecológico refleja las variaciones observadas en la población, no se per- cibirán claramente sin el uso de instrumentos analíticos que permitan estimar y/o predecir esas variaciones, consideran- do además las interacciones depredador-presa o huésped- parásito, que con frecuencia forman parte de los programas de manejo integrado de plagas. Lima et al. (2009) indicaron que una aproximación a los sistemas ecológicos se logra con la simplificación matemática; es decir, la creación de un modelo matemático con una estructura simple como proto- tipo de un sistema agroecológico, que facilite no solamente su utilización práctica, sino también su interpretación y va- lidación.
Venette et al. (2010) afirmaron que para construir mode- los que estimen con buena aproximación a las poblaciones de insectos, es necesario considerar información específica so- bre los factores externos e internos que inciden en su fluctua- ción a través del tiempo. En cuanto a los primeros, Hodgson et al. (2011) señalaron que las variables climáticas influyen en la duración de los ciclos de vida y supervivencia de los individuos, lo que resulta en un incremento o reducción de sus poblaciones.
Se han realizado estudios que utilizan diferentes tipos de modelado en algunas especies de insectos, buscando aproxi- mar el número de individuos considerando factores abió- ticos como variables independientes, como el de Acevedo et al. (2010), que con información de dos años, obtuvieron una ecuación para pronosticar la dinámica poblacional de adultos de Epilachna varivestis (Mulsant, 1850) (Coleopte- ra: Coccinellidae) en términos del tiempo, la edad fisioló- gica y la temperatura, logrando variaciones no mayores al 7% al validar dicho modelo; de igual forma, Ahumada et al. (2004) generaron un modelo de predicción del crecimiento poblacional de Culex quinquefasciatus (Say, 1826) (Diptera: Culicidae), incorporando variables como la temperatura, la precipitación y la altura sobre el nivel del mar, para com- prender cómo estos agentes influyen sobre los cambios po- blacionales de este díptero.
Se han construido modelos utilizando la técnica de la re- gresión con el fin de estimar poblaciones insectiles, cuyos predictores son precisamente variables climáticas. Por ejem- plo, utilizando datos de tres años y mediante un análisis de regresión múltiple, Murphy et al. (2012) determinaron que la distribución espacial de Circulifer tenellus (Baker, 1896) (Hemiptera: Cicadellidae) depende principalmente de las variables temperatura media, precipitación y velocidad del viento. En otra investigación, Villa y Villa (1996) a través de un modelo de regresión múltiple, concluyeron que la tem- peratura y la precipitación explican el 89% de la variación en el número de Dendroctonus mexicanus (Hopkins, 1905) (Coleoptera: Scolytidae), utilizando capturas durante un año. Empleando registros disponibles de año y medio, Cuéllar et al. (2012) encontraron una ecuación de regresión múltiple con un coeficiente de determinación R2 = 0,64, para explicar la relación entre la abundancia de individuos del mismo esco- lítido con respecto a la temperatura y la precipitación pluvial. Por su parte, Gopar y Losada (2004), en dos años de conteos, estimaron la densidad poblacional de cuatro especies de cur- culiónidos en alfalfa, a través de modelos de regresión múlti- ple en función de la temperatura, la precipitación y el fotope- riodo, obteniendo los mayores coeficientes de determinación para las especies Pantomorus auripes (Hustache, 1947) (R2 = 0,56) y Aramigus tessellatus (Say, 1824) (R2 = 0,67).
Existen pocos trabajos sobre modelos que aproximen la abundancia de adultos de coccinélidos para la región de El Bajío, por lo que analizar las variaciones poblacionales de H. convergens y O. v-nigrum a través del tiempo, será un aporte más a considerar en el manejo integrado de áfidos (insectos- presa); por lo anterior, los objetivos de este trabajo fueron: (1) proponer y estimar las funciones de crecimiento pobla- cional de estas dos especies en términos de la temperatura y la precipitación, (2) realizar la comparación entre si de las curvas de crecimiento para las dos especies, y (3) estimar los valores de la temperatura y precipitación, para los que se al- canza el número máximo de individuos adultos.
Materiales y métodos
Resultados y discusión

Temperatura media mensual (°C) en el periodo de 1993 a 2002 en Irapuato, Guanajuato, México.

Precipitación pluvial acumulada mensual (mm) en el periodo 1993 a 2002 en Irapuato, Guanajuato, México.

Fluctuación poblacional de adultos de H. convergens de 1993 a 2002 en Irapuato, Guanajuato, México.

Fluctuación poblacional de adultos de O. v-nigrum de 1993 a 2002 en Irapuato, Guanajuato, México.
Además, al estimar las curvas de crecimiento poblacional, se planteó como hipótesis que los dos modelos deben tener la misma forma funcional, se procedió entonces a la estimación del mejor modelo de regresión de manera univariada para cada variable dependiente (para cada especie), basándose en los valores de la prueba de F de los análisis de varianza de la regresión y los coeficientes de determinación parciales (R2). Una vez definido el modelo de mayor grado necesario para tener el mejor ajuste, representado por el polinomio de tercer grado: Y(T,PP) = b0 + b1*T + b2*PP + b3*Tβ + b4*PP2 + b5*TPP + b6*TPP2, se procedió a estimar los dos poli- nomios que modelaron las curvas de crecimiento poblacional para la cantidad de adultos de las dos especies bajo estudio, en función de los promedios mensuales de temperatura media (T) y precipitación pluvial acumulada (PP) durante diez años a través de las técnicas de regresión múltiple multivariada, las cuales también aplicaron Murphy et al. (2012) y Villa y Villa (1996) pero para el caso univariado, demostrando su eficiencia al estimar densidades poblacionales de insectos adultos.
La prueba utilizada en el análisis de varianza multivaria- do (MANOVA) de la regresión múltiple multivariada, fue el criterio de Roy, cuyo valor del estadístico de prueba, la aproximación de F y su probabilidad P se presentan en la Ta- bla 1, lo que indicó que hay evidencia suficiente en la muestra para rechazar la hipótesis nula, concluyendo que los polino- mios estimados para las curvas de crecimiento poblacional son adecuados, de forma similar a lo obtenido por Cuéllar et al. (2012) quienes, basándose en el coeficiente de deter- minación, lograron estimar en forma adecuada el número de adultos de D. mexicanus, pero mediante técnicas de regresión múltiple univariada.
Utilizando los algoritmos sugeridos por Rencher (2002), se calcularon las medidas de asociación multivariadas θ de Roy y ALH de Lawley-Hotelling, cuyos valores permiten concluir que ambos modelos pueden usarse como herra- mientas de estimación (Tabla 2), estas medidas son análogos multivariados de la medida de bondad de ajuste R2 utilizada en el caso univariado. Al comparar las medidas de asocia- ción obtenidas en este estudio (caso multivariado), con los coeficientes de determinación obtenidos por Gopar y Losada (2004) en cuatro especies de coleópteros curculiónidos (caso univariado), se concluyó que los modelos para estimar las poblaciones de H. convergens y O. v-nigrum tienen mayor grado de precisión, Rosa et al. (2011) indicaron que modelar un sistema biológico es muy difícil debido a la gran canti- dad de factores que intervienen, pero es posible hacer una generalización de los resultados obtenidos, en esta investi- gación se aprecia que el número de individuos capturados es variable, pero su presencia en los diferentes años es constante (Figs. 3 y 4), lo que posiblemente ayudó al ajuste obtenido en los modelos, obsérvese que en algunos años los individuos de O. v-nigrum superaron a los de H. convergens. Además, para una perspectiva univariada para el lector, en la tabla 3 se presentan los polinomios de tercer grado estimados, los valores de F para los análisis de varianza (ANOVA) y los correspondientes valores de R2 y el valor de la probabilidad P para cada especie.
Para comparar los modelos de las curvas de crecimien- to poblacionales polinómicos, se eliminó el parámetro b0 de cada polinomio, ya que cuando se comparan poblaciones de diferentes especies, tan sólo el b0 puede dar como resultado que los modelos sean diferentes, pues este parámetro repre- senta el efecto del número inicial de adultos presentes de cada especie, para esta comparación se utilizó la prueba de hipó- tesis sugerida por Graybill (1976), mediante el estadístico W = 1978,3 con probabilidad P = 0,0, concluyendo que las cur- vas de crecimiento polinómicas poblacionales son diferentes, ésto de alguna manera confirma la diferencia en los hábitos de alimentación de estas especies, O. v-nigrum es de hábi- tos arbóreos donde las condiciones son más estables y hay inóculos de las presas para alimentarse, en este trabajo las poblaciones de esta especie se incrementaron gradualmente a diferencia de H. convergens que mostraron una decremen- to poblacional al inicio y posteriormente se elevaron, ésta es una especie que se desarrolla en hábitats menos estables, ta- les como cultivos bajos o matorrales, los cuales desaparecen anualmente, o los incendios defolian las plantas y matan las presas, sobreviven entonces alimentándose de néctar y polen. Al respecto, Tarango et al. (2001) indicaron que O. v-nigrum es afectada cuando las poblaciones de áfidos son muy ba- jas en las nogaleras, por lo que es conveniente inocular las huertas para incrementar las poblaciones y tener resultados satisfactorios en el control; de la misma forma, esta catarinita es muy escasa en huertas de cítricos cuando la presencia de Diaphorina citri (Kuwayama, 1908) (Heteroptera: Psyllidae) es muy baja, pero se incrementa notablemente cuando sus poblaciones se elevan (Michaud 2001), lo que hace suponer en la presente investigación, que la variación observada en los diferentes años estuvo relacionada con la presencia o au- sencia de áfidos y, de manera indirecta el incremento de tem- peratura y precipitación pluvial elevaron las poblaciones de áfidos, sobre todo en árboles.
MANOVA de la regresión múltiple multivariada aplicando el criterio de Roy.
Medidas de asociación multivariadas de Roy y Lawley-Hotelling.
Modelos de regresión múltiple multivariada para el número de adultos de dos especies de coleópteros (Y), en términos de variables climáticas: temperatura (T) y precipitación (PP), con sus respectivos valores de la prueba de F, probabilidad P y coeficientes de determinación parciales R2.
Por último, en la tabla 4 se muestran los puntos críticos encontrados (T, PP), los cuales establecen las condiciones climatológicas de temperatura y precipitación bajo las cuales se presenta la máxima densidad poblacional de adultos para cada una de las dos especies de coccinélidos, esto se apre- cia gráficamente en las curvas de crecimiento estimadas para cada una de las dos especies estudiadas (Figs. 5 y 6).
Puntos críticos de temperatura y precipitación (T, PP) y valores máximos del número de adultos de dos especies de coccinélidos (Y).
Los resultados obtenidos permiten aproximar la cantidad de estas catarinitas, que como insectos benéficos podrán in- cluirse en la planeación de programas más precisos de ma- nejo integral de plagas, esta técnica la usaron eficientemente Northfield et al. (2008) con poblaciones de Frankliniella spp (Karny, 1910) (Thysanoptera: Thripidae) y la chinche Orius (Wolff, 1811) (Hemiptera: Anthocoridae) para modelar las posibles fuentes de incrementos poblacionales. Asimismo, estos resultados podrán incluirse en las alertas fitosanitarias establecidas en el estado de Guanajuato, considerando las condiciones de temperatura y precipitación para las que los enemigos naturales de los áfidos alcanzan su mayor inciden- cia y evitar en lo posible, la aplicación de insecticidas. La información también puede usarse como referencia en inves- tigaciones posteriores, tendientes a determinar el umbral eco- nómico de los pulgones, basándose en la presencia de estos coccinélidos.

Curva de crecimiento estimada para la H. convergens (Y1) en función de temperatura (T) y precipitación (PP).

Curva de crecimiento estimada para la O. v-nigrum (Y2) en función de temperatura (T) y precipitación (PP).
Conclusiones
Se estimaron las densidades de adultos de H. convergens (Y1) y O. v-nigrum (Y2) en términos de la temperatura y precipita- ción, a través de la obtención de curvas de crecimiento pobla- cionales, utilizando la técnica de regresión no lineal múltiple multivariada y proponiendo como modelos funciones polinó- micas de tercer grado definidas respectivamente por: Y1 (T,PP) = 1,058e+02 - 1,300e+01*T - 7,817e-01*PP + 4,058e-01*T2 + 1,037e-02*PP2 + 3,922e-02*TPP - 5,360e-04*TPP2 Y2 (T,PP) = 1,246e+02 - 1,761e+01*T - 8,349e-01*PP + 6,359e-01*T2+6,767e-02*PP2-1,166e-02*TPP-2,865e-04*TPP2
El rechazo de la prueba de igualdad entre las curvas de crecimiento de ambas especies, posiblemente se debió a que los hábitos de ingesta en un medio efímero en alimento para la primera especie y más estable para la segunda, influyeron en la diferencia.
A partir de estas curvas de crecimiento, se calcularon los valores de 13,7 °C de temperatura y 2,0 mm de precipitación, para los cuales se presentó el máximo número de individuos adultos de ambas especies.
