Five lactation curve models were compared for their fitness to milk data of
indigenous and crossbred cows using data from Bako Agricultural Research Center and Debre
Zeit Research Station. The fit of the models was compared using the model Rz value. The
overall mean RZ value was 81.5 t 0.12% and was significantly (p < 0.001) affected by model
type, data type, sire breed, parity, calving season and year, location and the interactions
of model type with these fixed effects. The mean Rz values for the Incomplete Gamma with
two parameters (IG ('n = 1)), Inverse polynomial (IP), Exponential, Incomplete Gamma (IG (b
* 1)) and Parabolic exponential models were 90.1 t 0.26, 90.2 t 0.26, 67.2 t 0.28, 74.7 t
0.26 and 77.8 t 0.26 percent, respectively. Higher values of RZ were obtained for the Debre
Zeit (80.5 %) herd than Bako. The sire breeds Friesian (81.6%) crosses had significantly
the highest RZ value compared to the other sire breeds. Cows that calved during Arfasa
(83.4%) had a higher fit than those in other calving seasons. Among the data types
considered, highest R2 was observed for monthly mean (85.6%) and monthly total (85.6 %)
milk data types. Besides, cows in the sixth parity (81.3 %) and those that calved in 1993
(84.5 %) had the highest Rz value. The highest Rz was obtained for IP fitted to monthly
mean (95.4 %) and monthly total (95.4 %) milk data followed by the IG (b = 1) fitted to
monthly mean (93.2 °h) and monthly total (93.2 %) milk data types. The IP followed by IG (b
= 1) fitted to milk data of the crossbreds had the highest R2. The reverse is true for the
indigenous breeds. Among the model type x parity, model type x season and model type x
location interaction groups, the IP followed by IG (b =1) had the highest Rz for all
parities, calving seasons and locations. From this study it can be concluded that both the
IP and IG (b=1) had the best fit and could be used to describe the lactation curves of both
zebu and crossbred cows under study.