Comparación de cuatro modelos de regresión para el estudio de factores de riesgo de hato para datos binomiales correlacionados / Comparison of four regression models for the study of herd risk factors for binomial correlated data

ABSTRACT

The objectives of this study were to compare four linear or logistic regression models and to determine their effects on the level of significance and parameter estimates, using the data from a study on seroprevalence of brucellosis in goats. Information on 5 114 does tested during 2002-2003 from 79 herds in the Bajio region in Michoacan, Mexico was used. The models were the prevalence of seropositive animals per herd (V1), analyzed by a general linear model (GLM), herds with at least one seropositive animal, analyzed by standard logistic regression (SLRH); V1 analyzed by standard logistic regression (SLR), assuming independence among results within a same herd (SLRA); and V1 analyzed by mixed LR, considering the herd as random effect (MLR). The risk factors included in the four models were the presence of abortions the year previous to the study, cleanness of the corral (hygiene) and length of lactation. The V1 variable transformed to arcsine-square root did not show a normal distribution. SLRH model (SLR assuming the herd as the unit of interest) and MLR were not compared because they were not nested models. MLR model adjusted the data better than the SLRA model. The deviance (-2LL) from model SLRH (70.6) was similar to their degrees of freedom (75), suggesting that the model adjusted the data very well. Levels of significance for the risk factors were different, depending of the model used. GLM and SLRH models showed significant effects (P < 0.02) only for the presence of abortions; SLRA model showed significant effect (P < 0.05) for the three risk factors, and MLR, effects of the presence of abortions and lactation length, but not for hygiene. The values for the odd ratios (OR) for the SLRA and MLR models were different; the narrowest confidence intervals corresponded to the SLRA model, and the widest to the SLRH model.

RESUMEN

Los objetivos de este estudio fueron comparar cuatro modelos de regresión lineal o logística y determinar sus efectos sobre los niveles de significancia y parámetros, utilizando los datos de un estudio de seroprevalencia de brucelosis en cabras. Se utilizó la información de 5 114 cabras en 79 hatos de la región del Bajío, en Michoacán, México, durante 2002-2003. Los modelos fueron la prevalencia de animales seropositivos por hato (V1), analizados mediante un modelo lineal general (MLG); hatos con al menos un animal seropositivo, analizados mediante regresión logística estándar (RLEH), V1 analizada mediante RLE, suponiendo independencia entre resultados dentro de un mismo hato (RLEA) y V1 analizada mediante RL mixta, considerando al hato como efecto aleatorio (RLM). Los factores de riesgo incluidos en los cuatro modelos fueron presencia de abortos el año anterior al estudio, limpieza del corral (higiene) y duración de la lactancia. La variable V1 transformada a arcoseno-raíz cuadrada no mostró distribución normal. El modelo RLEH (RL estándar considerando al hato como la unidad de interés) y RLM no se compararon por no ser modelos anidados. El modelo RLM ajustó mejor los datos que el modelo RLEA. La deviance (-2LL) del modelo RLEH (70.6) fue similar a sus grados de libertad (75), ello sugiere que este modelo ajustó estadísticamente bien los datos. Se encontraron niveles de significancia diferentes para los factores de riesgo, según el modelo estadístico utilizado. Los modelos MLG y RLEH mostraron efectos significativos (P < 0.02) sólo de presencia de abortos; el modelo RLEA mostró efecto significativo (P < 0.05) para los tres factores de riesgo, y el RLM, efectos de presencia de abortos y duración de la lactancia, pero no de higiene. Los valores de la razón de momios (OR) para los modelos RLEA y RLM fueron diferentes; los intervalos de confianza más estrechos correspondieron al modelo RLEA, y los más amplios, al RLEH.