1 PB
1 PB
1 PB
RESUMEN
Rev. Fac. Nal. Agr Medellín. Vol 48 Nos. 1 y 2, p.41 -68 . 1995 41
Arboleda, E.; Valencia, C. y Montoya, C. Comparación entre factores de ajuste ..
fue ron llevados los demás registros; luego se obtuvieron factores aditivos estimated constants by Harve's program and multiplicatives frol
mediante el uso de constantes estimadas por el programa Harvey y mean olleast squares.
multiplicativos a partir de la media de mínimos cuadrados. Una vez
ajustados los datos se corrió nuevamente el modelo inicial, con el fin de When data were adjusted it ran to initial model to quantify thé
cuantificar la varianza residual presente tanto para factores de ajuste dual variance which is present for multiplicative and additive adjus i
multiplicativos como aditivos, a los dos niveles propuestos. La compara facto rs for the two proposal le vels. Observing ratios of table F anú
ción se realizó mirando las relaciones de la tabla F y las varianzas dual varia nces between g roups and within them was mad
residuales entre grupos y dentro de ellos. companson.
Aunque no se presentaron diferencias marcadas, se observa que, en A lthoug h, significant di fferences were not present it was know tf
general los factores de ajuste aditivos tuvieron un mejor comportamiento general the additives adjustm ent lactors had a best behavior tha
que los multiplicativos, pues con el empleo de estos últimos queda una multiplicative ones b cause 01 a little variati n mQJJ..ed ~ ain
pequeña variación sin rem over; en tanto que con factores de ajuste aditi the lro'
vos la corrección se hace m ás precisa y no se altera la varianza dentro de bf
grupos ajustados. Así mismo para factores de ajuste multiplicativos, se
presenta una pequeña ventaja cuando se utiliza como base el nivel de 1
máxima producción, frente a la utilización de la media de mínimos cua I
drados estimada por H arvey.
INTRODUCCIO
42 Rev. Fac. Nal. Agr. Medellín. Vol 48 Nos. 1 y 2, p.41-68. 1995 Re". Fac. Nal. Agr.
Arboleda, E.; Valencia, C. y Montoya, C. Comparación entre factores de ajuste...
fueron llevados los demás registros; luego se obtuvieron factores aditivos estimated constants by Harve's program and multiplicatives from the
mediante el uso de constantes estimadas por el programa Harvey y mean of least squares.
multiplicativos a partir de la media de mínimos cuadrados. Una vez
ajustados los datos se corrió nuevamente el modelo inicial, con el fin de When data were adjusted it ran to initial model to quantify the resi
cuantificar la varianza residual presente tanto para factores de ajuste dual variance which is present for multiplicative and additive adjustment
multiplicativos como aditivos, a los dos niveles propuestos. La compara facto rs for the two proposallevels. Observing ratios oftable F and resi
ción se realizó mirando las relaciones de la tabla F y las varianzas dual variances between g roups and within them was made the
residuales entre grupos y dentro de ellos. companson.
Aunque no se presentaron diferencias marcadas, se observa que, en Although, significant differences were not present it was know that in
general los factores de aju te aditivos tuvieron un mejor comportamiento general the additi ves adjustm ent factors had a best behavior that the
que los multiplicativcw. lJues { 'n n Al ,,--~ - .~ "'-stos últimos queda una multipli ative ones because of a liule variation unremoved remain from
tore s de ajuste aditi the latest, in the meantime with the additive adjustment factor s correction
varianza dentro de became m ost accurate and the variarLce within adjusted groups was not
multiplicativos, se alt red. T hus for multiplicative adjustment factors a líttle advantage is
o base el nivel de presented when the maximum prodLL ction level is used as basis to face
ia de mínimos cua the use ol mean least quares e timated by Harvey p rogramo
INTRODUCCION
Los factores de ajuste fu r on creados con el fin de minimizar el elTor m elio
ambiental en la comparación de anim a les con difer ncias m ar adas n cuan to a edad,
T FACT ORS raza, año de parto, época de parto y o fOSo
YIELDS IN A
Se ha encontrado que n 1 caso de las vacaf'i Holstein los fac tores más influyentes
sobre la producción por lactan ia son la duración de la misma, la edad, el año de
parto, los días abí rtos, la época y el número d pa rto; sin emba rgo es im portante
erd "P aysandú"
realizar estudios en los cuales se pr ci e cuál facto r de aj uste produce un a mayor
Elena sidewalk,
seguridad en la reducción de estos efecto::;. La presente inve ti ación pretende a
ch.
través de la utilización de modelos fijos , determinar cuáles factores (aditivo y/o
multiplicativos ) reduc n más la variación medio-ambicntal.
tive advantages
s to produce by La estimación de las diferencias en el uso de factores de aj uste multiplicativos o
due to lactation aditivos permitirá decidir, al momento de hacer estimacion es gen ' ticas, cuál de los
by way of least dos produce una mayor preci Ión para la r educción de la variación d bido a otras
l environmental causas. omo consecuencia de lo antelior podrá hacerse una scl ección m ás adecuada
~ year; prevwus
d los individuos y un m j o!" uso de los factores de aj uste de acuerdo co n las
nd calvin time circunsLancia de variación propias del hato n studio.
ltiplicative and
subclass was El obj t ivo del presenL ' t rabajo fué del rm inar las v n ajas comp r ativas en tre
ds were carried los factores de ajuste multiplicativos y aditivo' para algunos efectos no gen 'ticos
btained using sobre la producción de leche.
1 Y2, p.41 -68. 1995 Rev. Fac. Na!. Agr. Medellin. Vol 48 Nos. 1 y 2, p.4 1-68. 1995 43
Arboleda, E .; Valencia , C. y Montoya , C. Comparación entre factores de ajuste ...
REVISION DE LITERATURA definida de producción al propio animal y en los multiplicativos lo que se aumeI
una parte proporcional de la producción del animal.
Los factores de ajuste , como se anotó previam ente, tienen como finalidad primor
dial corregir el error medio-ambiental ocasiona do por especie, raza, manejo-hato, Según Ochoa y Restrepo (1986), para obtener factores de corrección es nece
año, entre otros efectos que causan la sobre-estim ación del valor genético de algunos estimar las producciones promedio para cada grupo establecido; posteriormer
animales y la subestimación de otros . hace una relación entre los promedios de producción de los diferentes grupos, (
grupo al cual se quieren llevar todos los registros.
Falconer (1970 ) ilustra claram ente cómo el fenotipo es una expresión del genotipo
y del medio-ambiente; ello indica que a medida que se minimice el efecto de és te Otros autores han discutido sobre el uso de factores multir>licati O!;J y/o aditiva
último, la valoración fenotípica tendrá una m ayor aproximación a la varianza genética
de los individuos y como consecuencia los animales que presenten mejor calidad lo Searle y H
serán en virtud de la heren cia. hato x edad
Además, di
Cundiff; Willham y Pratt (1966 proponen por su parte, para defi nir el tipo de del hato, pru
factor a usar, comparar las m edias y las varianzas entre las subclases.
Miller; L.
También se h a dedicado u n amplio espacio al análisis del tipo de factor de aj uste estacionales
44 Rev. Fac. Na!. Agr. Medellín. Vol 48 Nos. 1 y 2. p41-68. 1995 Rev. Fac. Nal. Agr. Me
Arboleda, E.; Valencia, C. y Montoya, C. Comparación entre factores de ajuste...
Martínez; Lee y Linn (1990) recomiendan por su lado el uso de factores de ajuste
multiplicativos cuando se trata de corregir para la interacción edad-estación en
registros Holstein lecheros; ellos ajustan simultáneamente para ambos efectos.
.1995 Rev. Fac. Nal. Agr. Medellín. Vol 48 Nos. 1 y 2, p.41-68. 1995 45
Arboleda, E.; Valencia, C. y Montoya, C. Comparación e ntre factores de ajus te...
sus datos. Aunque los métodos aditivos son satisfactorios y por consiguiente no deben ner factores de ajuste; la gran mayoda de ellas, sin embargo, son similares entre
descartarse. pueden considerarse como adaptaciones de las tres principales enunciadas
Schmidt y Van Vleck (1976), las cuales son :
1. comparación grosera: ésta utiliza el promedio simple para todas las vaca~
INTERACCIONES cada edad. Se agrupan las vacas de primera lactancia, formando otro grl
con las de segunda lactancia y así sucesivamente; luego se procede a compa
El problema de las interacciones es que ellas pueden producir sesgos en los esti el promedio simple de cada grupo de lact . con resDAc:t.o ~~.. _l
mados de parámetros genéticos, haciendo qu e muchas veces se enmascaren algunos cido como bp
efectos y que estos u otros no sean considerados. cas rllás "
estas va
El tipo de interacción más a mpliamente debatida en la literatu r es de la edad equival.
con otros efectos medio-ambientales. Un vist zo a los probl em as básicos relativos a vadas.
los ajustes por edad en los registros de pr o ucción es presentado por Freeman (1973);
Hickman (1973 ); Miller (197 3); McDaniel (1973); el! s llaman la aten ción sobre la 2. compar~
necesidad de ajustar únicamente pa r a los efectos de la edad bajo las mismas condi misma v
ciones medio-ambientales en m adurez, sin elimin ar los fectos de selección y otros todas la!
factores. anterior
n o es del
valente e
Everet t, Taylor y Hammon d (1982 ) h allan igualmente dependen cia e n tre m es y que tien
luego de ajustar los registros por edad y mes de parto; a unque í la h a bí a hallado factores
46 Rev. Fac. Nal Agr. Medellín. Vol 48 Nos. 1 y 2, p.41-68. 1995 Rev. Fac Nal Agr. Medel
Arboleda, E.; Valencia, C. y Montoya, C. Comparación entre factores de aj us te ..
s datos. Aunque los métodos aditivos son satisfactorios y por consiguiente no deben ner factores de ajuste; la gran mayoría de ellas, sin embargo, son similares entre sí y
escartarse. pueden considerarse como adaptaciones de las tres principales enunciadas por
Schmidt y Van Vleck (1976), las cuales son:
l. comparación grosera: ésta utiliza el promedio simple para todas las vacas de
JTERACCIONES cada edad. Se agrupan las vacas de primera lactancia, formando otro grupo
con las de segunda lactancia y así sucesivamente; luego se procede a comparar
1<'1 - sesgos en los esti el promedio simple de cada grupo de lactancia con respecto a un grupo estable
mascaren algunos cido como base. Este método está sometido a desviaciones puesto que las va
cas más viejas serán generalmente un grupo selecto; por tanto el promedio de
estas vacas presenta desviaciones hacia producciones altas y los factores para
ura es de la edad equivalente maduro de vaca jóvenes mostrarán desviaciones hacia cifras ele
ásicos rela tivos a vadas.
Freeman (1973);
ten ción sobre la 2. comparación duplicada : se basa en la comparación de las producciones de la
s mismas condi misma vaca a edades diferentes; luego se promedian estas comparaciones para
elección y otros todas las vacas. Esto puede provocar una desviación en la dirección opuesta al
anterior método, puesto que la repetibilidad en el rendimiento de la lactación
no es del cien por ciento, ya que solamente cabe esperar que se repita la mitad
aproximadamente de la superiori.dad de una vaca. Al calcular por lo bajo el
incremento en la producción con el aumento de edad, los factores para el equi
valente adulto pueden ser muy pequeños para vacas jóvenes.
970), rev lan
o al h acer las Algunos investigadores han sugerido promediar los factores obtenidos por
o de evitar al ambos métodos, a fin de obtener factores de ajuste más exactos,
modelo fijo, se estaría presumiendo la completa repetibilidad para todas las observa MATERIALES YMETODOS
ciones,
Materiales
La discusión en torno a cuál o cuáles son las metodologías más apropiadas para la
obtención de factores de ajuste adecuados, sigue planteada; pero hoy en dia como Localización
puede observarse en la Tabla 1, existe la tendencia a emplear modelos mi,,'tos y el
Los da
método de semejanza m áxi ma al estimar parámetros no genéticos que afectan las
exp lota i,
producciones lácteas en las diferentes regiones del m undo, sin que su uso tan gene
versida
ralizado indique que sea el más preciso p ara tod s las situaciones,
corregimi.
19 5, a Ul
de 2600m
2500 m .mI
Tabla l. Algunos de los trabajos m ás n otables sobre fact ores de ajuste en ganado
lechero.
Descripción \
I
Autor Año Re~ión No, Registros- Model 1étodos Tipo de ajus te Efectos Se utí .
R.:tzil (1 usad os nxomenda do significativos
Chauhan 1988 A1pe;<; 168.4S0S-HR Mixto I\I L Ad itivo EdJ d ·Mes clasi
lJastu r.l-H.:lto
Martínez el al 1990 l3rasil 2488 1-HF Mixto ML-PC-CC Multiplicativo Edad- Es taci ón 1:$
2:33
Agudl'lo- 1991 Colombia 535-HF Fijo ,C Edad-Año-Días
Cardona ahi<'rtas-Epoca-
3: 46j
Número de 4: 59J
parto 5: 72 1
Djemali- 1992 Tu ncz 2,5,OOO-I-lF Fij o CC Resión-Hato- 6: 85':
Berger Mes-Año 7: ;:: ,
HF: Holstein Fr lSlan. J: .Jersey. G : Guernsey. A: Ayrshire. AMZ: Cebú lecher año d
AIZ: Criollo. S: Simen tal. HR: Holst in Rojo . ML: Sem janza máxima. PC: Compa
ración duplicada. GC : Compara ció n gros era. previd
1: $6j!
2: 61
3: 101
4:141
5:> 18
48 Rev, Fac, Na!. Agr. Medellín, Vol 48 Nos , 1 y 2, p.41-68, 1995 Rev, Fac, Na!. Agr.
Arboleda, E.; Valencia , C. y Montoya, C. Comparación entre factores de ajuste.. .
MATERIALES Y METODOS
odelo fijo, se estaría presumiendo la completa repetibilidad para todas las observa
ones. Materiales
La discusión en torno a cuál o cuáles son las metodologías más apropiadas para la Localización
tención de factores de ajuste adecuados, sigue planteada; pero hoy en dia como
ede observarse en la Tabla 1, existe la tenrl- - . ear modelos mixtos y el Los datos p ara la r ealización del presente estudio fueron tomados del centro de
étodo de semejanza máximR 0 1 - .' éticos que afectan las explotación lechera P aysan dú de la Facu ltad de Ciencias Agropecuaria. de la Uni
roducciones láct.p·~- que su uso tan gene- versidad NacionaJ de Colombia, Sede Medellín. Dicho Cen tro está ubi a do en el
corregim ien to Santa E len a del Municipio de Medellín; correspondiente, según Espina l,
1985, a un a zona ecológica de bosque húm edo Montano Bajo (bh-MB), con una altu ra
de 2600m. ·.n .m., temp ratura promedio de 12.fLC. y una precipitación promedia de
2500 m .m. anuales .
Descripción de materiales
Efectos Se utilizaron 55 regio tros d lactancia proveni ntes de 195 vaca s Holstein que
significativos tuvieron sus partos entre los años 197(;;-1 89 e n I h ato P ysandú .
Edad
P r lactancia se tomó la siguicnt información:
número de la vaca,
EJad-Mes
edad en cada uno de los partos. Esta se agrupó de 1-7, obedeeien do a la siguiente
clasificación:
1: ~32 ro se.
2: 33-45 m ses
3: 46-58 m ses
4: 59-71 m 3SCS
5: 72-84 meses
6: 85-97 m se
7: ~ 98 meses
a ñ o de parto ,
: ~60 dí a~
2: 61-100 días
3: 10 1-140 días
4: 14 1-180 días
5: >180 días
Para clasificar las novillas de primer parto dentro de las subclases de días abiertos Tabl a 2. gr pación de las ob' l-vaciones dentro de las difer ntes clases y subcl
previos, se tuvo en cuenta el período transcurrido desde el prim r calor s rvido hasta
el primer servicio fértil.
1: Diciembre a Marzo,
2: Abril y Mayo,
3: Junio a Septiembre,
4: Octubre y Noviem bre.
La lact ncia real fué corregida a 305 días empleando factor s de corrección
multiplicativos.
50 Rev. Fac. Nal. Agr. Medellín . Vol 48 Nos. 1 y 2, p.4 1·68. 1995 Rev. Fac.
Arboleda , E. ; Valencia, C. y Montoya, C. Comparación entre factores de ajuste..
Para clasificar las novillas de primer parto dentro de las subclases de días abiertos Tabla 2. Agrupación d las observaciones dentro de las diferentes clases y subclases.
evios, se tuvo en cuenta el período transcurrido desde el prim r calor servido hasta
primer servicio fértil.
SES SUB LASES No. DE
número de partos, OBSERVAOON -5
2 1::;4
,..',, 118
4 88
uerdo a la gr fica
5 52
7 20
1976 25
1977 29
1978 36
1979 41
1980 39
1981 50
1982 41
de correcci ón 19133 47
1984 51
1985 52
1986 46
1987 39
1988 27
1989 19
8. 1995
Rev, fac. Nal. Agr. Medellín, Vol 48 Nos. 1 y 2, p.41-68. 1995
51
Arboleda, E .; Va lencia, C. y Montoya, C. Comparación entre factores de ajuste ..
Continuación Tabla 2, Ji = media de mínimos cuadra dos, común a toda s las observaciones
A¡ = efecto de la edad de parto, i varia de 1-7,
CLASES SUBCLASES No. DE Bj efecto de la época de parto, j varía de 1-4,
OBSERVACIONES Ck efecto del año de parto, k varía de 1-15,
Número de parto 1 179
DI = número de parto, 1 vati de 1-8,
Fm previos días abiertos, m varía de 1-5,
2 136 b y/x efecto de la regre sión de la duración de 1 producción de lechE
3 102
sobre la duración de la lactancia (Y),
4 73
5 45
E¡jklmn = error experimental común a to das las observaciones, N(O, S2).
6 15
7 8 Par
Días abiertos 1 62 Y U'
2 174
3 123
Do
4 82
5 117
Y ijl
Epoca de parto 1 167
2 82 Lo'
3 203
Pa
4 106
me
Yij'
Y ¡j
METODOS
Y ij
Donde : La
descr;
Yijklmn = registro n de lactanci a r a l, a la edad de parto i, en la época j del año
k, con número de parto 1 y m días abiertos previos,
52 Rev. Fac. Nal. Agr Medellín. Vol ~8 Nos. 1 y 2, p . ~1-68. 1995 Rev. FaQ
Arboleda, E.; Valencia, C. y Montoya, C. Comparación entre factores de ajuste ..
·nuación Tabla 2. fA. media de mínimos cuadrados , com ún a todas las observaciones ,
CLASES
OB, ERV ACIONES Ck efecto del año de parto, k varía de 1-15,
de parto 1 179
Fm previos días abiertos, m varía de 1-5,
136
by/x efecto de la regresión d la du ración de la producción de leche (x),
Donde:
modelos:
Donde:
AB lJ· = A¡ + Bj + (AB)¡f efecto común t odos los registros en la ed d i, época j,
donde (AB)¡j r presenta la interacción edad j x época j,
Los demás efectos, como fueron denomin ados anteriormente . Todos los efectos
descritos en los modelos fueron considerados fijo~ .
a ño
. 1995 Rev. Fac. NaL Agr Medellín. Vol 48 Nos. 1 y 2, pA1·68 1995 53
Arboleda , E.; Valencia, C, y Montoya, C, Compara ción entre factores de ajuste ...
De acuerdo a la metodología empleada por Miller et al (1966); Chauhan (1988); Se ajustaron los datos por los factores multiplicativos y aditivos obteni<
Martínez Lee y Lin (1990), los factores de ajuste fueron obtenidos de la siguiente posteriormente se corrió el modelo anteriormente enunciado para lactancia corre€
forma : a 305 días teniendo en cup"·- 1- ri~G.-.. ..ont.,,,.,f. I'tn" P:;¡r:'l netflclar la variaci
removida en los regist ..
compararon las relaci
Factores aditivos aditiva y multiplicati'
otro, se considera qU\
Se tomó la constante de ajuste de los mínimos cuadrados de cada efecto meno la
constante de ajuste con pon diente la subclase que fué tomada como base a la q ue Posteriormente s~
serían llevados los regis tros. Esta se escogió inicialmente de acuerdo al n ivel de ajustadas por los fac1
producción m áxima. Luego se obtuvieron otros factores usando como base la media L a reducción en 1
de mínimos cuadrados d da por el programa Harvey. Para la subclase t omada como obs rva da para los r
base el factor aditivo fu é cero. significancia estadís
Factores multiplicativos
I
Se elaboraron basados en u na r lación entre la media de los m ínimos cuadrados RESULTADOS Y DISCUS¡
tomada como ba se y la media de los mínim os cua drados para los efectos a ajustar.
Factores no genéticos qu!
Donde:
Para determinar 1
Media de mínimos cuadrados base MMC b ) en cuenta en el pres
FMC¡=
Media de mínimos cuadrados del factor a ajustarClVIMC i )
Inicialmente se anal
(LR) como variable
(Tabla 3). Adem á
1
E l s ubíndice i indica la ubelas del efecto memoambienLal. coeficient de deter
lactancia ajus tad
Para la escogencia de la ba 'C se procedió n las dos formas d scnLa para los es uno de los cfectos~
factores aditivos; para dicha base el factor multiplica tivo fué 1.0 . Lo anterior concue
quien ,5 rec m iend
su trabaj o encon trai
Criterios de comparación I
omo pu ede ob
B sados en lo que expresa hauhan (1988), los criterios para juzgar la efectividad variación de la pro
de los factore de corrección, tant multiplicativos como aditivos, fueron los siguientes: días abiertos (DA), 1J
significati vo (P:-:; O.Oj
1. efectividad par remover la van ción,
también significan9
2. r ducción en la int racción de algunos efectos ,
5~ Re v. Fac. Nal. Agr. Medellín. Vol 48 Nos, 1 y 2. p.41-68. 1995 Rev. Fac . Nal. Agr. Medellinl
Arboleda, E .; Valencia, C. y Montaya, C . Comparación entre fa ctores de ajuste ...
t
tención de los factores de ajuste El procedimiento fué :
De acuerdo a la metodología empleada por Miller et al (1966); Chauhan (1988); Se ajustaron los datos por los factores multiplicativos y aditivos obtenidos ;
artínez Lee y Lin (1990), los factores de ajuste fueron- nht.enidos de la siguiente posteriormente se corrió el modelo anteriormente enunciado para lactancia corregida
a 305 días teniendo en cuenta los diferentes efectos. Para detectar la variación no
orma:
removida en los registros corregidos, se utilizaron las relaciones de la tabla F. Se
compararon las relaciones de F para los a nálisis de varianza de los registros ajustados
aditiva y multiplicativamente; cu a ndo las F son mayores para un método que para el
otro, se considera que es te es menos efi ciente para la corrección.
efecto menos la
o base a la que Posteriormente se compararon la s desviaciones es tándar para las distintas clases
r do al n ivel de ajustadas por los factores tanto m u ltiplicativos com o adit ivos.
o base la media L a r educción en la inter acción de los efectos d scrit os en el modelo no pudo ser
e tomada como observa da para los registros aj ustados, pu sto que estas interaccion es no presentaron
significancia estadística (P::s; O.05).
Lo anterior coin cide con los reportes de diferentes a u Lorcs tales como Schmidt y
rregidas Van VJeck (1976); Costa t al (1984 ); Ev ret t , Taylor y Hammond (1982 ; Wilmink
(1987); Chauh an (1988); Martínez, Lee y Lin (1990) y Dj emali y Berger (1992), quiene
en t rabajos similares a és t e h an reportado signifi L:a n cia, fundamentalmen te pa ra la s
mis mas causas de vaJi ación.
1-68.1995 Rev. Fac. Nal. Agr. Medellin. Vol 48 Nos. 1 y 2, p.4 1-68. 1995 55
Arboleda, E .; Valencia, C. y Montoya, C. Comparación entre factores de ajuste ...
Tabla 3. Análisis de varianza para lactancia real (1). Para la edad al parto, en el presente estudio, el efecto está dado básicamente
la alta producción de la edad 7 (6651.5Kg) frente a las demás.
FACTORES G .L. S.e. CM. F La significancia del efecto días abiertos (Anexo, Tabla 3) está representada pe
baja producción en la subclase 1 (513 2.8Kg) Y la alta producción para la subcla
DE 6 8265376.35 1377562.72 2.902 **
(5778.4Kg), no presentándose diferencia entre las subclases restantes (2,3 y 4), I
APA 14 42680018.25 3048572.73 6.423 ,..,..
sí de estas frente a la un o y a la cinco. Para el efecto época de parto se halló que
NP 5 26845816.79 5369163.36 11.312 ,..,..
subclases y 3 son dife rentes entre sí, siendo las de menor (5322.0Kg) y mayor]
DA 4 2961737.53 740434.38 1.560 N.s.
ducción (5555.2 K g · ~ I - ~ 1,.. Á_~_~ ~~ 1\4__\ _
Se analizan
tos, con los sigl
R2 = 72.2% X = 5586. 99kg CY= 1267.04 CV = 12.33%
P udiera ser
forma distint a '
plo produzcan j
Martínez Lee :
Tabla 4. Aná lisis de ari anza p ar a lactancia corregida a 305 días (Modelo 2). vaca s m adura~
diversa la pro~
seqtúa o de lJuv
I
FACTORES G.L. S.C C M. F se plantearon ~
que ninguna d
DE 6 7119058.10 1186509.68 2.581 ,..
h ato. Com co
APA 14 3 45~
7462.39 2466961.60 5.366 ,..,..
algunas intera
NP 5 25542405.78 5108481.16 11.112 ,..,..
, I
FACTORES DE A
ERROR 525 241364401.91 459741.72
Uso de factores a
R
1
= 49.7% X = 5304.53kg (J' = 927.84 v= 12.78% lnicialmen
par a obten r 1
los aju s tes por I
N.S.: No significativo modelo (2 ) pla
* : Significativo (P:-:;0.05)
56 Rev. Fac. Na!. Agr. Medellín. Vol 48 Nos. 1 y 2, p41-68 1995 Rev. Fac. Nal. Agr. M
Arboleda, E.; Valencia, C. y Montoya, C. Comparación entre factores de ajuste..
bla 3. Análisis de varianza p ara lacta n cia real (1) . Para la edad al parto, en el presente estudio, el efecto está dado básicamente por
la alta producción de la edad 7 (6651.5Kg) frente a las demás.
G.L. S.c. C.M. F La significancia del efecto días abiertos (Anexo, Tabla 3) está representada por la
baja producción en la subclase 1 (5132.8Kg) Y la alta producción para la subclase 5
6 8265376.35 2.902 ** (5778.4Kg), no presentándose diferencia entre las subclases restantes (2,3 y 4), pero
1A 6.423 ** sí de estas frente a la uno y a la cinco. Para el efecto época de parto se halló que las
11.312 ,..,.. subclases 1 y 3 son diferen tes entre sí, siendo las de menor (532 2.0Kg) y m ayor pro
1.560 N.s. ducción (5555. 2Kg) respectivamente ; es de aclarar que la época 1 (Die-Mar) pre se nta
3.712 * la menor precipita ción del año, en tanto que la 3 (Jun-Sept) es una época de lluvias
332.474 ** moderadas; las s u bclases 2 y 4, épocas de alta preci pitación, no ofrecen difer encias
frente a las dem ás .
Se anali za ron además a lgun a s posibles interacciones en los efectos antes descli
tos, con los siguien tes propósitos:
%
P udier a ser que las vacas de determinado grupo de edad fueran afectadas en
forma distinta por los diferentes años y épocas , es decir que vacas jóvenes por ejem
plo produzcan m ás leche en u n os a ñ os y las m ás vi jas en otros a ños o como reporta
Martínez, Lec y Lin (1 990 ), que las novillas respondan mej r a alguna s épocas y las
lo 2). vaca s m aduras a ot r as . P o,' último se pen só q ue las épocas podrían afectar n forma
diversa l a producción den tr'o de los dis tint os a ños analizados, o sea que las ép ocas de
seqlÚa o de lluvias sean m á s critica para unos añ os que para otros . Visto lo a nterior
F se plantearon las inter ac ione edad x a ño, edad x época y época x año, ver ificán dose
que ningun a de ellas p resenta incidencia obr la producción por lactan cia en el
2.581 ,.. h a to. Como consecuencia de estos resulta dos , no pudo s r estudiada la reducción ele
).366 ** algunas interacciones con los factores m ultip li cativos y/o aditivos.
.112 **
.330 ,..,..
.624 ,..
FACTORES DE AJUSTE
Inicia lmente se escogió el nivel de máxi m a pr oducción de cada e fecto como base
p ara obte ner los factores de ajuste de cada u n a de la s subclases. Luego de r alizar
los ajustes por los efe ctos hallados significativos, se procedió a correr nuevamente el
modelo (2) plantea do , L os resultados aparecen en la Tabla 5.
P uede observa rse que el R 2 para el modelo es < 0.01 y la relación de F p ara todos
los efectos analizados fué <0.01. Ello indie que la corrección por todos y cada uno de
los factores medio-am bientales fué hecha acertadamente y una - ez elabora do . ~te
nantes ajuste dichos factores ya no tienen ninguna influencia sobre la produccion; así, esta
rastes podrá deberse más a los efectos genéticos y las vacas que presenten superioridad
serán mejores en virtud de la herencia (Falconer, 1970).
Al observar la media de pr oducción en el modelo ajustado, (Tabla 5) ésta aparece TABLA 6. Aná lís ' de varianza para datos ajustados aditivamente a la prodt:
un poco elevada por cuanto se ajustó a los niveles de m áxima producción. E sto no es media.
inconveniente puesto que la corrección se hace con fines c mp arativos y para factor es
de corrección aditivos, los increm entos decrecimientos en producción de lech dados
por los efectos medio-ambientales, s as umen como los mismos para todas 1 s vacas FACTORES L S.e. c. M.
(Chauhan, 1988). Además a la s aca , on estos factor es, lo que se 1 s agTega es una
cantidad definida de produ cción (Dalton, 1980). Valga decir que este seria el promedio ED 6 11.82 1.97 o.e
de producción si los animales tuvieran sus p artos en la m ejor época del mejor año, APA 14 36.49 2.61 o.e
con un número de día s ab iertos previos óptimo y a una edad madura. NP 5 .97 1 70
""
DA
EP
Tabla 5. Análisis de varian za para datos ajustados aditiva mente al n ivel de máxima ERROR
producción. TOTAL
R2 ::::
FACTORES GL S.e. e. M . F
\
Con el propósito de comp a ra r los factores de ajuste aditivos obtenidos en esta Tabla 7. Análisis ~
investigación para el nivel de m áxima producción y las constantes de ajuste propuestas máxim a producció
por el programa de mínimos cuadrados de Harvey (p) (1988), se corrigieron los registros
por estas y se volvió a hacer un análisis de vrui nza par a el modelo (Tabla 6).
EFECTOS
En éste puede verse como al correr nuevamente el modelo, el R 2 y los valores de F
también dan 0.00 en tanto que la media (5466.3Kg) es similar a la r eal (5586.9Kg) y ED
la desviación estándar continúa siendo 658.27 lo cual indica que, al emplear factores APA
de ajuste aditivos, la varianza es independiente del valor usado como base para corregir NP
los regi stros y que, al adicionar una constante positiva o negativa a los valores no se DA
altera la varianza dent ro de un grupo ajustado (Cundiff; Willham y Pratt, 1966). EP
ERROR
TOTAL
58 Rev. Fac. Nal. Agr. Medellín. Vol 48 Nos . 1 y 2, p.41·68. 1995 Rev. Fac. Na! Ag r. Medellín .
Arboleda, E.; Valencia, C. y Montoya, C. Comparación entre factores de ajuste ...
Al observar la media de producción en el mode lo ajusta do, (Tabla 5) ésta aparece TABLA 6. Análisis de va rianza para datos ajustados aditivamente a la producción
n poco elevada por cuanto se aj ustó a los niveles de máxima producción. Esto no es media.
conveniente puesto que la corrección se h ace con fin es comparativos y para factores
e con-ección aditivos, los increm ento o decrecimient os en producción de leche dados
r los efectos medio-ambientales, se a sumen como los mismos p ara todas las vacas FACTORES GL S. C.M. F
hauhan, 1988). Además a las vacas , con estos factores lo que se les agrega es una
ntidad definida de producción (Daltcm OO{\\ · e este sería el promedio ED 6 11.82 1.97 0.000
e p dur."i' . 1_ época del mejor año, APA 14 36.49 2.61 0.000
fidura. NP 5 8.97 1.79 0.000
DA 4 41.09 10.27 0.000
EP 3 40 .36 1 .45 0.000
al nivel de máxima ERROR 525 241364401.91 459741.72
TOTAL 558 241364597.21
R! = <0.01% X = 5466.27 (J = 658.27 CV = 12.40%
F
tenidos en esta Tabla 7. Análisis de van oza para datos ajustado m ultiplicativamente al nivel de
uste propuestas máxima pr oducción.
on los registros
bla 6).
EFECTO Gl S.c. C.M. F
DS valores de F
I (5586.9Kg) Y ED 6 369230.1.4 61538 .36 0.052
plear factores APA 14 1095715.83 78265.42 0.067
para corregrr P 5 279361.65 55 872 .~ 3 0.047
valores no se DA 4 377319.06 94329 .77 0.080
att,1966). EP 3 3 5767.84 11 9.28 0.101
ERROR 525 617 7488 .22 1176333.11
TOTA L 558 621276931.45
l.
p.41-68. 1995 Rev. Fac. Nal. Agr Medellín. Vol 48 Nos. 1 y 2, p41-68 1995 59
Arboleda, E.; Valencia, C. y Montoya, C. Comparación entre factores de ajuste ..
Es de anotar también que este es el modelo ajustado que mayor desviación estándar factores de ajuste multiplicativos queda una pequeña variación sin remover desr
presenta (1056.12), esto ratifica las afirmaciones de Searle y Henderson, 1960 y de hacer la corrección con estos, en tanto que para los factores aditivos dicha rela
Cundiff, Willham y Pratt (1966) quienes expresan cómo al ajustar multiplicativamente es nula.
se amplía la varianza dentro de grupos.
Al comparar los factores multiplicativos obtenidos teniendo como base el niVE
Igual que para el caso anterior, se obtuvi r n fact ores de ajuste multiplicativos máxima producción y los obtenidos al tomar como base la media de mínimos cuadré
entre la media de m ínim os cuadrados y la m edí de la subclase a ajusta r. Una vez estimada, se pu e de observar cómo estos ú ltimos son más ineficientes, pues
ajustados los datos se h izo el análisis de varian z que aparece en la Tabla 8. El R 2 = relaciones de F son mucho m ayores p ara casi todos los efectos ( Tabla 9), lo (
1.1% queda n do un pequeño efecto residua1 e spe ialmente para el factor m edio indica una p ersist n cia en la variación que n o es removida por el factor de ajust
am biental días a biertos previos.SIendo riguro -os en la comparación de e te análisis
con el modelo inmediatamente anterior (Tabla 7) puede comprobarse cómo en éste
caso la r educción de la variación no fu é tan completa como en el caso anterior pu s en
I
este exi ten r elacione de F Jige>ram nte superiores. Ta bla 9. COn)
para factore~
R2 1. / % X = 5429.31 Kg = =
Tabla 10. Me~
(J 676.13 CV 12.76%
sin aj ust ar y a~
~
Es de anotar también que este es el modelo ajustado que mayor desviación estándar factores de ajuste multiplicativos queda una pequeña variación sin remover después
resenta (1056.12), esto ratifica las afirmaciones de Searle y Henderson, 1960 y de hacer la corrección con estos, en tanto que para los factores aditivos dicha relación
undiff, Willham y Pratt (1966) quienes expresan cómo al ajustar multiplicativamente es nula.
fe amplía la varianza dentro de grupos.
\
Al comparar los factores multiplicativos obtenidos teniendo como base el nivel de
I Igual que para el caso anteri or, se obtu vieron factor es de ajuste multiplicativos máxima producción y los obtenidos al tomar como base la media de mínimos cuadrados
~
ntre la media de mínim os u adrados y la m edia de la subclase a ajustar. U n a vez estimada se puede observar cómo estos últimos son más ineficientes, pues las
2
justados los datos se hiz el aná] ¡ ís de vari anza que ap arece en la Tabla 8. El R = relaciones de F son mucho m a yor es para ca si todos los efectos ( Tabla 9), lo cual
.1 % quedan do u n pequeño efecto residual especialmente p ar~ el factor m ~ ~i~ indica una persistencia en la variación que no es removida por el factor de ajuste.
~mbient al días abier s previo .SIendo l;gurosos en la com paraClOll de esLe anah sls
~on el modelo inmediatamente anterior (Tabla 7 ), puede comprobarse ómo en éste
~aso la r educción de la variación no fué tan completa como e so a n terior p ues e n
este existen r elaciones de ' ", " A Tabla 9. Comparación de las r laciones F a los n iveles de máxima y media producción
para factores de ajuste aditivos y mult iplica ti os.
Tabla 10. Medias y desviaciones estándar para prod ucción por 1 ctancia e n regi stros
sin ajustar y aj ust a dos a di tiva o multiplica tivamente.
68. 1995 Rev. Fac. Na!. Agr. Mede llín. Vol 48 Nos. 1 y 2, p.4 H8. 1995 61
I
Arboleda, E.; Valencia, C. y Montoya, C. Comparación entre factores de ajuste...
Ta bla 11. Desviaciones están dar para las subclases bas y los registros a ajustar. aditivos, concluyendo que estas no son significativas. Fávero et al (1991) trabaj
con cerdos Landrace en Brasil, concluyen que un tipo de ajuste multiplicativo
que mejor comportamiento tiene con sus datos, al medir el espesor de grasa de
SUBCLASE DS (O) aunque los dos métodos aditivos usados ofrecen, así mismo, resultados satisfact
y otro método multiplicativo es descartado por ofrecer medias y varianzas
Edad base 319 .00 diferentes a los anteriores" • .
Promedio demás edades 144.00 De acuerdo con Dalf
m ultiplicativos y no va
CONCLUSIONES
Con la base de datos analizada, los factor s de aj u · te aditivos ofrece n un mejor 1. De a cuerd con 1
comportamiento que los multiplicativos tenie n do com o r efer n cia la relación de F, aditivos resultad
para todos y cada uno de lo efectos e tu dia dos . E sto concuerda con os resultados factore s de ajustJ
obtenidos por Chauhan, 1988, quien halla pequeñas diferencias a favor de los factores I
aditivos, aunque podrían ser despreciables ; sin embargo recomienda su uso por la 2. Los fa ctores de a~
utilización práctica en las regiones Alpinas don de se desarrolló el tra b ajo. máxim a producc~
media de minim~
multiplicativamente (Tabla 10), cuando se usó como base la medí a de mínimos d a tos an !izada
cuadrados, no presentan difer n cias entre sí, observá n dose que ambos tipos de ajuste ajuste a di tivos c
bajan la varianza entre grupos si se les compara co n as lactancias sin ajustar. son muy m arcad1
Al tomar como base el nivel de máxima producción, la desviación estándar se 4. L a decisión sobri
incrementa para factores de ajuste multiplicativos y permanece igual para factores mul tipli cativos dJ
aditivos siendo mayor y menor respectivamente que la desviación estándar para datos sujeta a est
I
62 Rev. Fac. Na!. Agr Medellín. Vol 48 Nos. 1 y 2, p.41·68. 1995 Rev. Fac. Na!. Agr Medellín. Vol 48
I
Arboleda, E .; Valencia , C. y Montoya, C. Comparación entre factores de ajuste...
abla 11. Desviaciones estándar para las ubclases base y los registros a aj u s tar. aditivos, concluyendo que estas no son significativas. Fávero et al (1991) trabajando
con cerdos Landrace en Brasil, concluyen que un tipo de ajuste multiplicativo es el
que mejor comportamiento tiene con sus datos, al medir el espesor de grasa dorsal;
SUBCLASE D.s (O) aunque los dos métodos aditivos usados ofrecen, así mismo, resultados satisfactorios
y otro método multiplicativo es descartado por ofrecer medias y varianzas muy
Edad base 319.00 diferentes a los anteriores y a las medidas observadas.
Promedio demás ed ades 144.00 De acuerdo con Dalton (1980), la manera de escoger entre un factor de ajuste
Año base 107.55 aditivo y multiplicativo, es por medio de la comparación de las desviaciones estándar
' ,., _.--
Promedio rl del registro base y los registros a ajustar. En la presente investigación, podemos
observar (Tabla ll), cómo en algunos de los efectos, principalmente en aquellos donde
se hizo un aj uste multiplicativo y 1 s relacion 's de F fueron mayores, hay diferencias
marcadas en tre las desviaciones estándar del r gis tro base y de los registros a ajustar,
hallándose que la varianza dentro de gru pos siempre se incrementó con los factores
multiplicativos y no varió cuando se usaron factore aditivos.
CONCLUSIONES
mejor 1. De acu erdo con los criterios de comparación estableci'~os, los factores de ajuste
n de F, adi tivo s resultan ligeramente más eficientes en reducir la variación que los
ltados factores de ajuste multiplicati os.
ctores
par la 2. L . fa tores de ajuste m ultiplicativos obten idos al tomar como base el nivel de
máxima pr du ccion son más efectiv s que los obtenidos teniendo como base la
media de mínimos cuadrados.
va o 3. Par efectos prácticos de uso de facto res de aju ste, de acuerdo con la base de
dal s an aliza da en sta investigación, puede trabaj arse tanto con factores de
aj uste ditivos como multiplica tivos, puesto que 1, s diferencias entr ellos no
son muy marcadas .
os .
san BIBLlOGRAFIA
con
jor AGUDELO VELE Z, José Javier y CARDONA JIMENEZ, Jaime Alejandro. Algunos
la factore s no genéticos que afectan la producción por lactancia en un hato
uso Holstein . Medellín, 1991. 56 p. Tesis (Zootecnistas). Uni versidad Nacional
ln de Colombia. Facultad de Ciencias Agropecuarias.
ón
DALTON, D. S. Introducción a la gen ética anim a l práctica. Zaragoza: Acribia, 1980. McD
168 p.
HARVEY, W. R. Users guide for LSMLMW PC-1 versión; mixed modelleast-squares SCE
and maximum likelihood computer program, s.l.: s.n., 1988; 59 p .
OSTA N., P. R. et al. Fatores genéticos e de meio na produ~ao de leite do rebanho - - - o E stim ation in changes in herd environment. En: Journal ofDairy Science.
leiteiro da Universidade Federal de Vi~osa, Estado de Minas Gerais. En: Vol. 32, No. 8 (aug. , 1949); p. 706 .
Revista da Sociedade Brasileira de Zootecnia. Vol. 13, No. 3 (mai.-jun.,
1984); p. 334-346.
- - - et al. T he estimation of genetic nd environmental trends from record s
~
subjected to culling. E n: Biometrics. No . 15, 1959; p. 192. Citado por
CUNDIFF, L. V; WILLHAM, R. L. and PRATT, Charles A. Additive versus
MARTl NEZ, M. L.; LE E A. J. and LI , C. y. Adjustment factors by
multiplicative correction facto ~ f(\r ...-- . - • in beef cattle. En:
m a ximu n likelih ood, gTat;S compa risons, an d paired compa risons. En:
Jour-nal of A -- -, ~ . 983-987.
Jou rnal of D airy Science. Vol. 73, No. 3 (mar., 1990); p. 819-825.
airy records for
fDairy Science. HICKMAN, C. G. Her d-I vel methods for age adj ustm cnt ofmilk yi elds. En: Journal
of Dairy Science. Vol. 56, 0. 7 (jul. , 1 7 ); p. 94 7-951.
tics of friesian MARTINEZ, M. L. ; LEE, A J. and UN, C. y. Adj ustment fa tors by ma.xi mun
Science. Vol. bkeli h ood, 6'l'O SS comparisons, and paired omparisons. En: J ournal ofDai ry
SC! nc . Vol. 73, N o. 3 (mar., 1990); p. 819-825 .
e Antioquia :
tioquia. En: MILLER, Paul. A r cent stud of age adjustm nt. En : Joum al of Dairy Science.
, No. 1 (j un, Vol. 56, No. 7 (j ul. , 197 ); p. 952-958.
TAYLOR, J. M. and HEZEL, L. N . The gro th curve of pigs between 134 and 174 ANEXO 1.
days of age. En: Journal of Anima l Science. Vol. 14, No. 4 (nov., 1955); p.
1113-1139. Tabla 1. Factores de aj uste para el efecto edad al parto.
WILMINK J., B. M. Adjustment of lact ation for age a t calving in relation t o level of
production. En: Livestock Produ tion Sci nce . Vol. 16, N o. 4 (j un., 1987); p. Base: Producción máxima Med ia de mínimos cuadrad os
321-334
EDAD ADIT MUL'
1 (S32 meses) ;.> +411 1.08
2 (33-45 mI"
3 (46- 8
4 (59-71 n
5 (72-84 n
6 (85-97 n
7 (;~98 má
L as subc
Tabl a 2.
1 97~
1978
1977.
66 Rev. Fac. Nal. Agr. Medellín VOl 48 Nos 1 y 2. p.41 -68 . 1995
Arboleda, E.; Valencia, C. y Montoya, C. Comparación entre factores de ajuste ...
TAYLOR, J . M. and HEZEL, L. N. The growth curve of pigs b etween 134 and 174 ANEXO 1.
days of age. En: Journal of Animal Science. Vol. 14, No. 4 (nov., 1955); p.
1113-1139. Tabla l. Factores de ajuste para el efecto edad al parto.
I\V1LMINKJ., B. M. Adjustment ofl act ati n for age at cal ving in relation to level of
production. En: Livestock Production Scienc . Vol. 16, N o. 4 (jun., 1987); p. Base: Producción m áx ima Med ia de mínimos cuadrados
321-334
EDAD ADIT MULT ADIT MULT
1 (~32 meses) a +1597 1. 2 +411 1.08
2 (33-45 m es) a +1564 1.31 +378 1.07
3 (46-58 m es) a +1517 1.30 +331 1.06
4 (59-7 m es) a +1 52 1.2 +166 1.03
5 (72-84 mes) a +1129 1.20 - 57 0.99
6 (85-97 m es) a +1143 1.21 - 43 0.99
7 (~98 m ese. ) b O 1.00 -1186 0.82
66 Rev Fac Na!. Agr. Medelllo. Vo l 48 Nos. 1 y 2, p.41-68. 1995 Rev. Fac. Na!. Agr. Medellín. Vol 48 Nos. 1 y 2, pA1·68 1995 67
Arboleda , E .; Valencia, C. y Montoya, C.
1
Base Producción Máxima Media de mínimos cuadrados 1991
No. Parto ADIT MULT ADIT MULT cul ~
agr
,
1 +1034 1.21 +512 1.10
2
3
+326
O
1.06
1.00
+196
-522
0.97
0.91
dia~
4 +168 1.03 -354 0.94 par]
5 +351 1.06 -171 0.97
6 +1253 1.26 +731 1.15 par
de t
:emp¡
J
Tabla 5. Factores de ajuste para el efecto época de parto.
parl
Base Producción Máxima Media de mínimos cuadrados pro.
ciÓI.
Epoca ADIT MULT ADIT MULT
cal ,