Diversité des systèmes d'élevage bovin laitier et performances animales en région semi aride de Sétif

CHAPITRE IV

ANALYSE DES

PERFORMANCES

DES FEMELLES

LAITIERES

Chapitre IV. Performances animales à l'échelle troupeau et individu

La composante animale du système d'élevage constitue le niveau d'investigation principale du zootechnicien. Ce dernier considère l'animal en région méditerranéenne comme un système biologique placé dans un environnement contraignant et possédant des capacités d'adaptation qui vont lui permettre non seulement d'assurer sa propre vie mais également de maximiser ces fonctions de production (BOCQUIER et al. 2004).

Dans cette partie du travail nous étudions le comportement productif et reproductif d'un matériel animal sélectionné dans un milieu tempéré favorable, transféré vers un milieu plus difficile caractérisé par des contraintes liées à l'environnement et à la conduite. Dans cette approche, nous visons la caractérisation de l'effet des différents facteurs qui influent sur le niveau de performances du matériel animal dans son nouveau contexte d'élevage.

Notre étude analyse l'effet de six facteurs environnementaux (ferme, année et saison) et liés à l'animal (génération, parité et âge) sur les fonctions de reproduction et de production de lait de femelles de race Montbéliarde.

4.1. Performances de reproduction 4.1.1. Age au démarrage en production 4.1.1.1. Performances moyennes

Les génisses nées sous climat semi aride réalisent leur première saillie à l'âge de 779 jours (25 mois), alors que le premier vêlage est obtenu à l'âge de 1052 jours (35 mois).

La répartition en classes (figure 4.1a) montre que seule 10% des génisses entrent en reproduction à un âge inférieur de 18 mois, alors que 60% réalisent leur première fécondation à partir de l'âge de 2 ans. Les résultats ont montré aussi que 64% des femelles mettrent bas à un âge compris entre 2 à 3 ans. En effet, seul 2% des génisses vêlent avant l'âge de 2 ans. En Floride, aux Etats-Unis, 40% des génisses Holstein, qui semblent plus précoces, mettrent bas à moins de 23 mois d'âge (SIMERL et al. 1991) comparé à 27% en Californie (ETTEMA et SANTOS, 2004).

Tableau 4.1 : Performances moyennes des paramètres d'entrée en production

	Nombre	Min	u#177;ó (jours)	Max	Coeff. de variation
AMR (jours) APMB (jours)	267 409	301 608	7 79#177;185 1052#177;189	1411 1733	23,79% 17,95%

AMR : Age de la mise à la reproduction ; APMB : Age de la première mise-bas. a) Age à la mise à première saillie

45

75

60

30

15

0

6

7

30

46

63

69

25

10

7

5

Std. Dev = 184,92 Mean = 777

N = 271,00

AMR

b) Age à la première mise bas

120

100

40

20

80

60

0

26

53

75

111

63

41

13

10

Std. Dev = 188,82 Mean = 1052

N = 409,00

APMB

Figure 4.1 : Répartition des paramètres d'entrée en production en classe

4.1.1.2. Facteurs de variation

L'âge à la première mise bas est sensible aux facteurs environnementaux. Les résultats de l'analyse de la variance (tableau 4.2) indiquent que les facteurs ferme, année et saison de naissance affectent significativement l'âge d'entrée en production plus que les facteurs liés à l'animal (génération).

Tableau 4.2 : Carrés moyens des paramètres de fécondité obtenus à partir de l'analyse de variance

Facteurs de variation	DL	AMR	APMB
Ferme (F)	3	0,27 **	0,07 ns
Génération (G)	4	0,08 ns	0,05 ns
Année de naissance (AN)	16	0,12 **	0,11 ***
Saison de naissance (SN)	3	0,14 **	0,17 **
F*G	12	0,22 **	0,09 ***
F*SN	9	0,14 **	0,07 **
G*SN	12	0,05 ns	0,03 ns
Erreur	323	0,05	0,02

AMR : Age de la mise à la reproduction ; APMB : Age de la première mise-bas ; DL : degré de liberté ; ns : différence non significative ; ** : différence significative à p< 0,01 : *** : différence significative à p< 0,001.

4.1.1.2.1. Facteurs de l'environnement Effet de la ferme

Dans nos conditions d'élevage (tableau 4.3), la similarité des pratiques d'élevages des génisses reflète la faible variabilité de l'âge d'entrée en production, bien qu'il faut mentionner une légère supériorité de la ferme situé le plus au Sud. Dans cette dernière, durant la phase d'élevage des génisses, les pratiques alimentaires se basent sur la distribution du son de blé en plus de la ration de base. Cela reflète l'importance de la conduite de l'élevage des génisses en terme d'alimentation dans la phase d'élevage sur les paramètres d'entée en production (FIEZ, 1993).

Tableau 4.3 : Variabilité des paramètres d'entré en production selon la ferme

Ferme	AMR pió (n) jours		Erreur Std	APMB pió (n) jours		Erreur Std
F1	729#177;188^a	(100)	18,75	1 025#177;181	(213)	12,42
F2	805#177;152^b	(76)	17,47	1074#177;161	(94)	16,61
F3	812#177;191^b	(57)	25,33	1067#177;191	(63)	24,12
F4	776#177;155^ab	(32)	27,31	1042#177;156	(34)	26,77

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; AMR : âge à la mise à la reproduction ; APMB : âge à la première mise bas.

Effet de l'année de naissance

L'année de naissance affecte significativement les paramètres d'entrée en production
(p<0,01). Les différences sont de l'ordre de 209j pour l'âge de la mise à la première saillie et

de 220j pour l'âge à la première mise bas (Figure 4.2). On remarque aussi la similitude de l'évolution des deux paramètres selon l'année.

1600

1400

1200

1000

400

800

600

200

0

AMR APMB

année

Figure 4.2 : Evolution de l'âge d'entrée en production de 1982 à 2002

Effet de la saison de naissance

Le tableau 4.4 montre que les génisses nées en hiver expriment un retard concernant la maturité sexuelle en comparaison avec celles nées durant le reste de l'année. Ce retard peut aller jusqu'à 2 à 3 mois et semble lié aux disponibilités alimentaires après la naissance et au moment de la mise à la reproduction.

Tableau 4.4 : Variabilité des paramètres d'entré en production selon la saison de naissance

Saison de naissance	AMR piu (n) jours		Erreur Std	APMB piu (n) jours		Erreur Std
Hiver	810#177;143^a	(71)	16,93	1096#177;163 ^a	(106)	15,78
Printemps	786#177;212 ^ab	(69)	25,56	1048#177;177 ^b	(103)	17,39
Eté	755#177;1 74 ^ab	(56)	23,21	1023#177;174 ^b	(86)	18,74
Automne	742#177;173 ^b	(69)	20,80	1005#177;184 ^b	(108)	17,72

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; AMR : âge à la mise à la reproduction ; APMB : âge à la première mise bas.

3.1.1.2.2. Facteurs liés à l'animal Effet de la génération

La génération animale n'a pas d'effet significatif (p<0,05) sur l'âge d'entrée en production, la différence n'excède pas 50 jours (tableau 4.5). Néanmoins, cette tendance n'est pas la même dans l'ensemble des fermes. En effet, des résultats supérieurs sont observés dans l'ensemble des fermes chez la génération G1 importée au stade génisse pleine comparée aux femelles nées localement. La 2^ème génération présente un léger retard de l'âge à la première mise bas,

qui s'accentue pour les génération suivantes dans les fermes de la région du Centre (ferme 2 et 3). Dans les fermes situées au Sud (F1) et au Nord (F4) le retard est beaucoup plus important chez la génération 2 mais une nette amélioration est remarquée chez les générations G3 et G4. Les mêmes résultats sont observés par TRACH (2003) chez les vaches importées au Vietnam, dont l'augmentation de l'âge à la première mise - bas est de 110 jours entre la 1ère et la 2^ème génération.

Tableau 4.5 : Variabilité des paramètres d'entré en production selon la génération animale Génération AMR pió (n) jours Erreur Std APMB pió (n) jours Erreur Std

G1				1008#177;184	(67)	22,48
G2	764#177;161	(92)	16,72	1068#177;168	(121)	15,23
G3	788#177;188	(114)	16,71	1050#177;1 78	(135)	15,32
G4	735#177;218	(38)	35,33	1016#177;197	(52)	27,33
G5	760#177;193	(5)	86,45	1046#177;1 79	(6)	73,00

AMR : âge à la mise à la reproduction ; APMB : âge à la première mise bas. ES : erreur standard

1600

1400

1200

1000

400

200

800

600

0

F1 F2 Ferme F3 F4

G1 G2 G3 G4

Figure 4.3 : Variabilité de l'âge à la première mise bas entre fermes selon la génération animale

4.1.2. Paramètres de Fertilité des vaches 4.1.2.1. Performances moyennes

Les vaches laitières en région semi aride réalisent leur premier saillie à 98#177;63 jours post- partum. La fécondation est généralement obtenue après 1,43#177;0,99 services. La figure (4.4 b) montre un intervalle vêlage 1^ère saillie inférieur à 40 jours pour 13% des femelles et moins de 80 jours pour 50% des cas ; en revanche plus de 27% des premières saillies ont lieu après 120 jours post-partum. Au Canada, HAYES et al. (1992), rapportent que seul 0,8% des vaches sont inséminées avant 50j et 76% avant 100j post-partum.

La figure (4.4 a) montre que plus de 77% des fécondations sont obtenus après le 1^er service, alors que 11% des fécondations nécessite trois saillies ou plus. Ces résultats sont meilleurs par rapport à ceux enregistrés en insémination artificielle dans les pays tempérés. En effet, la fécondation après la 1^ère insémination est obtenue pour 49% des saillies en Irlande (BUCKLEY et al. 2003) et pour 60% au Canada (HAYES et al. 1992).

Tableau 4.6 : Performances moyennes de paramètres de fertilité

	nombre	Min	u#177;ó	Max	Coeff. de variation
IC IVS1 (jours)	1875 1290	1 7	1,43#177;0,99 98#177;63	10 418	69,23% 64,57%

IC : Indice coïtal ou nombre de services par conception ; IVS1 : Intervalle vêlage première saillie.

a) indice coïtal

1600

1400

1200

1000

1454

Std. Dey = ,99 Mean = 1,4

N = 1875,00

214

88 80

800

600

400

200

0

IC

b) Intervalle vêlage première saillie

400

500

300

200

100

0

164

444

334

171

82

41

25

Std. Dey = 63,46 Mean = 98

N = 1290,00

IVS1

Figure 4.4 : Répartition en classe des paramètres de fertilité des femelles

4.1.2.2. Facteurs de variation

Les carrées moyens des paramètres de fertilités obtenus à partir de l'analyse de variance sont résumés dans le tableau 4.7. En effet, tous les facteurs de variation sont hautement significatifs (p< 0,01) à l'exception de la saison de vêlage et la parité pour l'indice coïtal.

Tableau 4.7 : Carrés moyens des paramètres de fertilité obtenus à partir de l'analyse de la variance Facteurs Dl IC IVPS

Ferme (F)	3	0,67 *	13,94 ***
Année (An)	16	1,04 ***	1,66 ***
Saison de vêlage (SV)	3	0,50 ns	00,16 ns
Génération (G)	3	0,56 *	2,51 ***
Parité (P)	5	0,12 ns	00,84 *
Age (A)	4	3,17 **	1,32 **
F*An	48	0,74 ***	0,86 ***
F*SV	9	0,37 ns	1,52 ***
F*G	9	0,54 ***	0,44 ns
SV*An	48	0,29 *	0,41 ns
G*P	15	0,28 ns	0,31 ns
G*A	12	0,54 **	0,36 *
G*Sv	9	0,16 ns	0,85 **
P*Sv	15	0,27 ns	0,31 ns
Erreur	1066	0,17	0,292

IC : Indice coïtal ou nombre de services par conception ; IVS1 : Intervalle vêlage première saillie ; Dl : degré de liberté. ns : différence non significative ; * : différence significative à p< 0,05 ; ** : différence significative à p< 0,01 : *** : différence significative à p< 0,001.

4.1.2.2.1. Facteurs de l'environnement Effet de la ferme

Les paramètres de fertilités varient selon l'exploitation agricole et la région dans laquelle se situe. En effet, les meilleurs résultats sont rencontrés chez les exploitations pratiquant la monte naturelle comme mode de reproduction unique (F1, F2, F4). Les variabilités entre fermes sont expliquées aussi par le gradient d'aridité et son impact sur le comportement reproductif des animaux, suite probablement à l'effet des disponibilités en ressources alimentaires et du stress thermique.

Tableau 4.8 : Variabilité des paramètres de fertilité selon la ferme

Ferme	IC piu (n)		Erreur std	IVPS piu (n) jours		Erreur Std
F1	1,43#177;0,93 ^a	(854)	0,031	100#177;53 ^a	(565)	2,23
F2	1,33#177;0,75 ^ab	(488)	0,034	69#177;40 ^b	(349)	2,15
F3	1,49#177;1,00 ^a	(224)	0,067	113#177;60 ^a	(119)	5,52
F4	1,24#177;0,59 ^b	(291)	0,035	102#177;57 ^a	(231)	3,74

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; IC : Indice coïtal ou nombre de services par conception ; IVS1 : Intervalle vêlage première saillie.

Effet de l'année de vêlage

La figure 4.5 montre une amélioration continue de l'intervalle vêlage 1^ère saille au cours de la période de suivi entre 1987 et 2002, à l'exception de la période 1992-1994 où l'IVPS a connu une dégradation suite à des considérations encore mal connues. Toutefois, pour l'indice coïtal l'évolution fait ressortir deux aspects : une augmentation significative de ce paramètre passant de 1,19 en 1986 à 1,99 en 1990, suit un maintien au tour de la moyenne après 1990. Ces résultats reflètent l'attention portée par les gestionnaires des fermes à l'amélioration des paramètres de fertilité.

140

120

100

40

80

60

20

0

IVPS (jours) IC

année

2,2

2,0

0,8

0,6

0,4

0,2

0,0

1,8

1,6

1,4

1,2

1,0

IC

Figure 4.5 : Evolution des paramètres de fertilité au cours de la période de suivi

Effet de l'interaction entre ferme et année de vêlage

La variabilité de l'IVPS selon l'année dépend de la ferme (Figure 4.6). On peut noter trois types d'évolutions : une évolution régulière dans les fermes F1 et F2 avec des différences entre année n'excédant pas les deux mois ; une évolution irrégulière dans la ferme F3 avec des écarts entre année de plus de 6 mois liés probablement aux problèmes sanitaires rencontrés dans la période 94-96; une évolution caractérisée par trois période différentes dans la ferme F4 avec une amélioration de l'IVPS durant la première période allant de 87 à 91, dont l'amélioration est supérieur à 4 mois par rapport à la première année (1987), une chute durant la période 9 1-94 enregistrant une perte de 50 jours, puis une reprise durant la dernière période pour se situer autour de la moyenne générale.

270

240

210

180

150

120

90

30

60

0

F1 F2 F3 F4

Année

Figure 4.6 : Evolution de l'intervalle vêlage première saillie selon la ferme au cours de la période de suivi

Effet de la saison de vêlage

Les résultats obtenus montrent une évolution comparable des paramètres de fertilité au cours de l'année (tableau 4.9). L'indice coïtal suit une évolution décroissante mais non significative de l'hiver à l'été, suit une légère dégradation en automne, qui enregistre une augmentation de 0,17 points probablement à cause des disponibilités alimentaires à cette saison. Pour l'IVPS les variations sont faibles et ne dépassent pas 6 jours.

Tableau 4.9 : Variabilité des paramètres de fertilité selon la saison de vêlage

Saison de vêlage	IC piu (n)	Erreur std	IVPS piu (n) jours		Erreur Std
Hiver	1,46#177;0,88 (384)	0,045	94#177;53	(337)	2,88
Printemps	1,42#177;0,85 (397)	0,043	95#177;58	(353)	3,09
Eté	1,34#177;0,83 (291)	0,049	95#177;54	(262)	3,36
Automne	1,51#177;0,99 (368)	0,052	89#177;48	(312)	2,70

IC : Indice coïtal ou nombre de services par conception ; IVS1 : Intervalle vêlage première saillie.

Effet de l'interaction entre ferme et saison de vêlage

L'évolution des paramètres de fertilité au cours de l'année varie selon la ferme. En effet, l'effet de la saison de vêlage est remarquablement observé dans les fermes F3 et F4, qui enregistrent des variations importantes de stocks fourragers durant l'année.

210

180

150

120

30

90

60

0

F1 F2 F3 F4

Ferme

Hiver Printemps Eté automne

Figure 4.7 : Variabilité des paramètres de fertilité des fermes selon la saison de vêlage

4.1.2.2.2. Facteurs liés à l'animal

Effet de la génération animale

Selon le tableau 4.10, l'effet génération est significatif ; ce facteur agit plus sur l'IVPS
(p<0,001) que sur l'IC (p<0,05). L'indice coïtal le plus faible est observé chez la génération

G3 (1,31 points), exprimant une amélioration continue au cours des générations (G1 à G3). La

G4 manifeste en revanche une détérioration de ce paramètre.

Concernant l'IVPS, les générations nées localement réalisent des performances supérieurs à celles des vaches importées (p<0,001). En effet, on note une évolution décroissante (p<0.05) de ce paramètre de la G1 à la G4, gagnant ainsi 22 jours.

Tableau 4.10 : Variabilité des paramètres de fertilité selon la génération animale Génération IC piu (n) Erreur std IVPS piu (n) Erreur Std

G1	1,45#177;0,89 ^ab	(467)	0,041	103#177;53 ^a	(323)	2,96
G2	1,38#177;0,89 ^b	(602)	0,036	94#177;55 ^b	(417)	2,70
G3	1,31#177;0,75 ^c	(626)	0,030	87#177;52 ^b	(441)	2,48
G4	1,51#177;0,98 ^a	(130)	0,086	81#177;46 ^b	(71)	5,44

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ;IC : Indice coïtal ou nombre de services par conception ; IVS1 : Intervalle vêlage première saillie.

Effet de l'interaction entre génération et saison de vêlage

Les différences entre générations sont plus importantes en saison estivale, traduisant ainsi un
gradient d'adaptation de la G1 à la G4. Durant cette saison, les femelles de la G4 s'adaptent
mieux comparées aux femelles de G1, plus sensibles au stress thermique; par conséquent le

repos post-partum s'allonge de 2 mois pour la G1 par rapport à la G4. Les différences entre générations diminuent en se rapprochant de l'hiver.

180

160

140

120

100

40

80

20

60

0

Hiver Printemps Eté Automne

saison de vêlage

G1 G2 G3 G4

Figure 4.8 : Variabilité de l'intervalle vêlage première saillie entre génération et saison de vêlage

