Livestock Research for Rural Development 18 (4) 2006 | Guidelines to authors | LRRD News | Citation of this paper |
El propósito de este trabajo fue estimar factores multiplicativos de ajuste de la producción de leche a 240 y 270 días de lactación, edad adulta y época de parto en búfalas de la Costa Atlántica colombiana. Fueron empleados para la producción de leche a 240 días 2548 lactancias de 1059 búfalas, para producción de leche a 270 se utilizaron 1736 lactancias de 959 búfalas y para la estimación de los factores a edad adulta se utilizaron 3153 lactancias de 1150 búfalas. Se emplearon los registros de producción de la Hacienda Altamar, ubicada en el Municipio de Buenavista, en el Departamento de Córdoba. Los análisis estadísticos fueron realizados mediante la utilización de un modelo animal y el método de máxima verosimilitud restricta libre de derivadas. Los efectos fijos considerados en el modelo para el ajuste a los 240 y 270 días en leche fueron: grupo contemporáneo (año- época), clases de duración de la lactancia (meses de edad) y la covariable edad de la búfala al parto (efecto lineal y cuadrático). Para el ajuste de producción de leche a edad adulta fueron considerados como efectos fijos: época de parto y clases de edad de la búfala. Los efectos aleatorios para todos los análisis fueron: genético aditivo, de ambiente permanente y residual. Los factores de ajuste para 240 y 270 fueron eficientes con correlaciones entre producciones reales y ajustadas de 0.92. Los factores de ajuste para edad adulta fueron eficientes, ya que removieron la variación causada por el efecto de la edad y la época de parto.
Palabras Claves: Búfalas, factores de ajuste, mejoramiento genético, modelo animal, producción de leche
The purpose of this research was to estimate multiplicative adjustment factors for the milk yield at 240 and 270 days in milk (DIM), age of maturity and season of calving, in buffaloes in the Atlantic Colombian Coast. The number of lactations used for milk yield at 240 DIM was 2548 from 1059 buffaloes and 270 DIM was for 1736 lactations of 959 buffaloes. For effect of age of maturity, 3153 lactations of 1150 buffaloes of the Altamar farm from municipality of Buenavista in the department of Córdoba were used. The statistical analyses were done by the application of the animal model and the derivate-free restricted maximum likelihood method. For the adjustment of milk yield at 240 and 270 DIM, the fixed effects considered in the model were: contemporary group (year-season), and classes of lactation length (months of age) and the variable age of the buffaloes at calving (linear and quadratic effect). For milk yield at mature age the fixed effects considered in the model were: parity (year of calving), age and season of calving and lactation length (months of age). The random effects for all analyses were: genetic additive, permanent environment and residual effects. The adjustment factors for 240 and 270 DIM were significant with correlations between real and adjustment productions of 0.92. The adjustment factors of age of buffalo and parity season were significant.
Key words: Adjustment factors, animal model, buffalos, genetic improvement, milk yield
La selección aplicada sobre una o varias características productivas y reproductivas implica la escogencia de los individuos genéticamente superiores. Sin embargo, existen varios factores ambientales que influyen en el desempeño productivo y reproductivo de los animales (Schmidt y Van Vleck 1974. En virtud de esto, cuando es necesario comparar individuos por medio de registros de producción se requiere corregir los efectos ambientales que alteran el potencial productivo.
Algunos factores ambientales que influyen directamente sobre la producción de leche, pueden ser removidos utilizando el ajuste previo de los datos o pueden ser incluidos en los modelos estadísticos de evaluación genética. Los principales efectos ambientales controlados con factores de ajustes para producción de leche por lactancia, envuelven otras características de desempeño de la vaca como: duración del periodo seco, duración del periodo parto-concepción; efectos causados por el manejo ó nivel de producción de las haciendas como: sistemas de alimentación, sistemas y número de ordeños, efectos causados por el ciclo de vida del animal como: edad y número de partos de la hembra y circunstancias ambientales, como época de parto (Searle 1962; Miller 1973; Schmidt y Van Vleck 1974; Warwick y Legates 1980)
En Colombia no se tienen criterios definidos para la selección de búfalos lecheros y no se cuenta con estudios de curvas de lactancia y factores ambientales que influyen en la producción de leche de los animales, por lo tanto, el propósito de este trabajo fue estimar factores de corrección para la producción de leche de búfalas en la región de la costa atlántica colombiana a los 240 y 270 días de lactancia, edad adulta y época de parto.
Se utilizaron los registros de producción de leche de búfalas de la hacienda Altamar desde el año de 1999 hasta el 2004. El hato que se encuentra ubicado en el Municipio de Buenavista (Córdoba), y esta compuesto por búfalos mestizos con alto linaje de la raza Murrah, criados en pastoreo, y ordeñados manualmente en la mañana y en la tarde. La suplementación de la hacienda se lleva a cabo en la época de verano (diciembre a marzo), es realizada con heno de pasto Climacuna (Dichanthium annulatum).
Para el cálculo de las producciones hasta los 240 y 270 días en leche (PL240 y PL270), se utilizó el procedimiento recomendado por el Comité Internacional de registro animal (ICAR, 2002), según la expresión [1]
En que:
PL = producción acumulada de leche hasta 240 ó 270
días (PL240 ó PL270 días)
Ci = Producción de leche en el i-ésimo
control, i Î [1, n];
D1 = Días al primer control;
Di = Días en leche, i Î [2, n];
Ds = Dia de secado.
C1 = Producción de leche en el primer
control
Dn = Días al último control
Cn = Producción de leche en el último control
Las lactancias con mas de 240 días se agruparon al azar en 15 clases de días en leche (de 100 a 249 días con intervalos de 10 días) y se recalculó la producción de leche acumulada hasta el día de corte de cada clase (PD240), utilizando la expresión [2].
En que:
Pd =Producción de leche total al corte al azar;
C1 =Producción de leche en el primer
control;
Ci = Producción de leche en el i-ésimo
control, [1-u];
Cu =Producción de leche del último control
antes del corte al azar;
D1 = Intervalo entre el dato del parto y del primer
control;
Di = Dias en leche, i Î [2≤ u];
Dd= Días de leche donde se hizo el corte al
azar;
Du= Días en leche del último control, donde u ≤
d
De igual forma, se recalculó la producción de leche hasta el día de corte (PD270) de 18 clases, para lactancias con más de 270 días en lactancia (Tabla 1). Para la PD240 y PD270 el número de lactancias asignadas en cada clase varió entre 148 y 213 y 74 y 118 para un total de 2548 y 1736 lactancias de 1059 y 959 búfalas, respectivamente. (Tabla 1).
Tabla 1. Distribución de las lactaciones de búfalas en la hacienda Altamar en la Costa Atlántica Colombiana, en clases de duración de la lactación entre 100-240 días y 100-270 |
||||||
240 días en leche |
270 días en leche |
|||||
Clases |
Días en lactancia |
No |
|
Días en lactancia |
No |
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 |
100-109 110-119 120-129 130-139 140-149 150-159 160-169 170-179 180-189 190-199 200-209 210-219 220-229 230-239 240-249
|
179 184 157 161 152 176 171 161 213 148 162 161 172 170 181 |
460±128 508±122 571±147 587±126 642±143 676±152 716±160 758±182 799±165 846±182 899±206 926±186 957±219 1030±214 1034±213 |
100-109 110-119 120-129 130-139 140-149 150-159 160-169 170-179 180-189 190-199 200-209 210-219 220-229 230-239 240-249 250-259 260-269 270-279 |
118 107 74 85 95 77 93 107 92 98 112 86 98 93 85 100 109 107 |
452±104 505±112 544±150 595±148 644±156 734±188 716±176 784±179 796±205 866±190 899±190 945±201 985±182 985±206 1051±241 1107±227 1098±229 1158±242 |
No= número de lactaciones, = promedio y desviación estándar de la producción de leche |
Para la estimativa de las soluciones de cada clase para D240 y D270 se utilizó un modelo animal, con medidas repetidas. Fue utilizado el programa computacional MTDFREML (Boldman et al 1995), se consideró una matriz de parentesco con 5498 animales de ocho generaciones. Las hembras con registros fueron hijas de 149 machos reproductores. Se empleo el siguiente modelo:
yijklm = bi + dj + ak + pl + eijklm,
En donde:
yijklm = variables dependientes (producciones de
leche cortadas al azar).
bi = efecto fijo del i-ésimo grupo
contemporáneo (año, época de parto);
d j = efecto fijo de la j-ésima clase de
duración de la lactación;
ak = efecto aleatorio del k-ésimo animal, con
media 0 y varianza σ2 a ;
pl = efecto aleatorio del medio ambiente permanente, con
media 0 y varianza σ2 c
;
eijklm = efecto aleatorio residual, con media
0 y varianza σ2 e.
Se consideraron las épocas de parto así: verano (de Diciembre a Febrero), época de transición de verano a invierno (Marzo a Mayo), época de transición de verano a invierno (Julio a Septiembre) y época de lluvia (Septiembre a Noviembre).
Las estimativas para las clases de la duración de la lactación fueron obtenidas sumando la media general de cada clase a las soluciones de las clases de duración de las lactaciones, considerando las ecuaciones del modelo anteriormente mencionado. Estas estimativas fueron suavizadas usando un polinomio de segundo grado.
Los factores de corrección fueron estimados mediante la expresión:
Los factores de corrección estimados fueron aplicados a las producciones de leche para diferentes días en leche, obteniendo así, las producciones corregidas. Estas producciones fueron comparadas con las producciones observadas hasta los 240 y 270 días, con el objetivo de determinar la eficiencia de los factores de ajuste estimados mediante un análisis de correlación.
Para estimar los factores de ajuste a edad adulta a los 240 y 270 días en leche fueron utilizadas 3153 lactancias de 1150 búfalas respectivamente, estas producciones fueron distribuidas en clases de edad y época del parto de las búfalas.
Los factores de corrección fueron estimados mediante la utilización del siguiente modelo:
Yijklm= Hi+Ij +ak+pl+eijklm
Donde:
Yijklm= Producción de leche hasta 240 y 270
días.
Hi = Efecto fijo del i-ésimo grupo
contemporáneo (año del parto)
Ij = Efecto fijo de la j-ésima clase de edad -
época de parto.
ak = Efecto aleatorio del k -ésimo animal, con
media 0 e varianza σ2a ;
pl= Efecto aleatorio de ambiente permanente de la
búfala, con media 0 de varianza
σ2c, y
eijklm = Efecto aleatorio residual, con media 0 e
varianza σ2e
Los factores de ajuste fueron calculados mediante la expresión utilizada por Iloeje (1980), definida como:
Fj = Mjb / [Mjb -(Ajb -Aj)],
Donde:
Fj es el factor de ajuste multiplicativo para
la clase de edad - estación de parto j.
Mjb es la
media de producción de leche de la clase considerada como base
búfalas a edad adulta.
Ajb es la solución
obtenida en la clase de duración de la lactación escogida
como base y
Aj es la solución en la clase de edad -
y estación de parto j.
Los factores multiplicativos Fj fueron suavizados dentro de cada estación del año y utilizando el método descrito por Wiggans (1981).
Las producciones de leche ajustadas y no ajustadas a edad adulta fueron comparadas por la reducción de la variación asociada a los efectos de la edad y la estación de parto utilizando análisis de varianza, incluyendo efectos época - edad (grupo época y edad de la búfala), año de parto y búfala con el propósito de determinar la reducción de la variación asociada a la edad del animal (Searle y Henderson 1959; Miller et al 1968; Searle 1962).
Los factores de ajuste estimados para los 240 y 270 días en leche se encuentran en la tabla 2, se observó que los factores de ajuste a los 100 días fueron de 2.404 y de 2.465 respectivamente.
Tabla 2. Factores de corrección multiplicativos para el ajuste de producción de leche a 240 y 270 días en lactancia en búfalas de la costa atlántica colombiana |
||
Días en lactancia |
Factor de ajuste |
|
Para 240 días |
Para 270 días |
|
100-109 110-119 120-129 130-139 140-149 150-159 160-169 170-179 180-189 190-199 200-209 210-219 220-229 230-239 240-249 250-259 260-269 270-279 |
2.404 2.163 1.969 1.809 1.675 1.561 1.463 1.378 1.304 1.238 1.180 1.128 1.081 1.039 1.000 |
2.465 2.237 2.051 1.896 1.766 1.654 1.558 1.474 1.400 1.335 1.277 1.225 1.178 1.135 1.097 1.062 1.029 1.000 |
Las medias de producción de leche hasta y ajustada a los 240 días fueron de 1025 ± 213 y 1063 ± 243 kg, respectivamente, con una correlación entre las dos variables de 0.92. Las producciones de leche hasta y ajustada a los 270 días fueron de 1152 ± 229 y 1143 ± 259 kg, respectivamente, con una correlación de 0.92.
Lo anterior indica que los factores de ajuste fueron eficientes y confiables y pueden ser utilizados para ajustar lactancias incompletas. Sin embargo, la eficiencia del ajuste es menor en hembras con pocos días en leche, con correlaciones entre las producciones de leche hasta y ajustadas de 0.85 para 240 días y 0.79 para 270 días en la clase 1, pero esta eficiencia va mejorando a medida que la hembra tiene mas días en leche, con correlaciones entre las producciones de leche hasta y ajustadas de 0.99 para 240 días y 0.99 para 270 días en la clase 15 y 18, respectivamente.
En las tablas 3 y 4 se presentan los factores de ajuste para edad adulta. En hembras de 20 a 31 meses de edad, el factor de ajuste a edad adulta varió de 1.129 a 1.215 en la producción de leche hasta los 240 días, indicando que hembras de esta edad producen entre el 79 y el 83% de una hembra adulta, dependiendo de la época de parto (Tabla 3).
Tabla 3. Factores de corrección multiplicativos1 para el ajuste a edad adulta de la producción de leche hasta 240 días de búfalas de la costa atlántica colombiana |
|||
Edad de la búfala (meses) |
Factores de ajuste |
||
Diciembre-Febrero |
Julio-Septiembre |
Octubre-Noviembre |
|
20-31 32-38 39-44 45-51 52-59 60-67 68-75 76-85 86-95 96-108 109-125 126-299 |
1.129 1.093 1.065 1.044 1.029 1.020 1.016 1.018 1.025 1.038 1.056 1.082 |
1.103 1.103 1.071 1.047 1.028 1.014 1.005 1.000 1.000 1.004 1.012 1.026 |
1.215 1.169 1.133 1.107 1.088 1.077 1.072 1.074 1.082 1.098 1.121 1.152 |
1Clase Base: búfalas entre 86 y 95 meses de edad y ocurridos en el mes de Julio-Septiembre |
En el caso de la producción de leche hasta 270 días, las hembras de 20 a 31 meses producen entre un 83 y 91% de una hembra adulta (Tabla 4).
Tabla 4. Factores de corrección multiplicativos1 para el ajuste a edad adulta de la producción de leche hasta 270 días de búfalas de la costa atlántica colombiana |
|||
Edad de la búfala (meses) |
Factores de ajuste |
||
Diciembre-Febrero |
Julio-Septiembre |
Octubre-Noviembre |
|
20-31 32-38 39-44 45-51 52-59 60-67 68-75 76-85 86-95 96-108 109-125 126-299 |
1.094 1.066 1.045 1.030 1.019 1.013 1.011 1.013 1.021 1.033 1.049 1.072 |
1.118 1.084 1.057 1.035 1.019 1.008 1.001 0.998 1.000 1.006 1.016 1.031 |
1.167 1.131 1.102 1.081 1.067 1.058 1.056 1.058 1.067 1.082 1.103 1.132 |
1Clase Base: búfalas entre 86 y 95 meses de edad y ocurridos en el mes de Julio-Septiembre |
Estos factores fueron aplicados a las de producciones de leche hasta los 240 días, donde las producciones de leche sin ajustar y ajustada a edad adulta fueron de 954.18±214.93 y 1009.306 ±224.76, respectivamente (Tabla 5).
Tabla 5. Análisis de varianza para lactaciones hasta 240 días ajustadas y no ajustadas a edad adulta de búfalas de la costa atlántica colombiana |
|||||
Fuente de Variación |
GL |
Producción de leche |
Producción de leche ajustada a edad adulta |
||
CM |
Prob F |
CM |
Prob F |
||
Edad de la búfala y época de parto |
35 |
104086 |
0.0036* |
62026 |
0.24n.s. |
Año de parto |
7 |
743188 |
0.0001* |
811478 |
0.0001* |
Búfala |
1062 |
110714 |
0.0001* |
123751 |
0.0001* |
Error |
1713 |
46197 |
|
50517 |
|
Total |
2817 |
|
|
|
|
Coeficiente de determinación |
64.8 |
64.3 |
|||
Promedio de leche (kg) |
954 |
1009 |
|||
Coeficiente de Variación |
22.5 |
22.7 |
|||
GL= grados de libertad, CM=cuadrado medio, n.s.=diferencias estadísticas no significativas y *=diferencias estadísticas altamente significativas (P≤0.01) |
En los análisis de varianza para las producciones de leche sin ajustar se encontró que el efecto edad-época de parto fue altamente significativo (P≤0.01), y en la producción de leche ajusta fue no significativo (P≥0.05), indicando que los factores de ajuste removieron el efecto del grupo de edad de la búfala y época de parto en la producción de leche hasta los 240 días (Tabla 5).
También se aplicaron los factores de ajuste a las producciones de leche hasta los 270 días, donde las producciones de leche sin ajustar y ajustada a edad adulta fueron de 1027.874±237.01 y 1069.59 ±244.86, respectivamente (Tabla 6).
Tabla 6. Análisis de varianza para lactaciones hasta 270 días ajustadas y no ajustadas a edad adulta de búfalas de la costa atlántica colombiana |
|||||
Fuente de Variación |
GL |
Producción de leche |
Producción de leche ajustada a edad adulta |
||
CM |
Prob F |
CM |
Prob F |
||
Edad de la búfala y época de parto |
35 |
99733 |
0.0001* |
69433 |
0.17n.s. |
Año de parto |
7 |
997945 |
0.0001* |
1060714 |
0.0001* |
Búfala |
1062 |
133129 |
0.0001* |
144649 |
0.0001* |
Error |
1713 |
46197 |
|
59957 |
|
Total |
2817 |
|
|
|
|
Coeficiente de determinación |
64.7 |
64.6 |
|||
Promedio de leche (kg) |
1028 |
1070 |
|||
Coeficiente de Variación |
23.1 |
22.9 |
En los análisis de varianza para las producciones de leche sin ajustar se encontró que el efecto edad-época de parto fue altamente significativo (P≤0.01), y en la producción de leche ajusta fue no significativo (P≥0.05), indicando que los factores de ajuste removieron el efecto del grupo de edad de la búfala y época de parto (Tabla 6), en la producción de leche hasta los 270 días.
Las medias de producción de leche fueron bajas, comparadas con las encontradas en Brasil, donde la producción de leche corregida para los 240 y 270 días fueron de 1535 ± 566 y 1712 ± 347, respectivamente (Tonhati et al 2004), las encontradas en Italia, donde la producción de leche corregida para los 270 días fue de 2286± 492 kg (Rosati y Van Vleck 2002), en Pakistán, donde la producción de leche ajustada a 308 días de búfalos Nilli-Ravi fue de 2139 kg (Khan y Chaudry 2001). La mayor producción de leche en Italia, Brasil y otros países europeos se debe probablemente a las diferencias en la alimentación animal, duración de la lactancia y la constitución genética.
Se encontró que al aplicar los factores de ajuste para edad adulta se pudo remover el efecto del grupo de edad de la búfala. La mayor producción de leche fue alcanzada por búfalas que parieron en la estación que comprendía los meses de julio a septiembre, con una edad de 86 a 95 meses (Esta fue considerada la edad adulta). No se tuvo en cuenta la estación de Marzo a Junio en las tablas de factores de corrección para edad de la búfala (Tablas 3 y 4) por no existir datos suficientes para un análisis confiable, debido a que las búfalas tienden a parir en determinadas épocas del año.
No se encontraron investigaciones realizadas en los reportes de la literatura para la comparación y discusión de factores de ajuste a edad adulta en búfalos. Sin embargo, en bovinos se han encontrado edades adultas en ganado Holstein en Colombia entre los 76 a 86 meses de edad (Stanton 1990; Martinez 1993; Cerón-Muñoz et al 2003).
El uso de factores de ajuste permite predecir la producción de leche de las búfalas en forma eficiente.
Se recomienda el uso de factores de ajuste para duración de lactancia, para búfalas que aun no han terminado la lactancia, para predecir su producción total ó hembras que secaron por causas ajenas al desempeño productivo del animal ó salieron de la finca sin secarse.
El uso de factores de ajuste es para la finca analizada. Sin embargo, pueden ser utilizados para otras explotaciones con semejantes condiciones ambientales y de manejo.
Este Proyecto fue financiado por el "Fondo para apoyar trabajos de Grado en los programas de Pregrado" (Vicerrectoría de Investigación de la Universidad de Antioquia). Los autores agradecen a la hacienda Altamar S.A, por facilitar la información para el desarrollo de este trabajo y al centro de investigaciones de la Facultad de Ciencias Agrarias (CIAG) por el apoyo.
Boldman K, Kriese L, Van Vleck L, Van Tassell C and Kachman S 1995 A manual for use of MTDFREML. A set of programs to obtain estimates of variances and covariances. Lincon: Department of Agriculture Research Service. 120 p.
Cerón Muñoz M, Tonhati H, Costa C, Solarte C y Benavides O 2003 Factores de ajuste para producción de leche en bovinos Holstein colombiano. Revista Colombiana de Ciencias Pecuarias. Volume 16, pp. 26-32. http://kogi.udea.edu.co/revista/16/16-1-4.pdf
ICAR 2002 Section 2.4 Guidelines for buffalo milk
recording for low to medium and medium to high input production
systems. In: International agreement of recording practices.
International Commitee in Animal Recording www.icar.org. Accesado el 20 de mayo de
2005.
http://216.239.59.104/search?q=cache:sZwwoIKV6CoJ:www.icar.org/docs/technical_series/tec_series_4_bled_buffalo.pdf+Guidelines+for+buffalo+milk+recording+for+low+to+medium+and+medium+to+high+input+production+systems&hl=fr&client=firefox-a
Iloeje M 1980 Age-adjustment factors for Alpine, LaMancha, Nubian ,Saanen, and Toggenburg dairy goats, Journal of Dairy Science. Volume 63, pp 1309-1316.
Khan M and Chaudry H 2001 Prediction of lactaction yield from last- record day and average daily yield in Nili- Ravi Buffaloes. Bulletin of the FAO Inter-regional Cooperative Research Network on Buffaloes, Volume 15, pp 7-10.
Martinez M 1993 Determinación de factores de ajuste para ganado Holstein en la Sabana de Bogotá. Tesis de Maestría, Facultad de Ciencias Pecuárias, Universidad Nacional de Colombia, Bogotá. 86 p.
Miller R 1973 A recent study of age adjustment. Journal of Dairy Science, Volume 56, pp 952-959.
Miller R, Mc Daniel B and Myowman R 1968 Effects of errors in the edge adjustment of first lactactions. Journal of Dairy Sience, Volume 51, pp 378-384,
Rosati A and Van Vleck L 2002 Estimation of genetic parameters for milk, fat, protein and mozzarella cheese yield in the Italian river buffaloes population. Livestock Production Science, Volume 74, pp185-190.
Schmidt G and Van Vleck L 1974 Bases científicas de la producción lechera. Edición Acribia; España, 583 p.
Searle S 1962 Age and herd effects in New Zealand dairy cow records. Journal of Dairy Science Volume 45, pp 82-85.
Searle S and Henderson C 1959 Establishing age-correction factors related to the level of herd yield. Journal of Dairy Science, Volume 42, pp 824-835.
Stanton T 1990 Investigation of genotype by environment interaction for Holstein milk yield in Colombia, Mexico y Puerto Rico. Tesis de Doctorado, Department of Animal Science, Cornell University., Ithaca, 623p.
Tonhati H, Muñoz-Cerón M, Duarte J, Reichert R, Oliveira J and Lima a 2004 Estimates of correction factors for lactation length and genetic parameters for milk yield in buffaloes. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, Volume 56, pp 251-257 http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0102-09352004000200016&lng=en&nrm=iso
Warwick E and Legates J 1980 Cria y mejora del ganado. 3ra edición México: McGraw_Hill,. 623 p.
Wiggans G 1981 Smoothed Age- Season Adjustments Factors for Dairy Goat Lactation Milk and fat records. Journal of Dairy Science, Volume 64, pp 350-352.
Received 29 June 2005; Accepted 26 January 2006; Published 11 April 2006