Effet de la parité

Au cours de la carrière des femelles, on n'observe pas des différences significatives de l'IC, confirmant ainsi les observations de HAYES et al. (1992) chez la race Holstein au Canada. Toutefois une légère amélioration de l'IVPS est observée de la 1^ère à la 4^ème mise bas, qui gagne 15 jours.

Au début de la carrière des femelles encore en croissance, un allongement significatif de l'IVPS par rapport au 3^ème et 4^ème parités s'explique par une compétition accrue entre fonctions physiologiques, suit une amélioration continue, signe d'une adaptation progressive jusqu'à la 4^ème parité, et une rechute des performances pour se situer en fin de carrière à un niveau comparable à celui des primipares.

Tableau 4.11 : Evolution des paramètres de fertilité au cours de la carrière

Parité	IC piu (n)	Erreur std	IVPS piu (n) jours		Erreur Std
P1	1,48#177;0,99 (367)	0,052	99#177;54 ^a	(298)	3,13
P2	1,45#177;0,95 (315)	0,053	94#177;52 ^ab	(270)	3,17
P3	1,44#177;0,92 (250)	0,058	91#177;55 ^bc	(221)	3,70
P4	1,46#177;0,86 (196)	0,062	84#177;48 ^c	(169)	3,68
P5	1,29#177;0,60 (143)	0,050	93#177;62 ^abc	(125)	5,56
P6+	1,41#177;0,81 (184)	0,059	94#177;50 ^abc	(178)	3,75

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ;IC : Indice coïtal ou nombre de services par conception ; IVS1 : Intervalle vêlage première saillie.

Effet de l'âge des femelles

La fertilité des femelles varie selon l'âge de la vache (tableau 4.12). L'évolution de l'IVPS exprime des difficultés de reproduction au début de la vie productive, une amélioration jusqu'à l'âge de sept ans, puis une dégradation dès l'âge de 9 ans. En revanche, l'indice coïtal varie d'une manière irrégulière. Les vaches sont fécondées avec moins de 1,4 saillies durant les périodes d'âge 3 - 5 ans et 7 - 9 ans. L'IC s'élève en revanche à plus de 1,5 services par fécondation au début de carrière (moins de 3 ans) et au milieu (5 à 7 ans).

Tableau 4.12 : Variabilité des paramètres de fertilité selon les classes d'âge

Age	IC piu (n)	Erreur std	IVPS piu (n) jours		Erreur Std
A1	1,58#177;1,07 ^a (236)	0,070	104#177;55 ^a	(183)	4,05
A2	1,39#177;0,89 ^b (481)	0,041	94#177;54^b	(421)	2,62
A3	1,51#177;0,90 ^a (371)	0,047	87#177;50 ^b	(322)	2,78
A4	1,30#177;0,72^b (220)	0,048	90#177;56 ^b	(199)	3,95
A5	1,39#177;0,74 ^ab (135)	0,064	98#177;55 ^ab	(137)	4,70

A1: =3 ans; A2: 3 à 5 ans; A3: 5 à 7 ans; A4: 7 à 9 ans; A5: >9 ans.

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; IC : Indice coïtal ou nombre de services par conception ; IVS1 : Intervalle vêlage première saillie.

4.1.2.3. Impact de l'allongement de l'intervalle vêlage première saillie sur l'indice coïtal

La fertilité des vaches est fortement influencée par la durée du repos post-partum (Figure 4.9). Une faible fertilité traduite par un IC élevé (> 1,5) est observée si l'IVPS est court (<45j). Après 45 jours, la variabilité de l'IC est faible, moins de 0,1 points entre les différentes classes de repos post-partum (1,19 et 1,29). En Irlande, BERRY et al. (2003) rapportent une diminution des possibilités de fécondation à la première insémination dans le cas d'un IVPS court, et SMITH et al. (2002), ont observé une excellente détection de chaleur lorsque l'IVPS est en moyenne de 73j.

A B C D E

IVPS (jours post-partum)

2,60

2,30

2,00

1,70

1,40

1,10

0,80

0,50

a

a

a

a

b

Figure 4.9 : Effet de la durée de repos post-partum sur l'indice coïtal

lettres différentes exprimant des différences significatives (p< 0.05) ; A : IVPS< 45j, B : IVPS 45-90j, C : IVPS 90-120j ; D : IVPS 120-150j ; E : IVPS + 150j.

Effet de la génération

pour l'ensemble des générations, l'IC s'améliore au fur et à mesure que le repos post-partum augmente pour atteindre un optimum lorsque l'IVPS est compris entre 90j et 120j, au delà duquel la fertilité accuse des variations liées probablement à des problèmes de reproduction.

0,90

0,70

0,50

1,90

1,70

1,50

1,30

1,10

A B C D E

IVPS (jours post-partum)

G1 G2 G3 G4

Figure 4.10 : Effet de la durée de repos post-partum sur l'indice coïtal selon la génération A : IVPS< 45f, B : IVPS 45-90f, C : IVPS 90-120f ; D : IVPS 120-150f ; E : IVPS +150f.

Effet de la parité

Au cours de leurs carrières, les vaches inséminées tôt, à moins de 45j post-partum, rencontrent des problèmes de `Repeat breeding', et exigent plus de 1,5 saillies pour être fécondées. En effet, les meilleurs résultats sont observés pour un repos post-partum compris entre 45j et 1 20j.

1,80

1,60

1,40

1,20

1,00

P1 P2 P3 P4+

A B C D E

IVPS (jours post-partum

Figure 4.11 : Effet de la durée de reps post-partum sur l'IC selon la parité A : IVPS< 45f, B : IVPS 45-90f, C : IVPS 90-120f ; D : IVPS 120-150f ; E : IVPS +150f.

4.1.3. Paramètres de fécondité 4.1.3.1. Performances moyennes

La fécondité exprime dans le présent travail l'intervalle vêlage -fécondation (IVF) et intervalle entre mise bas (IMB). La moyenne de l'IVF dans les quatre fermes est de 125#177;87 jours, alors que celle de l'IMB est évaluée à 413#177;100 jours.

Selon la figure 4.12 a, 15% des fécondations ont lieu à moins de 50 jours post-partum et 50% à moins de 100 jours. Une partie des fécondations, estimée à 16%, n'a pu avoir lieu qu'après cinq mois post-partum. Au Canada, HAYES et al. (1992) ont noté que chez la race Holstein seul 4% des femelles sont fécondées avant 50jours post-partum et 57% avant 100 jours, alors qu'au Royaume Unis, 73% des vaches ont été observées gestantes entre 40 et 120 jours post- partum (ESSOLEMONT et KOSSAIBATI, 2001).

Concernant l'intervalle entre mise bas, 48% se situent à moins de 380 jours (12,5 mois), alors que 15% dépassent les 15 mois (figure 4.12 b).

Tableau 4.13 : Performances moyennes des paramètres de fécondité

	Nombre	min	u#177;ó	Max	CV %
IVF (j) IMB (j)	1321 1419	10 242	125#177;8 7 413#177;100	561 903	69,93% 24,12%

IVF : intervalle vêlage fécondation ; IMB : intervalle entre deux mise bas.

a) Intervalle vêlage fécondation

400

500

300

200

100

0

196

462

294

163

85

47

30

Std. Dey = 87,40 Mean = 125

N = 1321,00

IVF

b) intervalle entre mise bas

400

600

500

300

200

100

0

117

565

354

167

82

53

35

Std. Dey = 99,52 Mean = 413

N = 1419,00

IMB

Figure 4.12 : Répartition en classe de la fécondité des femelles

4.1.3.2. Facteurs de variation

Les paramètres de fécondité sont sensiblement affectés par les effets des facteurs de l'environnement et les facteurs liés à l'animal. En effet, la ferme, l'année, la génération et l'âge représentent les facteurs de variabilité les plus significatifs (p<0,001). Notons aussi l'existence de différences liées aux interactions entre facteurs (tableau 4.14).

Tableau 4.14 : Carrés moyens des paramètres de fécondité obtenus à partir de l'analyse de variance Facteurs Dl IVF IMB

Ferme (F)	3	9,49 ***	1,12 ***
Année (An)	15	2,82 ***	0,29 ***
Saison de vêlage (SV)	3	0,36 ns	0,07 ns
Génération (G)	3	6,53 ***	0,60 ***
Parité (P)	5	0,58 ns	0,10 *
Age (A)	4	1,47 **	0,22 ***
F*An	45	1,05 ***	0,12 ***
F*SV	9	0,64 ns	0,04 ns
F*G	9	1,25 **	0,16 ***
SV*An	45	0,37 ns	0,05 *
G*P	15	0,40 ns	0,04 ns
G*A	12	0,47 ns	0,04 ns
G*Sv	9	0,57 ns	0,06 ns
P*Sv	15	0,29 ns	0,02 ns
Erreur	1199	0,34	0,03

IVF : intervalle vêlage fécondation ; IMB : intervalle entre mise bas ; ns : différence non significative ; * : différence significative à p< 0,05 ; ** : différence significative à p< 0,01 : *** : différence significative à p<0,001.

4.1.3.2.1. Facteurs de l'environnement Effet de la ferme

Selon le tableau 4.15, des différences entre fermes sont observées. Pour les deux paramètres de fécondité les meilleurs records sont enregistrés à la ferme F2 dont les performances se rapprochent aux normes couramment admises en reproduction bovine. Dans les autres fermes, les vaches réalisent deux mise- bas en 14 mois en moyenne. Ces variations peuvent être expliquées par la localisation de la ferme sur le gradient de l'aridité, qui détermine le niveau de ressources alimentaires disponibles, et les pratiques de reproduction.

Tableau 4.15 : Variabilité des paramètres de fécondité selon la ferme

Ferme	IVF piu (n) jours		Erreur Std	IMB piu (n) jours		Erreur Std
F1	127#177;78 ^a	(594)	3,20	411#177;89 ^a	(648)	3,52
F2	93#177;61 ^b	(355)	3,26	373#177;61 ^b	(362)	3,19
F3	136#177;74 ^a	(127)	6,55	444#177;96 ^c	(151)	7,85
F4	130#177;79 ^a	(220)	5,31	411#177;86 ^a	(233)	5,62

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; IVF : intervalle vêlage fécondation ; IMB : intervalle entre deux mises bas.

Effet de l'année de vêlage

L'évolution des paramètres de fécondité selon l'année est illustrée dans la figure 4.13. On remarque une améliorations progressive sur 15 ans, qui s'est traduite par des gains de 87j et 100j respectivement pour l'IVF et IMB. Ces résultats attestent des efforts d'amélioration de la fécondité par la réduction du nombre des jours ouverts.

400

300

500

200

600

100

0

IVF (jours) IMB (jours)

année

Figure 4.13 : Evolution des paramètres de fécondité selon l'année

Effet de l'interaction entre ferme et année

Selon la ferme, trois profils d'évolution des paramètres de fécondité peuvent être décrites :

- Deux phases d'évolution à la ferme F1 : des intervalles prolongées observées durant la période 1987/1993 allant de 130 à 170 jours et de 430 à 470 jours respectivement pour l'IVF et IMB, suit un raccourcissement des intervalles pendant une 2^ème période (1994/2002) pour se situer à moins de 120 j et 390 j respectivement pour l'IVF et l'IMB.

- Une évolution régulière dans la ferme F2, les écarts entre les extrêmes ne dépasse pas
les 40 j pour les deux paramètres montrant une maîtrise de la conduite de l'élevage.

- Une évolution irrégulière dans les fermes F3 et F4 exprimant un manque de maîtrise. Les deux fermes en question possèdent de faibles superficies fourragères et accusent des variations inter -annuelles importantes des stocks fourragers.

450

400

350

300

250

650

600

550

500

700

F1 F2 F3 F4

année

Figure 4.14 : Variabilité de l'intervalle entre mise bas selon la ferme et l'année

300

270

240

210

180

150

120

90

30

60

0

F1 F2 F3 F4

Année

Figure 4.15 : Variabilité de l'intervalle vêlage fécondation selon la ferme et l'année

Effet de la saison de vêlage

La fécondité n'est pas affectée par la saison de vêlage. Néanmoins, une certaine différence bien que non significative est enregistrée, traduisant ainsi une évolution décroissante de l'été vers l'hiver et le printemps (tableau 4.16). La femelle peut gagner de 1 1j à 14j pour les deux paramètres en cas de vêlage d'été.

Tableau 4.16 : Variabilité des paramètres de fécondité selon la saison de vêlage Ferme IVF piu (n) Erreur std IMB piu (n) Erreur std

Hiver	122#177; 76	(351)	4,07	411#177;89	(377)	4,56
Printemps	124#177;83	(360)	4,35	408#177;93	(385)	4,72
Eté	113#177;68	(263)	4,18	397#177;78	(283)	4,62
Automne	114#177; 71	(322)	3,94	401#177;80	(344)	4,30

IVF : intervalle vêlage fécondation ; IMB : intervalle entre deux mises bas

Effet de l'interaction entre saison de vêlage et année

Des différences significatives (p<0,05) de l'évolution de l'IMB selon la saison et la compagne agricole sont observées. Deux périodes peuvent être illustrées (figure 4.16) : une variabilité importante entre saisons durant la période comprise entre 1987 et 1991, allant jusqu'à 5 mois en 1989 entre les vêlages du printemps et de l'été en faveur de ce dernier. La deuxième période débute en 1992 est caractérisée par des variations non significatives entre saisons, oscillant entre 6 et 45 jours.

450

400

550

500

350

300

600

Hiver Printemps Eté Automne

Année

Figure 4.16 : Evolution de l'intervalle entre mise bas entre saisons selon l'année de vêlage

4.1.3.2.2. Facteurs liés à l'animal Effet de la génération animale

Le tableau (4.17) montre une différence très hautement significative (p<0,001) des paramètres de fécondité liées à la génération animale. En effet, les générations nées localement réalisent des performances supérieures à celles des vaches importées. Pour les deux paramètres la G3 comparée à la G1 réalise les intervalles les plus courts, résultats d'une amélioration continue des performances de reproduction, qui atteignent respectivement pour l'IVF et l'IMB un gain de 34j et 40j.

Tableau 4.17 : Variabilité des paramètres de fécondité selon la génération animale
Génération IVF u#177;ó (n) jours			Erreur std	IMB u#177;ó (n) jours		Erreur std
G1	142#177;84 ^a	(342)	4,53	431#177;97 ^a	(356)	5,12
G2	112#177;69 ^b	(431)	3,34	401#177;83 ^b	(466)	3,84
G3	108#177;70^b	(441)	3,34	390#177;76 ^b	(477)	3,48
G4	115#177;74 ^b	(75)	8,49	401#177;78 ^b	(83)	8,55

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; IVF : intervalle vêlage fécondation ; IMB : intervalle entre deux mises bas

Effet de l'interaction entre ferme et génération

Les figures 4.17 et 4.18 montrent des évolutions différentes au cours des générations des paramètres de fécondité liées à l'effet troupeau. En effet, les variations les plus importantes sont observées à la ferme F1 et F3 avec deux tendances distinctes. Une amélioration à la ferme F1 (gain respectifs de 57j et 62j pour l'IVF et l'IMB) et une détérioration à la ferme F3 (perte de 96j et 80j pour l'IVF et l'IMB). Les fermes F2 et F4 enregistrent des écarts très réduits entre générations, compris entre 7 à 19 j pour l'IVF, et entre 13 à 20 j pour l'IMB. BONCZEK et al. (1992) rapportent le faible effet de la génération animale sur la variabilité des performances de reproduction chez la race jersiaise aux Etats-Unis.

300

250

200

150

100

50

0

F1 F2 Ferme F3 F4

G1 G2 G3 G4

Figure 4.17 : Variabilité de l'intervalle vêlage fécondation selon la génération et la ferme

400

500

300

200

600

100

0

F1 F2 Ferme F3 F4

G1 G2 G3 G4

Figure 4.18 : Variabilité de l'intervalle entre mise-bas selon la génération et la ferme

Effet de la parité

Au cours de la carrière des femelles l'amélioration de l'IVF est non significative bien qu'elle atteint 18 jours entre la 1^ère et la 4^ème mise -bas, suit un léger allongement en fin de carrière (tableau 4.18). En revanche, l'intervalle entre mise bas est significativement plus long en première parité comparée aux parités qui suivent (p<0,05) ; de la 1^ère à la 5^ème mise bas, les vaches gagnent en moyenne 28j d'intervalle, alors qu'après la 6^ème mise bas la diminution du taux de fécondité est liée au vieillissement des vaches.

Tableau 4.18 : Evolution des paramètres de fécondation au cours de la carrière des femelles

Parité	IVF pió (n) jours		Erreur std	IMB pió (n) jours		Erreur std
P1	127#177;81	(317)	4,56	418#177;95 ^a	(345)	5,11
P2	120#177;75	(285)	4,47	403#177;85 ^b	(301)	4,88
P3	118#177; 78	(232)	5,11	404#177;89 ^b	(250)	5,65
P4	109#177;67	(168)	5,13	397#177;81 ^b	(185)	5,96
P5	112#177;71	(127)	6,34	390#177;69 ^b	(133)	5,96
P6+	117#177;68	(158)	5,41	405#177; 75 ^ab	(175)	5,69

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; IMB : intervalle entre mise bas ; IVF : intervalle vêlage fécondation.

Effet de l'âge des femelles

Le tableau 4.19 montre qu'à un âge précoce les vaches manifestent des difficultés d'adaptation fonctionnelle de la reproduction traduite par des intervalles prolongées. Les paramètres de fécondité connaissent ensuite des améliorations continues pour se situer autours des normes à l'âge de 7 à 9 ans (tableau 4.19).

Tableau 4.19 : Evolution des paramètres de fécondation selon l'âge

Age	IVF piu (n) jours		Erreur Std	IMB piu (n) jours		Erreur Std
A1	137#177;86 ^a	(197)	6,11	431#177;101 ^a	(218)	6,85
A2	118#177;77 ^b	(446)	3,63	402#177;86 ^b	(476)	3,93
A3	114#177;71 ^b	(329)	3,89	399#177;83 ^b	(357)	4,37
A4	107#177;66 ^b	(198)	4,67	395#177;75 ^b	(210)	5,17
A5	123#177;73 ^ab	(121)	3,63	407#177;72 ^b	(127)	6,40

A1: =3 ans; A2: 3-5 ans; A3: 5-7 ans; A4: 7-9 ans; A5: >9 ans.

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; IMB : intervalle entre mise bas ; IVF : intervalle vêlage fécondation.

4.1.3.3. Effet de l'âge à la première mise bas sur les performances de reproduction

L'âge à la première mise- bas (APMB) influence significativement les performances post- partum de reproduction. Cette variabilité est remarquablement observée durant la première mise- bas (Tableau 4.20). En effet, on note de meilleurs résultats lorsque l'APMB est compris entre 36 et 40 mois. En comparant avec les vêlages de moins de 25 mois, les gains réalisés en première parité se situent à 2 et 3 mois respectivement pour l'IVF et l'IMB. A partir de la 2ème parité la variabilité diminue significativement.

Pour certains auteurs réduire l'âge d'entrée en production peut augmenter le nombre de naissance par vache (PIRLO, 2000). En situation semi aride algérienne, l'âge à la première mise bas n'a pas d'effet significatif sur la productivité numérique des femelles (tableau 4.20). en outre, plus les génisses mettant bas tard, à un âge de 36 à 40 mois, plus celles-ci peuvent avoir un veau en plus durant leurs carrières productives par rapport à celles ayant un âge à la première mise bas de moins de 25 mois. Dans le système d'élevage pratiqué en région semi aride caractérisé par la conservation des femelles jusqu'à un âge avancé, le gain lié à la réduction de l'âge à la première mise bas (coût d'élevage des génisses) est perdu par diminution du nombre de naissance par carrière, suite à la réduction de la durée de la vie productive, et par conséquent augmentation des frais d'amortissement des vaches.

Tableau 4.20 : Evolution des paramètres de reproduction durant la carrière selon l'age de la 1ère mise bas.

APMB	< 25 mois	25-30 mois	30-36mois	36-40 mois	> 40 mois	p
IMB
P1 '	497#177;154(10) c	460#177;128 (64) bc	433#177;109(137) ac	404#177;79(52) a	418#177;112(46) ab	0,024
P2	380#177;48(9) a	418#177;107(57) a	413#177;98(113) a	406#177;104(42) a	409#177;86(37) a	0,849
P3	390#177;69(15) ab	419#177;109(136) b	412#177;101 (2 69) b	405#177;78 (119) ab	389#177;89(100) a	0,168
IVF
P1 '	180#177;143(8) ab	178#177;134(63) b	144#177;96(122) a	116#177;67(48) a	128#177;94(44) a	0,016
P2	90#177;27(8) a	133#177;95 (52) a	138#177;93(108) a	118#177;96(42) a	129#177;82(38) a	0,531
P3+	113#177;72(13) ab	134#177;101(117) b	118#177;75(236) ab	120#177;80(117) ab	106#177;64(95) a	0,154
IC
P1	1,70#177;1,25(10)a	1,62#177;1,18(69) a	1,58#177;1,06(143) a	1,43#177;1,11(54) a	1,67#177;1,41 (49) a	0,830
P2	1,44#177;1,01(9) ab	1,93#177;1,44(58) b	1,51#177;1,13(117) a	1,41#177;1,05(46) a	1,48#177;0,89 (42) ab	0,131
P3+	1,07#177;0,26(15) a	1,5 7#177;1, 05(138)b	1,42#177;0,84(259) ab	1,43#177;0,85(128) ab	1,36#177;0,84 (112) ab	0,160
NNC	3,86#177;2,54(7) a	4,67#177;2,39 (55) a	4,68#177;2,51 (114) a	4,74#177;2,91 (47) a	4,44#177;2,95 (43) a	0.914

Lettres différentes sur la même ligne expriment une différence significative à p< 0,05 ; APMB : âge de la première mise bas ; IMB : intervalle entre mise bas ; IVF : intervalle vêlage fécondation ; IC : indice coïtal ; NNC : nombre de naissance par carrière.

4.1.4. Durée de gestation 4.1.4.1. Performance moyenne

La durée de gestation chez les vaches Montbéliardes en situation semi aride est en moyenne de 278#177;14 j, et une fréquence plus élevée en classe modale de 270 à 280j, qui regroupe plus de 53% des résultats. Les vaches de race Montbéliard se placent ainsi dans le groupe des races à durée de gestation moyenne (BOUJENANE, 1983).

Tableau 4.21 : Performances moyennes de la durée de gestation

nombre Min u#177;ó Max Coeff. de variation

DDG (jours) 1209 210 278#177;14 331 4,89%

DDG : durée de gestation

700

600

500

400

300

200

645

317

Std. Dev = 13,63 Mean = 278

N = 1209,00

84

36

100

0

DDG

Figure 4.19 : Répartition de la durée de gestation en classes

4.1.4.2. Facteurs de variation

La durée de gestation des vaches est soumise à l'effet des facteurs de l'environnement. La ferme, l'année et la saison de vêlage déterminent la variabilité de ce paramètre (p<0,001), alors que les facteurs liés à l'animal sont sans effet significatif (p<0,05) à l'exception de la génération animale et ses interactions avec quelques facteurs (tableau 4.22).

	Chapitre IV : Performances animales à l'échelle troupeau et individu
Tableau 4.22 : Carrés moyens de la durée de gestation obtenus à partir de l'analyse de la variance
Facteurs de variation	DL	DDG
Ferme (F)	3	0,03 ***
Année (An)	15	0,006 ***
Saison de vêlage (SV)	3	0,005 *
Génération (G)	3	0,004 *
Parité (P)	5	0,002 ns
Age (A)	4	0,002 ns
Sexe de porté (SP)	2	0,003 ns
F*An	45	0,004 ***
F*SV	9	0,003 *
SV*An	45	0,001 ns
F*G	9	0,004 **
G*P	15	0,003 *
G*A	12	0,003 *
G*SP	6	0,002 ns
G*SV	9	0,001 ns
P*SV	15	0,002 ns
Erreur	946	0,001

DDG : durée de gestation ; DL : degré de liberté ; ns : différence non significative ; * : différence significative à p< 0,05 ; ** : différence significative à p< 0,01 : *** : différence significative à p< 0,001.

4.1.4.2.1. Facteurs liés à l'environnement Effet ferme

Selon le Tableau (4.23), la durée de gestation est influencée par le facteur ferme (p<0.001). Celle-ci est plus longue de 3 à 5 jours chez les femelles de la ferme F2 par rapport aux autres, et semble affectée par le format des vaches.

Tableau 4.23 : Variabilité de la durée de gestation selon la ferme

Ferme	DDG u#177;ó (n) jours		Erreur std
F1	277#177;10 ^a	(550)	0,42
F2	282#177;11 ^b	(310)	0,64
F3	2 79#177;13 ^a	(113)	1,22
F4	277#177;10 ^a	(204)	0,65

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; DDG : durée de gestation

Effet de l'année de vêlage

De 1987 à 2002, on distingue trois période d'évolution de la durée de gestation (Figure 4.20) : une augmentation de la durée de gestation pour se situer à plus de 280j en 1990 ; suit un maintien de la durée de gestation à 280j durant la période 1991/1999 et un retour à un niveau supérieur à 280j après 1999. Ceci peut être lié aux disponibilités des ressources et à la proportion des différentes générations dans le troupeau.

310

305

300

295

290

285

280

275

270

265

260

année

DDG (jours)

Figure 4.20 : Evolution de la durée de gestation selon l'année de vêlage

Effet de l'interaction entre ferme et année

La figure 4.21 montre que seule la période 1994-1996 accuse des variations significatives entre fermes. Ces variations sont liées probablement aux problèmes sanitaires rencontrés dans la ferme F3 durant cette période.

305

300

295

290

285

280

275

270

265

F1 F2 F3 F4

année

Figure 4.21 : Evolution de la durée de gestation selon l'année dans quatre fermes

Effet de la saison de vêlage

La durée de gestation varie significativement selon la saison de vêlage. Elle est plus courte pour le vêlage d'été et plus longue pour le vêlage d'hiver. Nos résultats sont en accord avec ceux de BOUJENANE (1983) et GOYACHE et al. (2002) et en désaccord avec ceux de ECHTERNKAMP et GREGORY (1999) qui rapportent un effet non significatif de la saison de vêlage sur la durée de gestation.

Tableau 4.24 : Variabilité de la durée de gestation selon la saison
Saison DDG u#177;ó (n) jours			Erreur Std
Hiver	280#177;11 ^a	(308)	0,62
Printemps	2 79#177;11 ^ab	(340)	0,59
Eté	277#177;10 ^b	(238)	0,64
Automne	2 79#177;11 ^ab	(291)	0,63

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; DDG: durée de gestation

Effet de l'interaction entre ferme et saison

Selon la figure 4.22, la saison de vêlage affecte sensiblement la durée de gestation des femelles élevées dans les fermes F1 et F3, des meilleurs résultats sont observés pour les vêlages d'été à la F1 et les vêlages d'automne à la F3. Les écarts peuvent aller de 4 à 5 jours. Dans les fermes F2 et F4, la variabilité inter saison est de moins d'une journée.

300

295

290

285

280

275

270

265

260

F1 F2 F3 F4

ferme

Hiver Printemps Eté Automne

Figure 4.22 : Variabilité de la durée de gestation selon la ferme et la saison de vêlage

4.1.4.2.2. Facteurs liée à l'animal

Effet de la génération animale

Selon le tableau 4.25 la durée de gestation chez les vaches de 2^ème génération est significativement plus courte. Les vaches importées réalisent une durée de gestation intermédiaire et comparable à celle des femelles de G3 et G4.

Tableau 4.25	: Variabilité de la durée de gestation selon la génération animale
Génération	DDG u#177;ó (n) jours	Erreur std
G1	2 79#177;11 ^ab (312)	0,63
G2	2 78#177;10 ^a (396)	0,50
G3	280#177;11 ^b (396)	0,56
G4	281#177;10 ^b (66)	1,29

Effet de l'interaction entre ferme et génération

Dans les ferme F1 et F4, la durée de gestation diminue de la G1 à la G3, qui réalisent une baisse de 3 jours en moyenne, alors que dans la ferme F3 la durée de gestation s'allonge chez les générations nées localement comparée à la génération importée (Figure 4.23). Dans cette ferme les écarts entre la 1^ère et la 3^ème génération s'élèvent à plus de 10 jours. A la ferme F2 l'écart entre génération est d'une journée.

300

295

290

285

280

275

270

265

260

F1 F2 F3 F4

Ferme

G1 G2 G3 G4

Figure 4.23 : Variabilité de la durée de gestation selon la ferme et la génération animale

Effet de la parité

La variabilité de la durée de gestation selon l'ordre de vêlage est non significative et les écarts entre les extrêmes ne dépassent pas les 2 jours (tableau 4.26). Néanmoins, la durée de gestation est légèrement plus courte chez les primipares et plus prolongée durant le reste de la carrière. GOYACHE et al (2002) rapportent l'effet significatif de l'ordre de vêlage sur la durée de gestation, qui augmente de 0,5 jours par parité, et gagne à la 9^ème mise bas 5jours.

Tableau 4.26 : Evolution de la durée de gestation au cours de la carrière des femelles

Parité	DDG u#177;ó (n) jours		Erreur std
P1	2 78#177;10	(308)	0,55
P2	2 79#177;11	(269)	0,65
P3	2 79#177;12	(214)	0,79
P4	2 78#177;11	(148)	0,90
P5	280#177;11	(106)	1,07
P6+	280#177;11	(132)	0,96

DDG : durée de gestation

Effet de l'interaction entre parité et génération

Selon la génération animale, la durée de gestation évolue d'une manière différente au cours de la carrière. Les primipares des générations G1, G2 et G3 réalisent des durées de gestation plus courtes par rapport à celles de la génération G4. A partir de la 2^ème parité, les vaches de la G1 ont tendance de prolonger leurs durées de gestation, qui atteigne plus de 282 jours à la 3ème mise bas. Cette durée chute aux 4^ème et 5^ème parité à moins de 276 jours pour augmenter encore en fin de carrière. Concernant la 2^ème génération, la durée de gestation se maintien au dessus de 277 jours jusqu'à la 4^ème mise bas. En fin de carrière la durée de gestation se situe à plus de 280 jours. Chez la 3^ème génération, la durée de gestation varie dans un intervalle étroit, compris entre 279 et 280 jours durant toute la carrière. Pour la génération G4, la durée de gestation est au dessus de 281jours pendant les trois premières mise bas, accuse une chute au 4ème vêlage pour se situer à moins de 277j ours, puis un regain pour atteindre le niveau initial en fin de carrière.

282

280

278

276

G1 G2 G3 G4

284

274

272

270

0 1 2 3 4 5 6 7

Ordre de vêlage

Figure 4.24 : Evolution de la durée de gestation au cours du cycle vitale selon la génération animale

Effet de l'âge des femelles

La durée de gestation ne connaît pas des variabilités significatives selon les différentes classes d'âge. Bien que les variations soient minimes, ont note un prolongement de la durée de gestation de 2 jours entre 3 et 9 ans. L'effet de l'âge sur la durée de gestation est plus remarqué chez les races à viande où l'augmentation est de l'ordre de 3 jours entre le début et la fin de carrière (ECHTERNKAMP et GREGORY, 1999 ; PIEDRAFITA et al. 2000).

Tableau 4.27 : Variabilité de la durée de gestation selon l'âge des femelles

Age	DDG u#177;ó (jours)		Erreur std
A1	2 77#177;10	(191)	0,72
A2	2 79#177;11	(423)	0,52
A3	2 79#177;11	(295)	0,61
A4	279#177;11	(172)	0,86
A5	280#177;11	(96)	1,17

A1: =3ans; A2: 3-5ans; A3: 5-7ans; A4: 7-9ans; A5: >9zns ; DDG : durée de gestation. Effet de l'interaction entre l'âge et la génération

L'évolution de la durée de gestation selon l'âge diffère d'une génération à une autre (figure 4.25). La durée de gestation se situe au dessous de 280 j chez les générations G1, G2 et G3 qu'elle que soit la classe d'âge. En effet, la DDG est en moyenne de 277jours pour les vêlages de vaches âgées de moins de 3 ans, suit un allongement seulement chez les générations G1 et G3 pour se situer à 279j à l'âge de 9 ans. La génération G2 enregistre un raccourcissement de la DDG à l'âge de 5 ans, puis un allongement à l'âge avancé. Chez la 4^ème génération, la DDG est toujours au dessus de 282 jours sauf durant la période d'âge comprise entre 5 et 7 ans, où on enregistre une baisse pour se situer à moins de 277 jours.

272

270

G1 G2 G3 G4

286

284

282

280

278

276

274

0 1 2 3Age 4 5 6

Figure 4.25 : Evolution de la durée de gestation selon l'age et la génération animale

Effet du sexe du veau

Le sexe de la portée n'a pas d'effet significatif sur la durée de gestation (Tableau 4.28). En revanche la gestation gémellaire dure moins qu'une portée simple de sexe mâle, alors que la gestation dont la portée est de sexe féminin est intermédiaire.

Dans la littérature, l'effet du sexe du porté est controversé. Pour GOYACHE et al. (2002) et ECHTERNKAMP et GREGORY (1999), le sexe agit d'une façon significative sur la duré de gestation. La vache retarde la parturition plus d'une journée si elle porte un veau mâle. Pour d'autre (BOUJENANE, 1983 ; PIEDRAFITA et al. 2000), le sexe apparaît sans effet sur la durée de gestation.

Tableau 4.28 : Variabilité de la durée de gestation selon le sexe Nouveau né DDG u#177;ó (n) jours Erreur std

male	2 79#177;11 ^a (561)	0,47
femelle	2 79#177;10 ^ab (533)	0,42
jumeaux	2 76#177;11 ^b (44)	1,72

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; DDG : durée de gestation

Effet du sexe des naissances gémellaires

Concernant les jumeaux, la gestation d'une portée de sexe femelle est significativement la plus longue; la précocité de la parturition chez une portée gémellaire de sexes différents est de 6 à 8 jours par rapport à un porté de même sexe (Tableau 4.29). ECHTERNKAMP et GREGORY (1999) rapportent l'effet non significatif des sexes de porté gémellaire sur la durée de gestation.

Tableau 4.29 : Variabilité de la durée de gestation des jumeaux selon le sexe Jumeaux DDG u#177;ó (n) jours Erreur std

Femelle - femelle	280#177;12 ^a	(11)	3,73
Mâle - mâle	278#177;12 ^ab	(14)	3,19
Mâle - femelle	272#177;9 ^b	(19)	2,15

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; DDG : durée de gestation

4.2. Analyse des performances de production du lait à l'échelle troupeau et individu 4.2.1. Durée de lactation

4.2.1.1. Performance moyenne

Dans la région semi aride la durée de lactation moyenne est de 292 jours (Tableau 4.30), mais le coefficient de variation est assez élevé (23%). La répartition en classe de 40 jours d'intervalle montre que plus de 34% des lactations durent entre 270 à 310 jours (Figure 4.26), alors que 20% des lactations durent plus d'un ans (10%), ou bien moins de 7 mois (10%).

Tableau 4.30 : Performances moyennes de la durée de lactation

nombre Min u#177;ó Max Coeff. de variation

DDL (jours) 943 122 292#177;66 600 22,76%

DDL : durée de lactation

400

300

200

100

0

29

97

209

320

129

83

40

Std. Dey = 66,41 Mean = 292

N = 943,00

DDL

Figure 4.26 : répartition en classe de la durée de lactation

4.2.1.2. Facteurs de variation

L'analyse de la variance montre que la durée de lactation est sous l'influence de plusieurs facteurs, dont les effets sont différents (Tableau 4.31). Les effets ferme, année et génération animale sont les plus significatifs (p< 0,001). La durée de lactation varie aussi selon l'âge de l'animal (0,05), alors que l'ordre de parité et la saison de vêlage apparaissent sans effet sur la durée de lactation.

	Chapitre IV : Performances animales à l'échelle troupeau et individu
Tableau 4.31 : Carré moyen de la durée de lactation obtenu à partir de l'analyse de la variance
Facteurs de variation	DL	DDL (jours))
Ferme (F)	3	0,20 **
Année (An)	7	0,39 ***
Saison de vêlage (SV)	3	0,05 ns
Génération (G)	3	0,36 ***
Parité (P)	5	0,07 ns
Age (A)	4	0,13 *
F*An	21	0,13 ***
F*SV	9	0,11 **
SV*An	21	0,05 ns
F*G	9	0,05 ns
G*P	15	0,07 *
G*A	12	0,07 *
G*SV	9	0,03 ns
P*SV	15	0,05 ns
Erreur	713	0,04

DDL : durée de lactation ; DL : degré de liberté ; ns : différence non significative ; * : différence significative à p< 0,05 ; ** : différence significative à p< 0,01 ; *** : différence significative à p< 0,001.

4.2.1.2.1. Facteurs liés à l'environnement Effet de la ferme

La durée de lactation varie selon la ferme. La durée la plus longue est observée à la ferme F3 qui semble avoir plusieurs problèmes sanitaires conduisant à l'allongement de l'intervalle entre vêlage. La ferme F4 enregistre une durée intermédiaire et comparable à celle des autres fermes.

Tableau 4.32 : Variabilité de la durée de lactation selon la ferme

Ferme	DDL u#177;ó (n) jours		Erreur stand
F1	281#177;53 ^a	(414)	2,60
F2	288#177;5 7 ^a	(309)	3,24
F3	310#177;72 ^b	(93)	7,50
F4	297#177;75 ^ab	(115)	6,97

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; DDL : durée de lactation

Effet de l'année

La variabilité de la durée de production laitière selon l'année (tableau 4.33) montre une forte corrélation avec les variations du niveau de pluies. Les vaches prolongent la durée de lactation durant les années plus pluvieuses (1996, 1997, et 1998). En revanche, la chute des niveaux de précipitation durant les années 2000 et 2001 dont le cumul annuel n'a pas dépassé 350mm/an, semble avoir un effet sur le raccourcissement de la durée de lactation durant cette période.

Tableau 4.33 : Variabilité de la durée de lactation selon l'année de vêlage

Année	DDL piu (n) jours		Erreur Std
94	291#177;65 ^abcd	(12)	18,82
95	289#177;55 ^bc	(94)	5,63
96	287#177;60 ^b	(134)	5,20
97	309#177;64 ^a	(190)	4,63
98	302#177;5 7 ^ab	(167)	4,42
99	274#177;59 ^cd	(176)	4,43
2000	266#177;53 ^d	(101)	5,29
2001	260#177;46 ^d	(36)	7,65

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; DDL : durée de lactation

Effet de l'interaction entre ferme et année

Dans toutes les fermes, la période comprise entre 1995 et 1998 se caractérise par une durée de lactation plus longue que celle réalisée plus tard, confirmant ainsi les résultats observés précédemment. Au delà de 1998, le raccourcissement de la durée de lactation est beaucoup plus observé dans les fermes situées plus au Sud qui semblent plus sensibles au variations climatiques (Figure 4.27).

1993 1994 1995 1996 1997 1998 1999 2000 2001 2002
Année

360

340

320

300

280

260

240

220

200

F1 F2 F3 F4

Figure 4.27 : Variabilité de la durée de lactation selon la fermes et l'année Effet de la saison de vêlage

La variabilité de la durée de lactation selon la saison de vêlage n'est pas significative, bien que les lactations qui débutent en automne et en été sont plus courtes de 10 jours comparées à celles de l'hiver et du printemps (Tableau 4.34).

Tableau 4.34 : Variabilité de la durée de lactation selon la saison de vêlage

Saison de vêlage	DDL piu (n) jours		Erreur Std
Hiver	289#177;5 7	(224)	3,83
Printemps	292#177;62	(258)	3,85
Eté	289#177;62	(193)	4,43
Automne	283#177;60	(235)	3,92

4.2.1.2.2. Facteurs liée à l'animal Effet de la génération animale

La génération animale affecte significativement la durée de lactation. Selon le tableau 4.35, les vaches importées au stade génisse réalisent la durée la plus longue, l'écart est d'un mois de production en plus que les femelles de la génération G4. Pour les générations nées localement, on remarque la supériorité de la G3, qui réalise plus 5 jours par rapport à la G2 et 20 jours comparée à la G4.

Tableau 4.35 : Variabilité de la durée de lactation selon la génération animale
Génération DDL pió (n) jours			Erreur std
G1	309#177;58 ^a	(144)	4,87
G2	285#177;58 ^bc	(302)	3,35
G3	289#177;61 ^b	(341)	3,30
G4	271#177;56 ^c	(71)	6,59

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; DDL : durée de lactation

Effet de l'ordre de parité

Durant la carrière des femelles laitières, la durée de lactation est comparable jusqu'en 5ème parité, bien qu'elle soit légèrement plus longue en première lactation par rapport à la 2^ème et 3ème parité, et se situe à plus de 290 jours. La durée de lactation accuse une légère reprise durant la 4^ème et la 5^ème parité pour se situer à 284 jours (tableau 4.36). En fin de carrière la durée de lactation est significativement plus allongée que durant les 5 premières parité car seules les femelles hautement productrices sont conservés dans les élevages.

Tableau 4.36 : Evolution de la durée de lactation selon la génération animale

Parité	DDL pió (n) jours		Erreur Std
P1	290#177;60	(182)	4,43
P2	283#177;51	(198)	3,61
P3	284#177;68	(153)	5,47
P4	285#177;48	(126)	4,23
P5	285#177;61	(94)	6,33
P6+	300#177;68	(178)	5,11

DDL : durée de lactation

Effet de l'interaction entre génération et parité

Trois formes d'évolution de la durée de lactation au cours de la carrière sont observées (figure
4.28). Chez la génération G1, la durée de lactation augmente en rapport avec l'ordre de parité
et se traduit par un gain en fin de carrière de 31 jours. Pour les générations G2 et G3 la durée

de lactation se maintient au tour de la moyenne avec plus de stabilité pour la génération G3. La génération G4 connaît une régression significative de la durée de lactation à partir de la deuxième mise bas. Chez lez femelles de cette génération la durée de lactation est significativement plus longue au début de carrière (330 jours), elle chute considérablement pour se situer à moins de 240 jours à la 3^ème parité, perdant ainsi environ 3 mois de lactation.

340

320

300

280

260

240

220

200

0 1 2 3 4 5 6 7

Ordre de vêlage

G1 G2 G3 G4

Figure 4.28 : Durée de lactation au cours de la carrière selon la génération animale

Effet de l'âge des femelles

Le tableau 3.37 montre une variabilité de la durée de lactation liée à l'âge de l'animal. En effet, la durée de lactation est plus longue à un âge précoce, diminue jusqu'à l'âge de 7 ans, puis s'allonge en fin de carrière, suite à la conservation des meilleures laitières.

Tableau 4.37 : Variabilité de la durée de lactation selon l'age des femelles

Age	DDL (jours)		Erreur Std
A1	294#177;59 ^ab	(103)	5,85
A2	284#177;56 ^b	(300)	3,23
A3	284#177;59 ^b	(215)	4,00
A4	284#177;59 ^b	(146)	4,90
A5	305#177;70 ^a	(145)	5,81

A1: =3ans; A2: 3-5ans; A3: 5-7ans; A4: 7-9ans; A5: >9ans ; lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; DDL : durée de lactation..

Effet de l'interaction entre génération et âge de la femelle

La variabilité de la durée de lactation selon l'âge des femelles dépend de la génération
animale (Figure 4.29). Si les générations G1 et G3 suivent la tendance globale en prolongeant
les lactation en fin de carrière, Chez la génération G2, l'âge n'a pas d'effet sur la variabilité de

la durée de lactation, alors qu'une évolution différente est observée chez la génération G4. Celle-ci est caractérisée par une lactation plus longue en première parité (340 jours), suit une durée significativement plus réduite à l'âge adulte ; une perte de 100 jours est observée à l'âge de 7ans.

0 1 2 3Age 4 5 6

340

320

300

280

260

240

220

200

360

G 1 G 2 G 3 G 4

Figure 4.29 : Variabilité de la durée de lactation selon la génération et les classes d'âge

4.2.2. Production laitière réelle et standard 4.2.2.1. Performances moyennes

Les vaches Montbéliarde produisent en région semi aride 2741#177;1089kg de lait par lactation, ce qui est l'équivalent de 2794#177;893kg en 305 jours (Tableau 4.38). L'écart type élevé exprime en revanche une variabilité élevée. En effet, 40% des lactations sont inférieurs à 2500kg et 44% se situent entre 2500kg et 4000kg, alors que celles dépassant les 4000kg constituent 16% des lactations. En lactation de 305 jours 60% des vaches produisent moins de 3000kg et seules 9% produisent plus de 4000kg. En Turquie, TILKI et al. (2005) dans une étude faite sur 94 lactations rapportent que chez la race Brune Suisse, 57% des lactations sont inférieurs à 3000kg et seul 13% dépassent 4000kg.

Tableau 4.38 : performances moyennes des paramètres de production du lait

	nombre	min	u#177;ó	max	Coeff. de variation
LC (kg) LS 305j kg)	943 943	519 326	2741#177;1089 2794#177;893	7015 5896	39,74% 31,98%

LC : production par lactation réelle complète ; LS305j : production par lactation standard de 305j a) Production laitière réelle par lactation complète

200

150

100

50

0

25

85

131

1 75

191

142

83

47

23

22

Std. Dev = 1089,05 Mean = 2741

N = 943,00

PL

b) Production laitière de référence à 305 jours

300

266

221

200

158

141

100

Std. Dey = 893,25 Mean = 2794

N = 943,00

54

62

28

0

LS305J

Figure 4.30 : Répartition en classes des paramètres de production du lait

4.2.2.2. Facteurs de variation

L'analyse de la variance a montrée un effet hautement significatif des facteurs de l'environnement ainsi que ceux liés à l'animal sur la production du lait (p<0,001). Selon le tableau 4.3 9, tous les facteurs testés affectent sensiblement le niveau de production.

Tableau 4.39 : Carré moyen obtenu à partir de l'analyse de la variance

Facteurs de variation	DL	PL (kg)	LS 305j (Kg)
Ferme (F)	3	7,4*10^7 ***	6,58*10^7 ***
Année (An)	7	7,2*10^6 ***	3,83*10^6 ***
Saison de vêlage (SV)	3	6,85*10^6 **	5,84*10^6 ***
Génération (G)	3	1,75*10^7 ***	8,90*10^6 ***
Parité (P)	5	4,55*10^6 **	4,30*10^6 ***
Age (A)	4	6,67*10^6 ***	5,78*10^6 ***
F*An	21	3,31*10^6 ***	1,80*10^6 ***
F*SV	9	1,91*10^6 *	1,07*10^6 ns
SV*An	21	1,33*10^6 ns	1,17*10^6 ns
F*G	9	1,93*10^6 *	1,22*10^6 *
G*P	15	2,24*10^6 *	1,92*10^6 **
G*A	12	2,67*10^6 **	2,73*10^6 ***
G*SV	9	1,58*10^6 ns	1,13*10^6 ns
P*SV	15	0,9 7*10^6 ns	0,73*10^6 ns
Erreur	725	7,39*10^5	4,39*10^5

4.2.2.2.1. Effet des facteurs de l'environnement Effet de la ferme

Les vaches élevées dans le semi aride supérieur produisent 800kg de lait en plus par lactation comparée à celles exploitées au semi aride inférieur. Cette différence de production est liée à la durée de lactation et aux disponibilités en fourrages, notamment de meilleure qualité. Les animaux élevés à la ferme F1, située au Sud, enregistrent les résultats les plus faibles. L'effet du stress thermique notamment durant la période estivale et du raccourcissement de la durée de disponibilité en ressources alimentaires vertes durant le pâturage, semblent être la cause d'une telle baisse.

Tableau 4.40 : Variabilité des paramètres de production du lait selon la ferme

Ferme	PL u#177;ó (n) kg		Erreur Std	LS 305j u#177;ó (n) kg		Erreur std
F1	2196#177;816 ^a	(414)	40,09	2307#177;732 ^a	(414)	35,95
F2	3248#177;1042 ^c	(317)	58,51	3344#177;827 ^c	(317)	46,47
F3	3054#177;1094 ^cb	(93)	113,48	2918#177;674 ^b	(93)	69,94
F4	3040#177;1183 ^b	(119)	108,48	2922#177;798 ^b	(119)	73,11

Lettres différentes sur la même colonne exprime une différence significative à p< 0,05 ; PL : production par lactation réelle complète ; PS : production par lactation standard de 305j.

Effet de l'année de vêlage

La production par lactation varie significativement selon l'année (tableau 4.4 1). Cette variabilité est probablement liée aux variations climatiques, au changement de la proportion des femelles importées dans le troupeau, ainsi qu'à la disponibilité des fourragères. Les résultats obtenus manifestent des écarts entre années compris entre 400 à 600kg en faveurs des années plus pluvieuses (1996, 1997 et 1998). Pour l'année 1994, les données recueillies concernent essentiellement les lactations de la ferme F1 située plus au Sud. Cette dernière connaît des difficultés de production liées au stress thermique, à la réduction de la période de disponibilité de fourrages vert et la médiocrité de la qualité des fourrages conservés.

Tableau 4.41 : Variabilité des paramètres de production du lait selon l'année de vêlage

Année	PL u#177;ó (n) kg		Erreur Std	LS 305j u#177;ó (n) kg		Erreur std
94	2154#177;543 ^d	(12)	156,78	2161#177;534 ^d	(12)	154,14
95	2606#177;977 ^cd	(95)	100,19	2629#177;909 ^cd	(95)	93,24
96	2975#177;1132 ^ab	(138)	96,36	3013#177;838 ^a	(138)	71,30
97	3059#177;1212 ^a	(193)	87,21	2950#177;886 ^ab	(193)	63,79
98	2803#177;941 ^bc	(168)	72,60	2794#177;760 ^bc	(168)	58,63
99	2492#177;1051 ^cd	(178)	78,80	2 659#177;932 ^cd	(178)	69,83
2000	2494#177;1044 ^cd	(102)	103,33	2695#177;927 ^c	(102)	91,80
2001	2776#177;1010 ^abcd	(36)	168,39	3132#177;1010 ^a	(36)	168,26

Lettres différentes sur la même colonne exprime une différence significative à p< 0,05 ; PL : production réelle par lactation complète ; LS305j : lactation standard de 305jours

Effet de l'interaction entre ferme et année

La production laitière au cours des années varie selon la ferme (figure 4.31). Les fermes F1 et F2 enregistrent des variations de faible ampleur selon l'année, bien que leurs niveaux de production soient significativement différents du fait qu'elles soient situées dans deux étages bioclimatiques différents. Les deux fermes possèdent plus de prairies naturelles permanentes, respectivement 70 et 45ha, dont le niveau de production pour le pâturage et les stocks fourragers est assez stable quelque soit l'année. Dans les fermes F3 et F4, les niveaux de production sont plus influencés par les variations climatiques et leurs impacts sur les disponibilités fourragères interannuelles.

4000

3500

3000

2500

2000

1500

1000

1992 1994 1996 1998 2000 2002

Année

F1 F2 F3 F4

Figure 4.31 : Variabilité des niveaux de production par lactation selon la ferme et l'année

4000

3500

3000

2500

2000

1500

1000

1993 1994 1995 1996 1997 1998 1999 2000 2001 2002
Année

F 1 F 2 F 3 F 4

Figure 4.32 : Variabilité de la production standard (305j) selon la ferme et l'année

Effet de la saison de vêlage

L'analyse de la variance montre que la production du lait est significativement influencée par la saison de vêlage (tableau 4.42). Les vaches vêlant en hiver produisent en plus 112, 300 et 390 kg de lait par rapport à celles vêlant respectivement au printemps, en automne et en été. Les lactations qui démarrent en hiver et au printemps bénéficient durant la première phase de production, qui coïncide avec la période de pâturage, de ressources fourragères plus importantes et de meilleure qualité.

Les lactations d'automne sont intermédiaires, celle ci bénéficie de la repousse végétale des prairies après les premières pluies d'automne, alors que celles de l'été sont les plus faibles ; le démarrage des lactations coïncide avec la période de stress thermique et une qualité de fourrage plus faible, ce qui se traduit par la baisse significative du niveau de production.

Tableau 4.42 : Variabilité des paramètres de production du lait selon la saison de vêlage

Saison de vêlage	PL pió (n) (kg)		Erreur Std	LS 305j pió (n) kg		Erreur Std
Hiver	2950#177;1086 ^a	(228)	71,91	2998#177;941 ^a	(228)	62,33
Printemps	2838#177;1076 ^ab	(262)	66,48	2862#177;777 ^ab	(262)	47,99
Eté	2559#177;1048 ^c	(195)	75,06	2616#177;912 ^c	(195)	65,33
Automne	2650#177;1102 ^bc	(237)	71,59	2733#177;896 ^bc	(237)	58,20

Lettres différentes sur la même colonne exprime une différence significative à p< 0,05 ; PL : production réelle par lactation complète ; LS305j : lactation standard de 305jours

Effet de l'interaction entre ferme et saison de vêlage

Selon la ferme, trois niveaux d'évolution de la production laitière ont été constatés au cours de l'année (figure 4.33). En effet, dans les fermes F2 et F3, situées au Centre et au Nord, et disposant de prairies de qualité moyenne, l'évolution suit la tendance globale, les meilleurs lactations sont celles qui démarrent en hiver suivies par celles du printemps et d'automne. Dans la ferme F1 située au Sud, l'écart entre saison est très réduit, moins de 4% entre les saisons successives, mais les performances sont les plus faibles. La ferme F3 pratique de l'ensilage et la culture de Sorgho, les lactations d'été sont les plus élevées.

4500

4000

3500

3000

2500

2000

5000

1500

1000

500

0

F1 F2 Ferme F3 F4

Hiver Printemps Eté Automne

Figure 4.33 : Evolution de la production par lactation selon la ferme et la saison de vêlage

4.2.2.2.2. Facteurs liés à l'animal

Effet de la génération animale

Le tableau 4.43 montre des différences significatives des niveaux de production par lactation entre générations. Les vaches importées au stade génisse expriment plus d'aptitude de production laitière et réalisent plus de 3200kg par cycle de production. Les femelles de la génération G2 comparées aux vaches de G1 accusent une chute de production de plus de 700kg traduisant un changement de comportement productif. La génération G3 exprime une reprise de la production en réalisant une augmentation de 250 kg par rapport à la G2, signe d'une réadaptation de ces femelles au nouveau environnement. Une légère chute non significative de la production est observée chez la génération 4 comparée à la G3, moins de 190kg de lait, liée au raccourcissement de la durée de lactation.

Tableau 4.43 : Variabilité des paramètres de production du lait selon la génération animale Génération PL u#177;ó (n) kg Erreur Std LS305j u#177;ó (n) kg Erreur Std

G1	3272#177;1150 ^a	(146)	95,20	3178#177;918 ^a	(146)	75,95
G2	255 7#177;995 ^c	(302)	57,25	2 659#177;922 ^c	(302)	53,06
G3	2791#177;1107 ^b	(349)	59,26	2814#177;849 ^b	(349)	45,42
G4	2603#177;1021 ^bc	(73)	119,44	2779#177;825 ^bc (73)		96,53

Lettres différentes sur la même colonne exprime une différence significative à p< 0,05 ; PL : production réelle par lactation complète ; LS305j : lactation standard de 305jours

Effet de l'interaction entre ferme et génération

Les femelles des générations successives expriment des aptitudes de production laitière différentes selon la ferme (figure 4.34). Dans les fermes F1 et F2, on remarque une nette différence entre la génération importée (G1) et les locales. Les vaches nées localement (G2, G3 et G4) sont comparables entre elles et moins performantes que la G1 ; on enregistre dans ce cas une perte de plus d'un tiers de la production par rapport à la G1. Dans la ferme F3, la production laitière par lactation diminue significativement de l'ordre de 15% par génération. En revanche dans la ferme F4, située plus au Nord, l'effet génération n'est pas significatif et la dégradation de la production n'excédant pas les 4% par génération.

4000

5000

3000

2000

6000

1000

0

F 1 F 2 Ferme F 3 F 4

G 1 G 2 G 3 G 4

Figure 4.34 : Variabilité de la production par lactation selon la ferme et la génération

4000

3000

2000

6000

5000

1000

0

F 1 F 2 Ferme F 3 F 4

G 1 G 2 G 3 G 4

Figure 4.35 : Variabilité de la production par lactation de 305j selon la ferme et la génération

Effet de l'ordre de lactation

Durant leurs carrières, les femelles laitières augmentent constamment leur niveau de production par lactation (Tableau 4.44). Au début de carrière la production est relativement faible (<2600kg), signe d'un format réduit et de besoins de croissance encore élevés. A partir de la 3^ème lactation, la production accuse une augmentation remarquable, qui culmine en fin de carrière, le gain est de 400kg par rapport à la première lactation, sachant que seules les meilleures laitières sont conservées.

Tableau 4.44 : Evolution des paramètres de production du lait au cours de la carrière des femelles

Parité	PL u#177;ó (n)	kg	Erreur Std	LS305j u#177;ó (n) kg		Erreur Std
P1	2562#177;1119 ^a	(184)	82,48	2574#177;8 73 ^a	(184)	64,33
P2	2584#177;990 ^a	(201)	69,85	2664#177; 792 ^a	(201)	55,83
P3	2 744#177;1062 ^ab	(154)	85,55	2861#177;851 ^b	(154)	68,61
P4	2802#177;1132 ^ab	(127)	100,45	2877#177;975 ^b	(127)	86,48
P5	2914#177;1078 ^b	(95)	110,61	299 7#177;852 ^b	(95)	87,36
P6+	2958#177;1118 ^b	(182)	82,84	2937#177;960 ^b	(182)	71,14

Lettres différentes sur la même colonne exprime une différence significative à p< 0,05 ; PL : production réelle par lactation complète ; LS305j : lactation standard de 305jours

Effet de l'interaction entre génération et parité

Au cours de la carrière, la production du lait évolue d'une manière différente d'une génération à une autre (figure 4.36). En effet, chez les femelles importées la production est meilleure au début de carrière, diminue en deuxième lactation, 10% de production en moins, et reprenne en troisième lactation, durant la quelle la production représente le pic de la carrière (supérieur à 3500kg) et correspond à l'évolution de la production des races laitières. Au delà, la production diminue graduellement de 4% par lactation pour se situer à 3100kg en fin de carrière. La génération G2 se caractérise par un niveau de production plus faible jusqu'en 4ème lactation et la production reste inférieure à 2500kg. En revanche son niveau de production augmente significativement en fin de carrière, durant la quelle réalise 16% de production en plus. La production de la génération G3 est comparable à celle de la G2 en premier cycle de production (<2500kg), mais amorce une augmentation dès la 2^ème lactation, pour atteindre en fin de carrière un niveau de production équivalent à celui de la G1 (>3 100). L'amélioration est persistante et graduelle, et passe de 8% au début à 2% en fin de carrière. Chez la génération G4, la première lactation semble épuiser les capacités de production des femelles. A partir de la deuxième lactation, la production accuse une chute de 26%, et se maintien autour de 23 00kg.

Concernant la production à 305j, la figure 4.37 montre une nette différence des niveaux de production entre génération au début de leur carrière : un niveau élevé et supérieur à 3 500kg caractérise les femelles de la G1, un niveau moyen de 3000kg pour les vaches de la G4 et un niveau plus faible (2500kg) pour les laitières des générations G2 et G3. Au cours de la carrière, les niveaux de production par lactation de 305j diminuent chez les femelles de la G1, augmentent chez celles des générations G2 et G3 et se maintiennent autour de la moyenne pour la génération G4. En fin de carrière les niveaux de production sont comparables pour l'ensemble des générations et se rapprochent de la moyenne (3000 kg). En fin de carrière la variabilité de la production du lait entre générations est liée à la durée de lactation, car pour une durée de lactation constante les niveaux de production sont proches et les écarts ne dépassent pas 150kg.

4000

3500

3000

2500

2000

1500

1000

0 1 2 3 4 5 6 7

Numéro de lactation

G 1 G 2 G 3 G 4

Figure 4.36 : Evolution des niveaux de production par lactation selon la génération au cour de la carrière productive des femelles

1000

0 1 2 3 4 5 6 7

Numéro de lactation

3500

3000

2500

2000

1500

4000

G 1 G 2 G 3 G 4

Figure 4.37 : Evolution de la lactation standard selon la génération au cour de la carrière productive des femelles

Effet de l'âge des femelles

La production du lait augmente avec l'âge des femelles et est comparable avec l'évolution selon l'ordre de parité (tableau 4.45). En effet, on note deux niveaux de production : la production du lait est d'environ 2700kg jusqu'à l'âge de sept an, au-delà, la production augmente pour atteindre plus de 3 000kg. La lactation standard, se caractérise par une augmentation de la production jusqu'à l'âge de 9 ans, puis accuse une diminution. En fin de carrière, le niveau de production élevé n'est pas lié à l'augmentation de la production journalière par lactation mais plutôt à l'allongement de la durée de lactation, puisqu'on observe une diminution de la production ramenée à 305j.

Tableau 4.45 : Variabilité des paramètres de production du lait selon l'âge des femelles

Age	PL piu (n) kg		Erreur Std	LS305j piu (n) kg		Erreur Std
A1	2679#177;1158 ^a	(105)	113,05	2670#177;887 ^ab	(105)	86,52
A2	2576#177;1006 ^a	(303)	57,77	2650#177;815 ^a	(303)	46,84
A3	2741#177;1047 ^a	(216)	71,23	2859#177;902 ^b	(216)	61,38
A4	2968#177;1132 ^b	(148)	93,07	3093#177;987 ^c	(148)	81,11
A5	3009#177;1144 ^b	(149)	93,70	2889#177;834 ^b	(149)	68,36

A1: =3ans; A2: 3-5ans; A3: 5- 7ans; A4: 7-9ans; A5: > 9ans ; Lettres différentes sur la même colonne exprime une différence significative à p< 0,05 ; PL : production réelle par lactation complète ; LS305j : lactation standard de 305jours

Effet de l'interaction entre génération et âge de la vache

L'évolution de la production par lactation au cours de la vie productive des femelles est différente selon la génération animale. Les vaches importées au stade génisse accusent une réduction de 15% de leurs productions par lactation durant la période d'âge allant de 3 à 7ans. Au-delà, la production par lactation se maintient autour de la moyenne de cette génération (3300 kg). Les vaches de la génération G2 commencent leurs carrières par des lactations à faible niveau de production (2500 kg). A partir de l'âge de 7ans les femelles améliorent leurs niveaux de production pour atteindre 3000kg à l'âge de 9 ans. Les femelles laitières de la G3 démarrent leurs carrières à un niveau de production comparable à celui des G2 et G4, ensuite on observe une augmentation significative de la production pour atteindre un niveau comparable à la G1 à l'âge de 7 ans. Le niveau de production au début de la carrière de la G4 est supérieur à 3 000kg, il chute significativement et perd plus de 20% de la production par lactation à l'âge de 5 à 7ans. Cet âge constitue en terme de production la limite supérieure de la vie productive des femelles. Elles sont réformées par la suite à cause de la baisse de leur niveau de production.

Pour la génération G4, la chute de production est fortement liée à une réduction de la durée de lactation. La production par lactation de référence à 305j apparaît nettement supérieure aux différentes classes d'âges (figure 4.39) par rapport aux générations G2 et G3.

G 1 G 2 G 3 G 4

1000

0 1 2 3 4 5 6

Classe d'age

4000

3500

3000

2500

2000

1500

Figure 4.38 : Evolution des niveaux de production par lactation selon la génération et l'âge des femelles laitières

1000

3500

3000

2500

2000

1500

4000

G 1 G 2 G 3 G 4

0 1 2 3 4 5 6

Classe d'age

Figure 4.39 : Evolution des lactations à 305j selon la génération et l'âge des femelles laitières

4.2.2.3. Effet de l'âge à la première mise bas sur les paramètres de production du lait D'après le tableau 4.46, nous remarquons que l'âge à la première mise bas influe sensiblement sur les performances post-partum. En effet, la variabilité est plus significative à

partir de la deuxième lactation. En revanche, la durée de lactation est nettement indépendante de l'âge à la première mise bas.

Les performances de production du lait sont significativement supérieures lorsque les vaches vêlent à un âge compris entre 30 et 40 mois. Les vêlages précoces réduisent considérablement le niveau de production, alors que les vêlages tardifs raccourcirent significativement la durée de vie productive. Toutefois, dans notre situation la baisse de l'âge à la première mise bas de 3 6-40 mois à moins de 30 mois conduit à une perte de plus de 1000kg de lait durant les trois première lactation lactations et de 2200kg si les génisses vêlent la première fois à un âge inférieur à 25mois.

Tableau 4.46 : évolution des paramètres de production du lait durant la carrière selon l'age de la 1ère mise bas

APMB	< 25 mois	25-30 mois	30-36mois	36-40 mois	> 40 mois	p
DDL
P1	300#177;46 ^a (8)	292#177;80 ^a (24)	302#177;65 ^a (62)	294#177;60 ^a (35)	304#177;57 ^a (21)	0,950
P2	281#177;68 ^a (7)	2 74#177;50 ^a (23)	286#177;54 ^a (63)	291#177;8 7 ^a (39)	282#177;43 ^a (31)	0,869
P3	229#177;58 ^a (5)	292#177;66 ^ab (56)	294#177;78 ^a (164)	287#177;63 ^ab (88)	290#177;63 ^ab (107)	0,360
PL
P1	2276#177;902 ^a(8)	2346#177;1217^a(24)	2693#177;1132 ^a (62)	2751#177;1118 ^a (35)	2640#177;967 ^a (21)	0,562
P2	2159#177;937 ^ab(7)	2203#177;783 ^b(23)	25 74#177;867 ^ab (63)	2834#177;1325 ^a (39)	2625#177;887 ^ab (31)	0,132
P3+^*	1420#177;570 ^b (5)	2598#177;1178^a (5 6)	2846#177;1140 ^a 164)	2778#177;1011 ^a(88)	2761#177;923 ^a (107)	0,036
LS
P1	2214#177;642 ^a (8)	2353#177;876 ^a(24)	2661#177;849 ^a (24)	2 725#177;852 ^a (35)	2639#177;869 ^a (21)	0,311
P2 *	2166#177;605 ^bc (7)	2348#177;658 ^c(23)	2662#177;723 ^abc(63)	2876#177;990 ^a (39)	2 771#177;784 ^ab (31)	0,037
P3+ *	1 681#177;591 ^c (5)	2613#177;869 ^b(56)	2885#177;839 ^a (164)	2860#177;942 ^ab (88)	2848#177;784 ^ab(107)	0,010

Lettres différentes sur la même ligne expriment une différence significative à p< 0,05 ; APMB : âge à la 1^ère mise bas ; DDL : durée de lactation ; PL : production par lactation ; LS : Lactation standard à 305j.

4.2.3. Durée de tarissement 4.2.3.1. Performance moyenne

La durée de tarissement en situation semi aride est en moyenne de 88 jours avec un coefficient de variation élevé, équivalent à 63%, témoignant de la diversité des pratiques de tarissement. La répartition en classes (figure 4.40) montre que 30% des tarissements durent moins à deux mois, et la majeure partie des femelles sont taries durant un à quatre mois (72%). Notant aussi que certaines femelles sont taries pendant plus de 5 mois (8%).

Tableau 4.47 : Performances moyennes de la durée de tarissement

Nombre Min u#177;ó jours Max Coeff. de variation

DDT (jours) 573 00 88#177;56 424 63,44

DDT : durée de tarissement.

250

200

150

100

50

0

95

166

209

64

19

Std. Dey = 55,87 Mean = 88

N = 573,00

DDT

Figure 4.40 : Répartition des durées de tarissement en classes

4.2.3.2. Facteurs de variation

Les carrés moyens obtenus à partir de l'analyse de la variance montre l'effet hautement significatif de la ferme et de l'année sur la variabilité de la durée de tarissement (tableau 4.48). En revanche, la saison de vêlage et les facteurs liés à l'animal, à l'exception de la génération animale, n'ont pas d'effet significatif sur la variabilité de la durée de tarissement.

Tableau 4.48 : Carrés moyens obtenus à partir de l'analyse de la variance

Facteurs de variation	DL	DDT (jours))
Ferme (F)	3	4,43 ***
Année (An)	5	2,11 **
Saison de vêlage (SV)	3	0,96 ns
Génération (G)	3	1,98 *
Parité (P)	5	0,66 ns
Age (A)	4	1,22 ns
F*An	15	3,02 ***
F*SV	9	0,33 ns
SV*An	15	0,43 ns
F*G	9	0,25 ns
G*P	15	1,23 ns
G*A	12	1,38 *
G*SV	9	0,67 ns
P*SV	15	0,95 ns
Erreur	406	0,61

DDT : durée de tarissement ; DL : degré de liberté ; ns : différence non significative ; * : différence significative à p< 0,05 ; ** : différence significative à p< 0,01 ; *** : différence significative à p< 0,001.

4.2.3.2.1. Facteurs liés à l'environnement

Effet de la ferme

La durée de tarissement est plus courte au Centre et au Nord (F2, F3 et F4) comparés au Sud (F1). En effet, la ferme F2 enregistre la durée la plus courte avec 15 jours de moins par rapport à la ferme F1. Dans les deux autres fermes la durée de tarissement est intermédiaire.

Tableau 4.49 : Variabilité de la durée de tarissement selon la ferme Ferme DDT u#177;ó (jours) Erreur Std

F1	91#177;44 ^a	(242)	2,83
F2	76#177;39 ^b	(216)	2,67
F3	87#177;42 ^ab	(36)	7,04
F4	81#177;5 7 ^ab	(71)	6,78

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; DDT : durée de tarissement

Effet de l'année de vêlage

Une différence significative des durées de tarissement est observée entre années (p<0,05). Les
périodes de tarissement les plus courtes (tableau 4.50) sont enregistrées durant les années les

plus pluvieuses (1996, 1997 et 1998). On observe aussi une relation entre la durée de tarissement et le niveau de production du lait. En effet, les années caractérisées par une production laitière élevée sont celles qui présentent les périodes de tarissement les plus courtes. Durant les années favorables les lactations persistent plus, conduisant à la réduction de la durée de tarissement.

Tableau 4.50 : Variabilité de la durée de tarissement selon l'année

Année	DDT u#177;ó (n) jours		Erreur Std
95	8 7#177;47 ^ab	(69)	5,61
96	85#177;4 7 ^bc	(108)	4,48
97	76#177;45 ^c	(142)	3,76
98	77#177;45 ^bc	(109)	4,31
99	94#177;37 ^a	(88)	3,91
2000	94#177;34 ^a	(40)	5,36

Lettres différentes sur la même colonne expriment une différence significative à p< 0,05 ; DDT : durée de tarissement

Effet de l'interaction année ferme

Dans l'ensemble des fermes, la phase 1996-1998 présente les périodes de tarissement les plus courtes. Selon la figure 4.41, la maîtrise de la durée de tarissement est beaucoup plus élevée dans les fermes F1 et F2.

F 1 F 2 F 3 F 4

150

130

110

90

70

50

30

1993 1994 1995 1996 1997 1998 1999 2000 2001 2002
Année

10

170

Figure 4.41 : évolution de la durée de tarissement selon l'année et la ferme

Effet de la saison de vêlage

Les variations climatiques et des ressources alimentaires entre saisons n'ont pas d'effet sur la
variabilité de la durée de tarissement. Selon le tableau 4.51, les différences ne sont pas

significatives mais les mises-bas de printemps ont une durée de tarissement légèrement plus courte par rapport aux vêlages d'hiver et de l'été (6 et 2 jours respectivement).

Tableau 4.51 : variabilité de la durée de tarissement selon la saison de vêlage

Saison de vêlage	DDT u#177;ó (n) jours	Erreur std
Hiver	8 7#177;43 (156)	3,45
Printemps	81#177;43 (169)	3,29
Eté	84#177;50 (109)	4,81
Automne	85#177;43 (131)	3,77

DDT : durée de tarissement

4.2.3.2.2. Facteurs liée à l'animal Effet de la génération

La génération animale influe significativement sur la durée de tarissement (tableau 4.52). La génération G2 enregistre la durée la plus longue avec plus de 10 jours d'écart par rapport aux autres générations. Celles -ci enregistrent en revanche des durées de tarissement comparables entre elles.

Tableau 4.52	: variabilité de la durée de tarissement selon la génération animale
Génération	DDT u#177;ó (jours)		Erreur Std
G1	79#177;48 ^a	(91)	5,10
G2	90#177;42 ^b	(182)	3,13
G3	80#177;46 ^a	(211)	3,16
G4	81#177;36 ^ab	(43)	5,48

DDT : durée de tarissement

Effet de l'ordre de vêlage

Selon le tableau 4.53, la durée de tarissement est indépendante à la parité, les écarts enregistrés sont de l'ordre de 2 à 8 jours. Toutefois, durant les trois premières lactations, la durée de tarissement est inférieure à 85j ours, augmente légèrement durant la 4^ème , et diminue en fin de carrière pour se situer à moins de 82j ours.

Tableau 4.53 : Variabilité de la durée de tarissement selon l'ordre de lactation

Parité	DDT u#177;ó (jours)		Erreur Std
P1	83#177;40	(135)	3,44
P2	82#177;44	(130)	3,85
P3	85#177;4 7	(99)	4,71
P4	86#177;46	(71)	5,50
P5	90#177;43	(59)	5,61
P6+	81#177;50	(71)	5,92

DDT: durée de tarissement

Effet de l'âge des femelles

L'âge des femelles est sans effet sur la variabilité de la durée de tarissement (tableau 4.54). En effet, une légère réduction non significative de la durée de tarissement est observée entre 3 et 7 ans d'âge, suit une légère augmentation en fin de carrière.

Tableau 4.54 : Variabilité de la durée de tarissement selon l'âge de l'animal
Age DDT u#177;ó (jours)			Erreur Std
A1	89#177;40	(76)	4,54
A2	83#177;44	(208)	3,03
A3	78#177;44	(133)	3,84
A4	92#177;48	(87)	5,12
A5	81#177;47	(61)	6,03

A1: =3ans ; A2: 3-5 ans; A3: 5-7ans; A4: 7-9ans; A5: >9ans ; DDT : durée de tarissement.

4.2.3.3. Durée de Tarissement et production du lait

L'analyse de corrélation entre la durée de lactation précédente et la durée de tarissement, ainsi qu'entre la durée de tarissement d'une part et la durée de lactation suivante, la production réelle et standard de la lactation suivante d'autre part sont illustrées dans le tableau 4.55.

En effet, une corrélation négative est observée entre la durée de lactation précédente et la durée de tarissement, plus la durée de lactation se prolonge l'animal a tendance à réduire sa durée de tarissement et maintenir l'intervalle entre mise bas dans des niveaux moyens. Par contre, la durée de tarissement est sans effet sur la durée et le niveau de production des lactations ultérieures.

Dans nos systèmes de production, la durée de tarissement est peu maîtrisée en terme de période ou de pratiques. En outre, les niveaux de production faible liés à d'autres facteurs limitent l'effet de la préparation au vêlage pendant la période de tarissement.

Tableau 4.55 : Corrélation entre la durée de tarissement et les paramètres de production du lait

	DDLP	DDLS	PPLS	LSS
DDT	-0,20** (572)	0,04 ns (570)	-0,04 ns (569)	-0,07 ns (565)
p	0,000	0,382	0,299	0,104

DDT : durée de tarissement ; DDLP : durée de lactation précédente ; DDDS : durée de lactation suivante ; PPLS : production par lactation suivante ; LSS : lactation standard suivante ; ns : non significatif ; ** hautement significatif ; p : probabilité.

4.3. Etude de la forme des courbes de lactation, des facteurs de variation et prédiction de la production

4.3.1. Périodes critiques

L'évolution de la production quotidienne du lait au cour de l'année exprime trois périodes (Figure 4.42) : (i) une progression de la production quotidienne durant la période Décembre- Mai avec un taux d'accroissement de 5% par mois. Le pic de production est atteint au mois de Mai, et la production moyenne à ce stade est de 11kg par jour, (ii) une régression continue équivalente à 9% par mois est observée de Mai à Août, le niveau de production à la fin de cette phase se situe à 8 kg/jour. (iii) La période allant de Août à Novembre est la plus critique et la production se maintien autour de 8 kg/jour, ce qui ne représente que l'équivalent de 70% de la production au pic.

Cette évolution traduit l'effet important des facteurs de l'environnement sur le niveau de production journalier. La production est plus élevée durant la période favorable allant de fin de l'hiver jusqu'au début de l'été. Lors de cette période le pâturage des prairies contribue pour une grande partie à l'augmentation de la production. La période allant de l'été vers l'automne est moins favorable à la production du lait. Elle se caractérise par des températures élevées et une faible disponibilité en fourrages de qualité.

11,00

10,00

9,00

8,00

Période Critique

7,00

6,00

5,00

4,00

Jan Fév Mar Avr Mai Juin Juit Août Sep Oct Nov Déc
Mois

Figure 4.42 : Période critique de la production du lait durant l'année

4.3.1.1. Périodes critiques selon la ferme

Les différences des disponibilités de ressources fourragères sur le plan qualitatif et quantitatif, sont liées au niveau de l'aridité. Les résultas montrent des variations de forme d'évolution de la production laitière journalière entre fermes. La ferme F1 située plus au Sud se caractérise par des faibles variations de productions entre saisons, le pic du mois de Mai est plus effacé comparativement aux autres fermes. Cette région se caractérise par une durée de la période de vert plus réduite et une qualité des fourrages utilisés plus médiocre. Dans les fermes du Centre et Nord, les variations sont liées aux ressources alimentaires et leurs mode d'utilisation. Dans les fermes du Centre (F2 et F3) la production maximale plafonne durant la période allant de Mars à Mai, toutefois le niveau de production est plus élevé à la ferme F2. Celle-ci dispose de plus de ressources en prairies (45 Ha). La ferme F4 située plus au Nord, se caractérise par une augmentation progressive de la production journalière de lait à partir de Janvier et atteint un pic au début de l'Eté, suit une diminution continue jusqu'à la fin de l'année.

14,00

13,00

12,00

11,00

10,00

4,00

9,00

5,00

8,00

6,00

7,00

Jan Fév Mar Avr Mai Juin Juit Août Sep Oct Nov Déc
Mois

F1 F2 F3 F4

Figure 4.43 : Périodes critiques de la production du lait durant l'année selon la ferme

4.3.1.2. Périodes critiques selon l'année

Deux remarques peuvent être dégagées de l'analyse de la figure (4.44). L'évolution de la production journalière moyenne par vache au cours de l'année montre des différences entre campagnes. Il existe deux pics de production laitière, l'un principale est situé au printemps, durant les mois de Avril et Mai, alors qu'un pic secondaire est observé certaines années uniquement et est situé en mi automne, correspondant à la période de reprise de la végétation dans les prairies. La phase de production minimale se situe en été quelque soit l'année.

13,00

12,00

11,00

10,00

4,00

9,00

8,00

6,00

5,00

7,00

Mois

1995 1996 1997 1998 1999 2000 2001

Figure 4.44 : Périodes critiques de la production du lait durant l'année selon la l'année

4.3.2. Aspect général de la courbe de lactation

A partir du contrôle laitier mensuel, nous avons représenté la courbe générale de lactation pour une durée moyenne de lactation de 10 mois (Figure 4.45). La lactation commence à un niveau de production relativement moyen (10,50 kg/jour), augmente par la suite pour atteindre le pic de production au 2^ème mois post-partum (12 kg/jour) et décline graduellement selon un rythme de baisse mensuel de 5% jusqu'au 6^ème mois de lactation. Du 6^ème au 8ème mois de lactation, la baisse de production est plus de 8% par mois. A la fin de la lactation, la production se maintien autour de 8kg/jour.

0 1 2 3 4 5 6 7 8 9 10 11

Mois post-partum

5,00

4,00

12,00

11,00

10,00

9,00

8,00

7,00

6,00

Figure 4.45 : Evolution de la courbe de lactation

4.3.3. Facteurs de variation de la courbe de lactation

La forme générale de la courbe de lactation varie selon plusieurs facteurs : la ferme, la saison, la génération animale et l'ordre de lactation. Les différences concernent le niveau de production au pic de lactation et la persistance de la production maximale.

4.3.3.1. Effet de la ferme

La figure 4.46 montre des formes d'évolution de la courbe de lactation différentes selon la ferme. Cette variabilité est liée au niveau de l'aridité, qui agit sur la productivité et la qualité des ressources fourragères. Dans la ferme située au Sud (F 1), la production au pic est la moins importante par rapport aux autres fermes. Le déclin est plus rapide, le taux de décroissance se rapproche de 10% par mois. La ferme F2 enregistre une production au pic et une persistance plus élevées (13kg) ; la courbe de lactation se caractérise par la persistance au pic pendant quatre mois, les différences durant cette phase sont inférieurs à 1%. Après le 6^ème mois de lactation la production régresse graduellement à un taux de 11% par mois. Dans les autres fermes on observe l'existence d'un pic de production et un niveau de lactation au pic plus prononcés dans la ferme F4, alors que la ferme F3 se caractérise par une courbe plus plate, une absence de pic de lactation, se traduisant par une évolution de la production plus uniforme au cours de la lactation.

F1 F2 F3 F4

13,00

12,00

11,00

10,00

9,00

8,00

7,00

6,00

5,00

4,00

0 1 2 3 4 5 6 7 8 9 10 11

mois post-partum

Figure 4.46 : Variabilité de l'évolution de la courbe de lactation selon la ferme

4.3.3.2. Effet de la saison

Selon la figure 4.47, on note deux niveaux de production. Un premier relativement élevé pour les lactations débutant durant la période allant de Décembre à Mai. L'analyse des courbes montre que les lactations de Décembre - Janvier se caractérisent par plus de persistance, bien que celles de Février - Mars expriment une production maximale plus élevée au pic de lactation. En revanche, un deuxième niveau de production, plus bas, caractérise les lactations de l'été et de l'automne. Pour les lactations de début d'été (Juin - Juillet) la courbe de lactation est plus plate et uniforme quelque soit la période de mise bas, bien que la production au pic est plus importante.

14,00

13,00

12,00

11,00

10,00

4,00

9,00

8,00

6,00

5,00

7,00

0 1 2 3 4 5 6 7 8 9 10 11

Mois post-partum

Déc/Jan Fév/Mars Avr/Mai Juin/Juit Aout/Sep Oct/Nov

Figure 4.47 : Variabilité de la forme de courbe de lactation selon la saison de vêlage

4.3.3.3. Effet de la génération animale

L'évolution générale de la courbe de lactation est semblable pour l'ensemble des générations, qui réalisent un pic de production durant le 2^ème mois de lactation (Figure 4.48). En revanche, une forte distinction est observée entre vaches importées et celles nées localement concernant le niveau de production au pic. Les femelles de la génération G1 produisent 2kg de lait en plus au pic que les générations nées en Algérie. A partir du 7^ème mois de lactation les niveaux de production de l'ensemble des générations ont tendance à se rapprocher, la différence est de 1 kg/jour.

13,00

12,00

11,00

10,00

9,00

8,00

7,00

6,00

5,00

4,00

14,00

G1 G2 G3 G4

0 1 2 3 4 5 6 7 8 9 10 11

Mois Post-partum

Figure 4.48 : Variabilité de la courbe de lactation selon la génération

4.3.3.4. Effet de l'ordre de lactation

Selon l'ordre de lactation, la figure 4.49 dégage trois niveaux de production. (i) Un niveau de production plus faible est observé durant la première lactation. Le pic de production se maintien à un plafond du 2^ème au 3^ème mois (10,65 et 10,62 kg). La régression par la suite est de 5% par mois jusqu'au 7^ème mois de lactation. (ii) Un niveau moyen est rencontré durant la deuxième parité, la production au pic est supérieure de 1kg de lait par jour par rapport à la première lactation. La production par la suite diminue de 7% par mois, approchant les 10 kg par jour au 5^ème mois. Au delà, la courbe prend la forme d'un escalier en trois plateaux : 5^ème - 6ème _mois, 7^ème -8^ème mois et 9^ème et 1 0^ème mois de lactation. (iii) Les courbes de lactation des mise bas ultérieures sont en revanche superposées jusqu'au 5^ème mois. Une régression est plus rapide uniquement pour la 3^ème lactation entre le 5^ème et le 9^ème mois de lactation.

5,00

4,00

0 1 2 3 4 5 6 7 8 9 10 11

Mois post-partum

13,00

12,00

11,00

10,00

9,00

8,00

6,00

7,00

P1 P2 P3 P4+

Figure 4.49 : Variabilité de la courbe de lactation selon le numéro de lactation

4.3.4. Paramètres de prédiction de la production laitière

La difficulté de suivre les performances de production du lait dans le monde paysan, nous a conduit à chercher des formules pour la prédiction de la production laitière à partir d'un nombre minimum de contrôles. Dans ce chapitre, nous cherchons à préciser les formules selon les facteurs de variation, principalement : l'ordre de la lactation et la saison de son démarrage qui nous semblent des variables pertinents de la production laitière en milieu semi aride.

4.3.4.1. Corrélation entre paramètres

Avant de retenir les éléments de prédiction de la production, nous avons réalisé au préalable une analyse de corrélation entre la production à 305j et la production réelle avec la production journalière de chaque contrôle. Le tableau 4.56 montre des corrélations significatives entres les contrôles rapprochés. En effet, la production à 305j est significativement corrélée (r² >0,70) à la production journalière des cinq mois qui suivent le pic. Par ailleurs, la production totale par lactation est en relation directe non seulement avec la production journalière des différents contrôles mais aussi avec la durée de lactation.

Tableau 4.56 : Corrélation entre la production par lactation et à 305j et les différents contrôles mensuels.

PM1 PM2 PM3 PM4 PM5 PM6 PM7 PM8 PM9 PM10 P305J DDL PPL

PM1 1 ,687** ,531** ,414** ,327** ,240** ,184** ,139** ,114** ,129** ,464** ,162** ,427**

PM2 1 ,773** ,571** ,459** ,372** ,267** ,183** ,207** ,188** ,618** ,158** ,537**

PM3 1 ,742** ,596** ,502** ,345** ,231** ,225** ,206** ,704** ,183** ,612**

PM4 1 ,792** ,666** ,461** ,332** ,239** ,240** ,749** ,177** ,639**

PM5 1 ,789** ,587** ,424** ,301** ,255** ,772** ,155** ,647**

PM6 1 ,760** ,576** ,443** ,359** ,773** ,204** ,680**

PM7 1 ,735** ,593** ,476** ,711** ,242** ,651**

PM8 1 ,751** ,578** ,616** ,302** ,620**

PM9 1 ,761** ,601** ,334** ,620**

PM10 1 ,580** ,353** ,617**

P305J 1 ,330** ,881**

DDL 1 ,698**

PPL 1

PM : production du contrôle au mois 1, 2 ......10 ; DDL : durée de lactation ; PPL : production par lactation réelle complète ; P305j : production par lactation standard de 305j ; ** Corrélation significative à p< 0,01.

4.3.4.2. Prédiction de la production standard à 305j

En utilisant les contrôles laitiers de la période allant du 2^ème au 6^ème mois de lactation, la production standard de 305j peut être prédite à partir d'un seul, deux ou trois contrôles successives ou non. Les formules retenues sont celles qui présentent un coefficient de précision élevé.

Utilisation d'un seul contrôle

A partir d'un seul contrôle la prédiction de la production standard à 305j est faite avec un coefficient de précision relativement limité (Tableau 4.57). Les contrôles utilisés les mieux prédictifs, sont la production au 5^ème ou au 6^ème mois de lactation. Les formules de l'estimation sont les suivantes : P305j=1 159,43+165,26 PC5 (r²=0,61).

P305j=1260,14+165,20 PC5 (r²=0,61).

Tableau 4.57 : Prédiction de la production standard à 305j avec un seul contrôle

variable	R^2	ddl	F	Sig	b1	b0	Err. Std de l'estimation
PC4	0,57	818	1074,19	0,000	157,35	1139,99	571,59
PC5	0,61	816	1289,50	0,000	165,26	1159,43	540,35
PC6	0,61	810	1257,74	0,000	165,20	1260,14	541,70

PCX : production au contrôle du mois X ; r² : coefficient de précision ; ddl : degré de liberté ; b1 et b0 : coefficients de régression.

La figure 4.50 représente la relation entre la production à 305j estimée à partir de la formule citée précédemment et celle calculée réellement à partir des 816 lactations. En effet, avec un coefficient de précision de 61% les points représentant les lactations se répartissent d'une façon plus ou moins homogène autour de la ligne tracée par la formule de prédiction.

4000

5000

3000

2000

6000

7000

1000

0

0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000

Lactation standard estimée

Figure 4.50 : Relation entre les lactations standards estimée avec le model d'un seul contrôle et les lactations réelles (PC5)

Facteurs de variation

La formule de prédiction de la production à 305j est modulée selon l'ordre de parité et la saison de début de lactation. Selon le tableau 4.58, la production est prédite généralement par le biais des contrôles de trois mois quelque soit la saison ou la parité. Ces contrôles sont ceux du 4^ème , 5^ème et 6^ème mois de lactation.

Dans le cas des lactations débutant durant la période comprise entre Octobre et Novembre, un seul contrôle permet d'estimer la production à 305j avec une précision élevée (r²>0,76). Pour ces lactations le 5^ème et le 6^ème moins coïncident avec la période printanière, caractérisée par une production laitière élevée (ressources fourragères disponibles de qualité et températures tolérables) ce qui favorise une meilleure persistance de la production.

Tableau 4.58 : Effet de la parité et de la saison sur les variations des formules de prédiction de la roduction à 305j à partir d'un seul contrôle

Facteur		Cont	r^2	ddl	F	Sig	b1	b0	Err. Std
	Parité1	Pc6	0,69	163	354,07	0,000	185,67	991,32	492,43
		Pc5	0,67	163	338,75	0,000	179,12	944,90	505,13
		Pc4	0,66	163	325,73	0,000	180,17	856,64	511,79
	Parité2	Pc5	0,53	163	182,19	0,000	163,28	1215,25	531,37
		Pc6	0,47	163	145,19	0,000	149,16	1389,35	560,66
		Pc4	0,46	163	140,37	0,000	145,27	1263,64	566,82
	Parité3	Pc5	0,58	120	167,71	0,000	170,10	1099,61	556,72
		Pc6	0,54	120	140,00	0,000	159,10	1381,85	586,12
		Pc4	0,53	120	137,06	0,000	164,79	1075,42	588,97
	Parité4+	Pc6	0,61	361	561,49	0,000	160,21	1329,69	557,85
		Pc5	0,58	362	490,82	0,000	154,77	1289,63	582,46
		Pc4	0,54	364	431,42	0,000	148,29	1257,00	607,15
	Déc-Jan	Pc5	0,64	153	267,92	0,000	162,48	1076,79	575,44
		Pc6	0,61	153	236,69	0,000	159,95	1260,38	598,06
		Pc4	0,59	153	216,41	0,000	165,75	997,51	614,24
	Fév-Mars	Pc4	0,55	121	148,64	0,000	141,98	1251,32	545,89
		Pc5	0,52	121	132,53	0,000	155,39	1311,94	562,97
		Pc3	0,48	121	110,26	0,000	143,41	1113,79	589,46
	Avr-Mai	Pc6	0,57	167	220,33	0,000	179,38	1263,25	528,69
		Pc3	0,54	168	193,00	0,000	166,89	1009,35	556,71
		Pc5	0,51	167	170,98	0,000	178,93	1241,41	565,97
	Juin-Juit	Pc6	0,63	127	217,50	0,000	179,62	1113,43	472,18
		Pc5	0,61	129	197,69	0,000	164,55	1062,98	489,10
		Pc4	0,55	129	156,58	0,000	161,62	1020,51	523,11
	Aôut-Sep	Pc6	0,67	90	180,92	0,000	169,40	1160,90	462,42
		Pc4	0,56	93	119,99	0,000	169,53	1086,44	537,67
		Pc5	0,48	92	84,38	0,000	167,41	1176,90	590,71
	Oct-Nov	Pc5	0,76	144	456,37	0,000	178,66	914,75	468,03
		Pc6	0,74	142	407,92	0,000	173,65	949,91	487,01
		Pc4	0,65	144	269,40	0,000	163,30	1012,18	564,02

PCX : production au contrôle du mois X ; r² : coefficient de précision ; ddl : degré de liberté ; b1 et b0 : coefficients de régression.

Utilisation de deux contrôles

Le coefficient de précision s'élève dans ce cas à plus de 0,73 (tableau 4.59). Pour les trois formules retenues, la connaissance de la production journalière au 6^ème mois de lactation est essentielle pour une meilleure prédiction. La deuxième variable concerne les productions journalières au 2^ème , 3^ème ou 4^ème contrôle. Les deux formules retenues présentent les coefficients de précision (r²) les plus élevés et l'erreur standard la plus faible sont :

P305j = 664,36 + 91,09(PC3) + 1 19,9 1(PC6).

P305j = 582,99 + 83,52(PC2) + 134,46(PC6).

Comparé au graphique précédent, La répartition des productions laitières à 305j estimées par rapport à celles réellement calculées est plus homogène (Figure 4.51).

Tableau 4.59 : Prédiction de la production standard à 305j avec deux contrôles

Variable1	Variable2	r^2	ddl	F	Sig	b1	b2	b0	Err. Std
PC3	PC6	0,73	809	1112,91	0,000	91,09	119,91	664,36	447,14
PC2	PC6	0,73	809	1066,49	0,000	83,52	134,46	582,99	454,14
PC4	PC6	0,70	809	962,84	0,000	86,63	106,35	887,57	471,04

PCX : production au contrôle du mois X ; r² : coefficient de détermination ; ddl : degré de liberté ; b1, b2 et b0 : coefficients de régression.

4500

4000

5500

5000

3500

3000

2500

2000

6500

6000

1500

1000

500

0

0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 5500

Lactation standard estimée (kg)

Figure 4.51 : Relation entre les lactations standards prédites avec le model d'un seul contrôle et réelles (PC2 et PC6)

Facteurs de variation

L'analyse du tableau 4.60 montre que pour toutes les parités et toutes les saisons, la production journalière au contrôle du 6^ème mois de lactation est nécessaire pour la prédiction de la production laitière à 305j. Le deuxième contrôle s'effectue au 2^ème ou au 3^ème mois de lactation selon la parité ou la saison. En revanche, durant la période Août - Novembre, les contrôles du 4^ème ou 5^ème mois sont plus efficaces pour une meilleure précision.

Tableau 4.60 : Effet de la parité et de la saison sur les variations des formules de prédiction de la roduction à 305j à partir de 2contrôles

Facteur		1Cont	2Cont	r^2	ddl	F	Sig	b1	b2	b0	Err. Std
	Parité1	Pc3	Pc6	0,78	163	281,41	0,000	89,25	134,13	521,62	415,80
		Pc4	Pc6	0,77	163	273,52	0,000	97,77	111,47	689,45	42041
		Pc2	Pc6	0,77	163	265,55	0,000	83,68	149,52	440,44	425,22
	Parité2	Pc2	Pc6	0,63	161	136,07	0,000	90,99	127,85	556,69	442,16
		Pc3	Pc6	0,62	161	133,45	0,000	107,72	101,79	658,73	475,04
		Pc2	Pc5	0,62	161	130,93	0,000	72,42	135,53	659,86	479,53
	Parité3	Pc3	Pc6	0,72	118	150,19	0,000	110,64	116,88	498,74	461,16
		Pc2	Pc6	0,71	118	147,59	0,000	103,84	128,63	421,43	464,05
		Pc2	Pc5	0,67	119	123,28	0,000	81,33	133,43	509,41	493,89
	Parité4+	Pc3	Pc6	0,72	360	469,09	0,000	80,81	120,30	775,86	470,25
		Pc2	Pc6	0,71	360	429,96	0,000	73,13	132,94	713,03	485,10
		Pc4	Pc6	0,70	360	414,16	0,000	78,53	107,96	960,47	491,51
	Déc-Jan	Pc3	Pc6	0,71	152	184,05	0,000	96,90	104,43	691,15	517,67
		Pc4	Pc6	0,69	152	167,10	0,000	90,57	96,71	845,41	535,42
		Pc5	Pc6	0,69	152	166,30	0,000	98,65	78,97	967,72	535,29
	Fév-Mars	Pc2	Pc5	0,70	120	138,56	0,000	90,18	134,93	354,93	449,82
		Pc3	Pc6	0,70	120	138,31	0,000	112,65	116,36	460,59	450,11
		Pc2	Pc6	0,70	120	137,29	0,000	108,77	144,74	252,23	451,27
	Avr-Mai	Pc3	Pc6	0,74	166	230,99	0,000	106,10	124,12	568,36	415,20
		Pc2	Pc6	0,71	166	198,31	0,000	89,68	129,43	619,68	438,66
		Pc3	Pc5	0,65	166	155,35	0,000	106,73	109,05	684,56	476,55
	Juin-Juit	Pc2	Pc6	0,76	126	199,72	0,000	96,79	136,92	392,52	382,33
		Pc3	Pc6	0,74	126	174,65	0,000	93,72	138,61	492,77	401,99
		Pc4	Pc6	0,72	126	165,37	0,000	85,96	121,98	730,53	410,08
	Aôut-Sep	Pc4	Pc6	0,73	89	121,44	0,000	76,82	121,00	875,77	417,80
		Pc3	Pc6	0,73	89	118,86	0,000	65,50	142,22	755,50	421,10
		Pc2	Pc6	0,69	89	99,39	0,000	38,17	156,48	879,76	448,67
	Oct-Nov	Pc5	Pc6	0,80	141	275,87	0,000	112,20	73,97	837,81	433,91
		Pc4	Pc6	0,79	141	262,86	0,000	70,32	118,65	775,23	442,30
		Pc3	Pc5	0,79	143	265,09	0,000	45,74	151,32	698,48	441,99

PCX : production au contrôle du mois X ; r² : coefficient de précision ; ddl : degré de liberté ; b1, b2 et b0 : coefficients de régression.

Utilisation de trois contrôles

Selon le tableau 4.61, la différence de précision entre l'utilisation de deux ou trois contrôles n'est pas significative (0,73 vs 0,75). On remarque que l'information du 6^ème mois de lactation est toujours essentielle pour un niveau de prédiction élevée. Les deux formules présentant l'erreur standard les plus faibles sont résumées comme suite :

P305j = 532 + 61,47PC2 + 53,35PC4 + 106,34PC6 (r²>0,75).

P305j = 542 + 69,44PC2 + 61,06PC5 + 91,00PC6 (r²>0,75).

Tableau 4.61 : Prédiction de la production standard à 305j avec trois contrôles

V1	V2	V3	r^2	ddl	F	Sig	b1	b2	b3	B0	Err. Std de l'estimation
PC2	PC4	PC6	0,75	808	822,19	0,000	61,47	53,35	106,34	532,34	430,46
PC2	PC5	PC6	0,75	808	819,81	0,000	69,44	61,06	91,00	542,42	430,93
PC3	PC5	PC6	0,75	808	814,86	0,000	75,32	51,31	86,87	637,43	431,91

PCX : production au contrôle du mois X ; r² : coefficient de détermination ; ddl : degré de liberté ; b1, b2, b3 et b0 : coefficients de régression.

4500

4000

3500

3000

2500

2000

6500

6000

5500

5000

1500

1000

500

0

0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 5500

Lactation standard estimée (kg)

Figure 4.52 : Relation entre les lactations standards prédites avec le model d'un seul contrôle et réelles (PC2, PC4 et PC6)

Facteurs de variation

L'augmentation du nombre de contrôle favorise l'amélioration de la précision de l'estimation de la production. Elle est toujours supérieure durant la première lactation et pour des démarrages de production en été et en automne (r²>0,80). L'analyse du tableau 4.62 confirme la nécessité de connaître la production au 6^ème mois de lactation. Selon la saison et la parité, les autres contrôles utilisés varient du 2^ème au 5^ème mois.

Tableau 4.62 : Effet de la parité et de la saison sur les variations des formules de prédiction de la roduction à 305j à partir de 3contrôles

Facteur		C1	C2	C3	r^2	ddl	F	Sig	b1	b2	b3	b0	Err. Std
	Parité1	Pc3	Pc5	Pc6	0,80	163	216,10	0,000	69,96	63,52	92,44	496,54	393,47
		Pc2	Pc4	Pc6	0,80	163	215,02	0,000	55,87	68,81	109,31	411,04	394,26
		Pc2	Pc5	Pc6	0,80	163	213,00	0,000	66,05	71,07	98,03	414,88	395,75
	Parité2	Pc2	Pc5	Pc6	0,67	160	106,78	0,000	76,27	72,27	77,57	504,13	448,36
		Pc2	Pc3	Pc6	0,66	160	101,52	0,000	55,93	62,21	108,71	455,60	455,92
		Pc2	Pc4	Pc6	0,66	160	101,17	0,000	73,59	49,44	100,02	502,47	456,43
	Parité3	Pc2	Pc4	Pc6	0,74	117	110,31	0,000	78,63	52,81	105,25	361,47	445,68
		Pc2	Pc3	Pc6	0,74	117	110,17	0,000	56,48	64,50	117,92	344,63	445,91
		Pc2	Pc5	Pc6	0,74	117	109,39	0,000	87,93	59,65	87,43	375,22	447,07
	Parité4+	Pc2	Pc4	Pc6	0,74	359	332,55	0,000	51,92	51,88	106,33	647,98	460,01
		Pc3	Pc5	Pc6	0,74	359	331,23	0,000	69,92	39,18	94,63	743,00	460,68
		Pc3	Pc4	Pc6	0,73	359	330,47	0,000	59,71	36,64	106,34	748,20	461,07
	Déc-Jan	Pc3	Pc5	Pc6	0,73	151	133,46	0,000	71,53	55,81	73,16	674,61	502,77
		Pc3	Pc4	Pc6	0,72	151	127,92	0,000	70,49	42,15	90,13	653,16	510,53
		Pc2	Pc3	Pc6	0,71	151	125,10	0,000	26,39	73,80	106,84	631,26	514,61
	Fév- Mars	Pc2	Pc4	Pc6	0,77	119	131,48	0,000	79,07	67,17	105,78	173,34	395,60
		Pc2	Pc5	Pc6	0,75	119	119,03	0,000	96,55	77,67	82,02	144,35	410,82
		Pc2	Pc3	Pc6	0,74	119	112,63	0,000	62,64	65,88	126,51	160,58	419,36
	Avr-Mai	Pc2	Pc3	Pc6	0,75	165	165,27	0,000	40,55	76,27	117,21	472,70	404,76
		Pc3	Pc5	Pc6	0,74	165	158,82	0,000	96,99	33,66	105,79	528,28	410,82
		Pc3	Pc4	Pc6	0,74	165	153,54	0,000	101,62	8,23	121,62	553,94	415,99
	Juin-Juit	Pc2	Pc5	Pc6	0,80	125	163,83	0,000	78,45	61,63	98,70	338,17	352,97
		Pc2	Pc4	Pc6	0,78	125	144,98	0,000	74,01	43,45	118,65	368,70	370,37
		Pc2	Pc3	Pc6	0,77	125	138,21	0,000	71,34	38,04	131,93	330,17	330,17
	Aôut-Sep	Pc3	Pc4	Pc6	0,74	88	85,06	0,000	37,61	49,73	122,45	743,50	410,80
		Pc2	Pc4	Pc6	0,73	88	80,25	0,000	6,52	72,70	121,39	843,03	419,79
		Pc4	Pc5	Pc6	0,73	88	80,05	0,000	76,88	-0,13	121,04	875,98	420,17
	Oct-Nov	Pc3	Pc5	Pc6	0,82	140	206,02	0,000	38,29	95,34	97,22	663,33	414,80
		Pc2	Pc5	Pc6	0,81	140	198,54	0,000	35,77	99,99	70,18	616,50	421,08
		Pc4	Pc5	Pc6	0,81	140	196,73	0,000	41,71	79,56	70,35	766,83	422,64

PCX : production au contrôle du mois X ; r² : coefficient de précision ; ddl : degré de liberté ; b1, b2, b3 et b0 : coefficients de régression.

4.3.4.3. Prédiction de la production réelle par lactation

Contrairement à la production de 305j, l'utilisation des contrôles laitiers mensuels pour la prédiction de la production par lactation ne donne pas des résultats significativement précis (r²<0,60). En revanche, l'introduction dans la formule de prédiction de la variable durée de lactation augmente significativement le coefficient de la précision de la quantité de lait produite par lactation. La précision atteint 78% en cas d'utilisation d'un seul contrôle et dépasse 84% en cas d'utilisation de deux contrôles.

Condition 1 : Sans connaître la durée de lactation Utilisation d'un seul contrôle

Selon le tableau 4.63, un seul contrôle ne permet pas de prédire avec précision acceptable la quantité du lait produite par lactation réelle complète. Néanmoins, le meilleur résultat est obtenu avec l'utilisation du contrôle du 6^ème mois de lactation. Avec un coefficient de détermination r² =0,46, la formule d'estimation est la suivante : PPL = 1073,13 + 182,06PC6

Tableau 4.63 : Prédiction de la production par lactation avec un seul contrôle

variable	r^2	ddl	F	Sig	b1	b0	Err. Std
PC4	0,41	818	563,18	0,000	169,06	982,50	848,12
PC5	0,42	816	587,32	0,000	173,09	1050,65	838,58
PC6	0,46	810	695,40	0,000	182,06	1073,13	802,53

PCX : production au contrôle du mois X ; r² : coefficient de précision ; ddl : degré de liberté ; b1 et b0 : coefficients de régression.

Facteurs de variation

Il est difficile de connaître la production par lactation complète avec un seul contrôle (r²<0,60), bien que la précision diffère d'une parité à une autre et d'une saison à l'autre. Dans la plupart des cas le contrôle au 6^ème mois est le plus fiable pour une meilleure prédiction.

Tableau 4.64 : Effet de la parité et de la saison sur les variations des formules de prédiction de la roduction par lactation à partir d'un seul contrôle

		Cont	r^2	ddl	F	Sig	b1	b0	Err. Std
	Parité1	Pc4	0,53	164	184,69	0,000	209,31	574,70	789,62
		Pc5	0,51	164	173,79	0,000	203,76	717,63	802,22
		Pc6	0,51	164	166,95	0,000	206,84	818,10	798,89
	Parité2	Pc5	0,35	163	86,90	0,000	169,78	1107,04	800,04
		Pc3	0,34	163	84,67	0,000	180,83	748,66	803,63
		Pc6	0,34	163	82,88	0,000	161,33	1234,27	802,67
	Parité3	Pc6	0,34	119	62,43	0,000	158,41	1319,12	873,97
		Pc5	0,31	119	54,24	0,000	155,97	1174,12	897,65
		Pc4	028	119	45,71	0,000	148,72	1178,61	920,46
	Parité4+	Pc6	0,51	361	380,14	0,000	184,04	1076,34	778,79
		Pc5	0,42	362	263,56	0,000	165,62	1161,60	850,57
		Pc4	0,41	364	255,65	0,000	162,03	1084,43	861,75
	Déc-Jan	Pc6	0,58	153	207,65	0,000	167,61	1080,26	669,10
		Pc5	0,53	153	171,24	0,000	159,29	1021,53	705,66
		Pc4	0,50	153	154,18	0,000	165,13	911,17	724,99
	Fév-Mars	Pc3	0,32	121	57,61	0,000	160,67	893,94	913,66
		Pc4	0,32	121	56,62	0,000	147,07	1194,56	916,21
		Pc5	0,29	121	48,50	0,000	156,60	1304,20	937,89
	Avr-Mai	Pc6	0,41	167	116,36	0,000	203,94	1049,46	827,12
		Pc3	0,39	168	106,01	0,000	189,59	761,21	853,32
		Pc5	0,33	167	82,56	0,000	193,63	1115,39	881,34
	Juin-Juit	Pc6	0,45	127	104,88	0,000	194,08	917,92	747,17
		Pc5	0,42	129	92,65	0,000	178,02	880,96	772,93
		Pc4	0,38	129	78,87	0,000	175,01	833,42	798,12
	Aôut-Sep	Pc6	0,54	90	106,52	0,000	224,40	718,53	798,34
		Pc4	0,42	93	65,88	0,000	214,75	693,47	919,18
		Pc5	0,37	92	54,07	0,000	216,57	781,47	954,59
	Oct-Nov	Pc5	0,64	144	251,28	0,000	199,71	664,00	705,05
		Pc6	0,62	142	232,71	0,000	193,32	718,59	717,84
		Pc4	0,57	144	194,24	0,000	187,37	721,69	762,18

PCX : production au contrôle du mois X ; r² : coefficient de précision ; ddl : degré de liberté ; b1 et b0 : coefficients de régression.

Utilisation de deux contrôles

Avec l'utilisation de deux contrôles, la précision s'élève à 0,55. Les contrôles les plus prédictifs sont ceux effectués au 2^ème et 6^ème ou 3^ème et 6^ème mois de lactation. Les formules de prédiction sont :

PPL = 435,76 + 97,45PC3 + 133,60PC6 (r²>0,55).

PPL = 350,15 + 89,18PC2 + 149,14PC6 (r²>0,55).

Tableau 4.65 : Prédiction de la production par lactation avec deux contrôles

Variable1	Variable2	r^2	ddl	F	Sig	b1	b2	b0	Err. Std
PC3	PC6	0,55	809	497,88	0,000	97,45	133,60	435,76	733,09
PC2	PC6	0,55	809	485,26	0,000	89,18	149,24	350,15	738,27
PC4	PC6	0,52	809	441,58	0,000	86,47	123,31	701,22	757,09

PCX : production au contrôle du mois X ; r² : coefficient de détermination ; ddl : degré de liberté ; b1, b2 et b0 : coefficients de régression.

Facteurs de variation

Deux remarques peuvent être dégagées de l'analyse du tableau 4.66. (i) La prédiction de la production nécessite toujours la détermination de la production journalière du 6^ème moins de lactation. (ii) et la prédiction est meilleure (r2=>0,60) pour la première et les dernières lactations ainsi que pour celles qui démarrent pendant l'automne et l'hiver (Septembre - Janvier).

Tableau 4.66 : Effet de la parité et de la saison sur les variations des formules de prédiction de la roduction par lactation à partir de 2 contrôles

Facteur		1Cont	2Cont	r^2	ddl	F	Sig	b1	b2	b0	Err. Std
	Parité1	Pc4	Pc6	0,59	161	117,89	0,000	126,77	110,62	426,66	727,40
		Pc3	Pc6	0,59	161	117,35	0,000	111,83	142,26	229,56	728,39
		Pc3	Pc5	0,58	161	112,54	0,000	105,92	138,48	233,59	748,60
	Parité2	Pc2	Pc6	0,47	161	70,85	0,000	106,15	136,48	262,87	721,95
		Pc3	Pc6	0,45	161	66,59	0,000	119,74	108,68	422,14	732,33
		Pc2	Pc5	0,44	162	62,27	0,000	91,16	134,85	407,94	747,14
	Parité3	Pc2	Pc6	0,46	118	49,80	0,000	104,76	127,68	350,19	798,04
		Pc3	Pc6	0,43	118	45,13	0,000	97,96	121,04	537,25	815,74
		Pc4	Pc6	0,39	118	37,86	0,000	76,60	113,76	893,89	845,78
	Parité4+	Pc3	Pc6	0,60	360	266,11	0,000	86,59	141,27	482,91	709,80
		Pc2	Pc6	0,58	360	252,07	0,000	78,32	154,83	415,90	415,90
		Pc4	Pc6	0,56	360	233,09	0,000	74,90	134,20	724,18	737,62
	Déc-Jan	Pc2	Pc6	0,65	152	143,12	0,000	75,21	136,76	522,86	606,98
		Pc3	Pc6	0,64	152	134,98	0,000	83,30	119,89	590,96	618,58
		Pc4	Pc6	0,63	152	126,46	0,000	76,18	114,42	731,21	631,47
	Fév-Mars	Pc3	Pc6	0,47	120	52,83	0,000	126,63	128,80	170,93	812,83
		Pc2	Pc6	0,46	120	51,62	0,000	120,66	160,85	-43,66	817,25
		Pc2	Pc5	0,46	120	41,83	0,000	103,66	133,08	203,73	855,61
	Avr-Mai	Pc3	Pc6	0,53	166	93,31	0,000	119,71	141,82	265,44	741,46
		Pc2	Pc6	0,48	166	76,73	0,000	85,89	156,10	433,06	778,99
		Pc3	Pc5	0,45	166	67,47	0,000	128,20	109,68	446,51	802,60
	Juin-Juit	Pc3	Pc6	0,54	126	73,73	0,000	111,47	148,87	179,83	688,04
		Pc2	Pc6	0,53	126	71,91	0,000	99,46	153,29	177,17	692,65
		Pc4	Pc6	0,53	126	66,45	0,000	90,50	136,56	514,81	707,12
	Aôut-Sep	Pc3	Pc6	0,58	89	60,94	0,000	74,67	193,41	256,33	770,67
		Pc4	Pc6	0,58	89	60,25	0,000	81,27	173,20	416,89	773,22
		Pc2	Pc6	0,56	89	55,56	0,000	42,63	209,97	404,51	791,13
	Oct-Nov	Pc4	Pc6	0,68	141	147,08	0,000	92,77	120,77	488,16	666,12
		Pc5	Pc6	0,67	141	144,59	0,000	132,09	75,98	586,62	669,96
		Pc3	Pc5	0,66	143	139,89	0,000	54,51	167,12	406,28	681,72

PCX : production au contrôle du mois X ; r² : coefficient de précision ; ddl : degré de liberté ; b1, b2 et b0 : coefficients de régression.

Utilisation de trois contrôles

De l'analyse du tableau 4.67, on observe une faible différence de précision dans le cas d'utilisation de deux ou trois contrôle pour la prédiction de la production par lactation (0,55 vs 0,56). En plus de la production au 6^ème contrôle, la formule de prédiction inclus deux autres variables, la production au 2^ème , 3^ème ou 4^ème contrôle. Les formules utilisées sont :

PPL = 285,97	+ 47,10PC2	+ 61,97PC3	+ 133,91PC6.
PPL = 303,44	+ 68,83PC2	+ 49,21PC4	+ 123,29PC6.

Tableau 4.67 : Prédiction de la production par lactation avec trois contrôles

V1	V2	V3	r^2	ddl	F	Sig	b1	b2	b3	b0	Err.
											Std
PC2	PC3	PC6	0,56	808	347,11	0,000	47,10	61,97	133,91	285,97	724,20
PC2	PC4	PC6	0,56	808	343,36	0,000	68,83	49,21	123,29	303,44	726,42
PC3	PC4	PC6	0,56	808	337,52	0,000	79,04	30,91	121,76	423,23	729,90

PCX : production au contrôle du mois X ; r² : coefficient de précision ; ddl : degré de liberté ; b1, b2, b3 et b0 : coefficients de régression.

Facteurs de variation

Une faible différence de précision de la prédiction est observée entre l'utilisation de deux ou trois contrôles. En revanche, le contrôle du 6^ème mois est toujours obligatoire pour une précision fiable, les deux autres contrôles sont influencés par le numéro de lactation et la saison de démarrage de la lactation (Tableau 4.68).

Tableau 4.68 : Effet de la parité et de la saison sur les variations des formules de prédiction de la roduction par lactation à partir de 3 contrôles

Facteur		C1	C2	C3	r^2	ddl	F	Sig	b1	b2	b3	b0	Err. Std
	Parité1	Pc3	Pc5	Pc6	0,62	160	85,21	0,000	88,18	77,86	91,15	198,82	710,71
		Pc3	Pc4	Pc6	0,61	160	84,58	0,000	67,88	78,65	107,94	218,01	712,34
		Pc2	Pc4	Pc6	0,61	160	83,11	0,000	52,14	99,75	108,60	166,86	716,17
	Parité2	Pc2	Pc3	Pc6	0,48	160	50,12	0,000	71,42	61,62	117,52	162,73	712,96
		Pc2	Pc4	Pc6	0,48	160	49,71	0,000	89,82	46,40	110,36	211,98	714,41
		Pc2	Pc5	Pc6	0,48	160	49,51	0,000	95,10	54,24	98,73	223,43	715,09
	Parité3	Pc2	Pc3	Pc6	0,46	117	33,39	0,000	82,05	30,93	122,54	313,37	798,81
		Pc2	Pc4	Pc6	0,46	117	33,12	0,000	96,85	16,58	120,34	331,36	800,30
		Pc2	Pc5	Pc6	0,46	117	33,10	0,000	99,89	18,26	115,07	336,04	800,43
	Parité4+	Pc2	Pc3	Pc6	0,60	359	181,67	0,000	34,41	61,10	141,03	367,43	705,16
		Pc3	Pc4	Pc6	0,60	359	179,20	0,000	72,91	23,75	132,23	464,98	708,07
		Pc2	Pc4	Pc6	0,60	359	177,53	0,000	60,37	43,91	132,31	360,84	710,05
	Déc-Jan	Pc2	Pc3	Pc6	0,66	151	97,26	0,000	53,88	36,13	124,81	468,71	603,86
		Pc2	Pc4	Pc6	0,66	151	96,80	0,000	61,40	29,56	121,79	489,76	604,81
		Pc2	Pc5	Pc6	0,66	151	96,06	0,000	66,57	24,00	120,60	515,64	606,34
	Fév-Mars	Pc2	Pc3	Pc6	0,49	119	38,76	0,000	67,05	76,57	139,66	-150,18	796,03
		Pc2	Pc4	Pc6	0,49	119	38,12	0,000	95,67	56,52	128,07	-110,04	799,32
		Pc2	Pc5	Pc6	0,47	119	35,66	0,000	113,16	47,66	122,36	-109,87	812,24
	Avr-Mai	Pc2	Pc3	Pc6	0,53	165	62,18	0,000	16,70	107,42	138,89	226,05	742,60
		Pc3	Pc5	Pc6	0,53	165	62,01	0,000	115,77	14,19	134,01	248,55	743,14
		Pc3	Pc4	Pc6	0,53	165	61,91	0,000	124,16	-8,77	144,50	279,78	743,44
	Juin-Juit	Pc3	Pc5	Pc6	0,56	125	53,93	0,000	90,86	65,46	108,65	113,51	671,85
		Pc2	Pc5	Pc6	0,56	125	52,73	0,000	79,67	66,48	112,06	118,54	676,11
		Pc2	Pc3	Pc6	0,55	125	50,88	0,000	52,07	70,83	144,00	61,07	682,81
	Aôut-Sep	Pc3	Pc4	Pc6	0,58	88	41,23	0,000	48,82	46,10	175,09	245,21	769,17
		Pc2	Pc3	Pc6	0,58	88	40,37	0,000	-19,05	89,51	193,70	304,86	773,91
		Pc3	Pc5	Pc6	0,58	88	40,22	0,000	71,84	9,25	188,85	242,65	774,77
	Oct-Nov	Pc4	Pc5	Pc6	0,69	140	104,08	0,000	62,99	82,79	70,50	479,41	653,42
		Pc3	Pc5	Pc6	0,69	140	103,54	0,000	44,34	112,56	68,16	384,58	654,60
		Pc2	Pc5	Pc6	0,69	140	102,74	0,000	47,01	116,05	71,00	295,79	656,34

PCX : production au contrôle du mois X ; r² : coefficient de précision ; ddl : degré de liberté ; b1, b2, b3 et b0 : coefficients de régression.

Condition 2 : Avec la connaissance de la durée de lactation Utilisation d'un seul contrôle

L'introduction de la durée de lactation comme variable dans la formule de prédiction de la production par lactation augmente significativement la précision (tableau 4.69). En effet, pour un seul contrôle la précision est largement plus élevée comparée à celle obtenue sans l'intégration de la durée de lactation (0,78 vs 0,46). Les contrôles de 5^ème et 6^ème mois sont les mieux prédictifs. Les formules utilisées sont :

PPL= -1680,42 + 10,12DDL + 147,87PC5

PPL= -1465,84 + 9,59DDL + 150,66PC6.

Tableau 4.69 : Prédiction de la production par lactation avec un seul contrôle

Variable	r^2	Ddl	F	Sig	b1	DDL	b0	Err. Std
PC4	0,76	818	1297,10	0,000	140,79	9,97	-1652,32	539,67
PC5	0,78	816	1422,59	0,000	147,87	10,12	-1680,42	514,15
PC6	0,78	810	1425,98	0,000	150,66	9,59	-1465,84	514,72

PCX : production au contrôle du mois X ;DDL : durée de lactation ; r² : coefficient de précision ; ddl : degré de liberté ; b1 et b0 : coefficients de régression.

Facteurs de variation

Le tableau 4.70 montre que la précision de l'estimation varie en rapport avec la parité et la

saison. En effet, durant les trois premières lactation la précision est supérieure à 80%, mais au delà la précision se situe autour de la moyenne. La prédiction de la production est meilleure pour les lactations débutant entre le mois Août et le mois de Novembre) utilisant les contrôles du 5^ème et 6^ème mois de lactation qui coïncident avec les périodes favorables de production.

Tableau 4.70 : Effet de la parité et de la saison sur les variations des formules de prédiction de la roduction par lactation à partir d'un seul contrôle et la durée de lactation

Fact eurs		Cont	r^2	ddl	F	Sig	b1	DDL	b0	Err. Std
	Parité1	Pc4	0,82	164	368,83	0,000	155,77	9,69	-1765,70	491,74
		Pc6	0,82	164	356,47	0,000	158,88	10,11	-1771,82	490,13
		Pc5	0,81	164	355,46	0,000	151,71	9,82	-1697,40	499,20
	Parité2	Pc5	0,80	162	330,79	0,000	140,97	11,12	-1870,85	440,70
		Pc6	0,78	161	278,35	0,000	129,81	10,97	-1689,86	468,85
		Pc4	0,77	162	271,57	0,000	124,70	11,00	-1788,41	476,28
	Parité3	Pc5	0,80	119	244,38	0,000	157,98	10,74	-1991,37	480,63
		Pc4	0,77	119	197,34	0,000	150,78	10,74	-1987,07	522,79
		Pc6	0,74	118	170,64	0,000	139,80	9,82	-1385,84	549,31
	Parité4+	Pc6	0,78	360	631,21	0,000	154,88	8,78	-1240,06	526,37
		Pc5	0,75	361	543,63	0,000	145,24	9,67	-1495,70	559,01
		Pc4	0,73	361	483,05	0,000	136,90	9,42	-1423,82	588,41
	Déc-Jan	Pc5	0,79	152	284,26	0,000	151,48	9,09	-1510,21	473,38
		Pc6	0,77	152	253,07	0,000	151,17	7,92	-1022,70	495,31
		Pc4	0,73	152	207,19	0,000	151,74	8,57	-1399,07	533,92
	Fév- Mars	Pc4	0,82	120	269,10	0,000	119,67	11,44	-1903,10	475,99
		Pc5	0,81	120	250,98	0,000	130,55	11,65	-1910,27	489,62
		Pc3	0,77	120	205,56	0,000	117,13	11,04	-1851,10	529,84
	Avr-Mai	Pc6	0,76	166	265,51	0,000	172,69	10,03	-1686,73	527,37
		Pc3	0,74	167	238,76	0,000	155,44	10,02	-1864,52	556,39
		Pc5	0,74	166	236,68	0,000	174,03	10,73	-1934,76	550,64
	Juin-Juit	Pc6	0,81	126	263,52	0,000	166,00	8,80	-1427,86	445,23
		Pc5	0,78	126	223,27	0,000	149,58	8,81	-1425,62	480,08
		Pc4	0,75	128	191,41	0,000	146,49	8,93	-1495,25	509,14
	Août-Sep	Pc6	0,85	89	251,91	0,000	145,21	9,46	-1444,06	459,65
		Pc4	0,85	92	249,87	0,000	142,41	10,27	-1703,43	476,29
		Pc5	0,80	91	183,25	0,000	130,25	10,58	-1632,72	539,39
	Oct-Nov	Pc5	0,84	143	371,94	0,000	162,15	8,89	-1535,64	470,69
		Pc6	0,83	141	332,65	0,000	157,45	8,99	-1535,34	489,35
		Pc4	0,77	143	239,90	0,000	146,63	8,85	-1417,54	561,69

Utilisation de deux contrôles

La précision de l'estimation atteint 84% avec l'utilisation de deux contrôles. Selon le tableau 4.7 1, cette précision est obtenue en utilisant les contrôles du 2^ème ou du 3^ème mois en plus de celui du 6^ème mois de lactation. Les formules retenues sont les suivantes :

PPL (kg) = -1920,72 + 9,25DDL + 83,51PC3 + 110,26PC6.

PPL (kg) = -2008,53 + 9,29DDL + 76,98PC2 + 123,3 3PC6.

Tableau 4.71 : Prédiction de la production par lactation avec deux contrôles

V1	V2	r^2	ddl	F	Sig	b1	b2	DDL	b0	Err. Std
PC3	PC6	0,84	808	1461,76	0,000	83,51	110,26	9,25	-1920,72	432,16
PC2	PC6	0,84	808	1421,68	0,000	76,98	123,33	9,29	-2008,53	437,25
PC4	PC6	0,83	808	1298,42	0,000	78,50	97,83	9,44	-1762,71	454,12

PCX : production au contrôle du mois X ; DDL : durée de lactation ; r² : coefficient de précision ; ddl : degré de liberté ; b1, b2 et b0 : coefficients de régression.

Facteurs de Variation

Quelque soit la parité ou la saison, la précision de l'estimation est élevée (>84%). Elle varie de 83% à 86% selon la parité en faveur de la première et les dernières parités ; la précision est de 84% à 88% selon la saison en faveur de la saison estivale allant du juin à septembre (Tableau 4.72). Les contrôles utilisés sont ceux de 2^ème et 3^ème mois d'une part et du 5^ème et 6ème mois d'autre part.

Tableau 4.72 : Effet de la parité et de la saison sur les variations des formules de prédiction de la roduction par lactation à partir de 2 contrôles et la durée de lactation

Facteurs		Cont1	Cont2	r^2	ddl	F	Sig	b1	b2	DDL	b0	Err. Std
	Parité1	Pc3	Pc6	0,86	160	326,56	0,000	80,99	114,94	9,52	-2045,39	429,20
		Pc4	Pc6	0,86	160	325,52	0,000	90,90	92,80	9,50	-1895,12	429,78
		Pc2	Pc6	0,86	160	321,00	0,000	77,87	126,49	9,85	-2216,70	432,37
	Parité2	Pc2	Pc6	0,83	160	263,36	0,000	71,05	115,40	10,19	-2133,60	407,50
		Pc2	Pc5	0,83	161	261,62	0,000	51,37	122,70	10,57	-2117,87	411,23
		Pc3	Pc6	0,83	160	256,97	0,000	82,32	95,50	10,31	-2073,57	411,68
	Parité3	Pc3	Pc5	0,85	118	216,76	0,000	77,11	144,57	10,67	-2405,74	427,48
		Pc2	Pc5	0,85	118	215,74	0,000	69,14	126,76	10,43	-2400,47	428,33
		Pc3	Pc6	0,85	117	214,36	0,000	105,27	99,30	9,99	-2273,82	426,99
	Parité4+	Pc3	Pc6	0,85	359	651,55	0,000	77,50	117,50	8,51	-1700,44	440,78
		Pc2	Pc6	0,84	359	616,19	0,000	71,89	128,68	8,60	-1798,79	451,25
		Pc4	Pc6	0,82	359	553,85	0,000	69,54	108,99	8,67	-1536,83	471,67
	Déc-Jan	Pc3	Pc6	0,84	151	262,99	0,000	84,33	100,87	8,05	-1570,48	414,45
		Pc2	Pc6	0,83	151	249,30	0,000	67,98	123,88	7,63	-1450,59	423,81
		Pc5	Pc6	0,83	151	239,98	0,000	91,28	75,19	8,43	-1427,60	430,57
	Fév-Mars	Pc2	Pc5	0,87	119	266,63	0,000	74,72	114,91	11,06	-2540,44	402,81
		Pc2	Pc6	0,87	119	259,06	0,000	92,17	123,68	10,50	-2485,95	407,89
		Pc3	Pc6	0,87	119	247,55	0,000	93,00	100,86	10,40	-2257,56	415,98
	Avr-Mai	Pc3	Pc6	0,84	165	297,88	0,000	100,21	122,19	9,54	-2210,70	427,92
		Pc2	Pc6	0,83	165	271,59	0,000	85,90	124,84	10,03	-2303,28	444,81
		Pc2	Pc5	0,81	165	234,41	0,000	87,47	120,01	10,55	-2453,78	472,52
	Juin-Juit	Pc2	Pc6	0,87	125	282,72	0,000	89,02	130,10	8,61	-2039,62	364,72
		Pc3	Pc6	0,85	125	240,89	0,000	81,44	134,39	8,36	-1847,96	390,79
		Pc4	Pc6	0,85	125	231,59	0,000	73,92	119,75	8,57	-1695,14	397,39
	Aôut-Sep	Pc4	Pc6	0,88	88	222,76	0,000	82,03	93,44	9,47	-1750,77	406,96
		Pc3	Pc6	0,88	88	210,16	0,000	65,55	119,00	9,34	-1822,40	417,52
		Pc2	Pc6	0,86	88	180,21	0,000	37,33	132,94	9,41	-1708,86	446,37
	Oct-Nov	Pc3	Pc5	0,86	142	284,72	0,000	46,23	135,14	8,74	-1717,19	444,13
		Pc5	Pc6	0,86	140	275,53	0,000	101,92	68,57	8,58	-1532,69	446,94
		Pc3	Pc6	0,85	140	264,79	0,000	52,09	128,58	8,82	-1752,36	454,59

PCX : production au contrôle du mois X ;DDL : durée de lactation ; r² : coefficient de précision ; ddl : degré de liberté ; b1, b2 et b0 : coefficients de régression.

DISCUSSION

GENERA LE