RICA : Typologie des exploitations en fonction de leur autonomie alimentaire. Cas Concrets de Statistiques 2025
2025-01-12
Introduction
Le jeu de donnée est les résultats comptables de 3985 exploitations issus du Réseau d’Information Comptable Agricole (RICA)
Nous chargeons dans un premier temps le jeu de données simlifié, en ajoutant certaines colonnes (‘OTEXE’, ’ Type d’exploitation’,‘PITOT’ correspondant à l’intraconsommation ) à partir du jeu de données complet de 988 variables. Nous ajoutons aussi une colonne UGB volailles, à partir des colonnes UGB Total et UGB Bovin, Ovins, Caprins et Porcins (soustraction), ainsi que les nouveaux noms de régions.
Présentation générale des données
Nous regardons dans un premier temps le nombre d'exploitations en fonction de(s) l'espèce(s) élevée(s) (nous comptons qu'une exploitation élève un type d'animal à partir de 5 UGB pour cet animal) : nous obtenons le tableau suivant pour les exploitations élevant une ou deux espèces.| Bovins | Ovins | Porcins | Volailles | Caprins | |
|---|---|---|---|---|---|
| Bovins | 2850 | 201 | 144 | 324 | 56 |
| Ovins | 541 | 8 | 54 | 14 | |
| Porcins | 353 | 30 | 0 | ||
| Volailles | 770 | 9 | |||
| Caprins | 145 |
Nous observons donc une plus grande majorité d'élevage bovin qui peuvent être spécialisés ou mixtes. Les élevages spécialisés sont majoritaires. Pour information, pour les nombres exploitations élevant plus de 5
UGB de trois espèces distinctes, les valeurs non nulles sont les
suivantes :
- Ovins, Volailles, Caprins : 2 exploitations
- Bovins, Ovins, Volailles : 31 exploitations
- Bovins, Volailles, Caprins : 3 exploitations
- Bovins, Ovins, Caprins : 6 exploitations
- Bovins, Ovins, Porcins : 4 exploitations
- Bovins, Ovins, Volailles, Caprins : 1
exploitation
La répartition des exploitations en fonction du nombre d’UGB Total
sur les exploitations donne le graphe suivant:
Nous observons deux choses : d'abord, que la variable UGBTO n'est pas une variable réellement quantitative, mais une variable qualitative avec certains choix de classes. Par exemple, les valeurs de UGBTO=400 correspondent à la moyenne des valeurs d'UGBTO comprises entre 375 et 425. Pour régler ce problème, nous créons des classes d'UGB qui pourront servir en tant que variables qualitatives lors de futures ACP. Aussi, certaines des exploitations possèdent moins de 10 UGB. Nous décidons de les exclure de notre analyse (pourquoi?compléter)
On peut décomposer la répartition des exploitations en fonction du nombre d’animaux qu’ils possèdent en fonction des différentes espèces : on voit alors que la répartition globale est fortement impactée par les exploitations bovines, majoritaires :
Le nombre d’UGB par exploitation se trouve souvent entre 50 et 150 UGB. Ceci est dû à l’influence des exploitations comportant des bovins (71% des exploitations) et au fait qu’elles soient souvent spécialisées en bovins. On remarque également une forte densité d’exploitations ayant plus de 500 UGB. Ceci est lié aux exploitations porcines (et dans une moindre mesure avicoles et bovines) qui comportent souvent plus de 500 UGB.
On aensuite regardé le nombre moyen d’UGB (total) pour chaque type d’exploitation (*moyenne pour chaque type d’exploitation représentée sur le graphe en rouge).Pour rappel les rectangles des boxplot sont délimités par les premier et troisième quartile:
On observe ici aussi que les exploitations à dominante élevage de petits ruminants possèdent peu d’UGB (en moyenne un peu moins de 100). Parmi les exploitations bovines, les exploitations mixtes (lait+viande) se distinguent des exploitations spécialisées par un nombre d’UGB plus important : environ 210 contre 130 UGB. Les exploitations de polyélevage sont en moyenne plus grandes (plus d’UGB), et ce, encore plus si elles sont à orientation granivores. De même, les exploitations à dominante élevage de granivore possèdent en moyenne un nombre élevé d’UGB (310 UGB). Les exploitations mixtes grandes cultures et herbivores sont à un niveau d’UGB intermédiaire entre les exploitations bovines et de petits ruminants. Les exploitations ayant des ateliers de cultures majoritaires (viticultures, fruitières, céréales…) ont peu d’UGB (une cinquantaine en moyenne).
On s’est ensuite intéressés aux régions pour voir si le nombre d’UGB leur était lié, pour les UGB Totaux et les différentes espèces (différents UGB):
Les régions situées dans le grand ouest (Bretagne et Pays de la Loire) ont des exploitations qui comportent le plus d’UGB ce qui est cohérent car ces régions ont la plus grande concentration en élevages bovins et porcins ainsi que de nombreux élevages de volailles. Les élevages bovins de l’ouest sont souvent des élevages laitiers et ont donc un nombre d’UGB élevés contrairement aux élevages bovins de régions montagneuses ou du centre France qui ont des effectifs plus restreints. Nous avons également vu plus tôt que les élevages porcins et avicoles ont des nombres d’UGB supérieurs à 500, de ce fait les régions comportant ces types d’élevage ont un nombre d’UGB élevé.
Pour les élevages avicoles, la Bretagne comporte les élevages ayant le plus d’UGB (de 90 à plus de 500). Les autres régions ont des élevages de tailles plus restreintes : de taille moyenne (jusqu’à 275 UGB environ) pour les régions Normandie, Pays-de-la-Loire et Hauts-de-France et de petites tailles pour les autres régions (moins de 200 UGB).
Comme précisé précédemment, les régions comportant des élevages porcins et avicoles ont un nombre d’UGB très élevé (elles comportent entre 100 et plus de 500 UGB). Hormis dans les régions Auvergne-Rhône-Alpes et Hauts-de-France où les exploitations ont entre 100 et 375 UGB. Ainsi, le nombre d’UGB porcins reste assez uniforme sur toute la France, indépendamment des régions.
Les UGB ovins ont une répartition différente selon les régions. Les plus grandes exploitations se trouvent dans les régions PACA et Grand Est. Ce ne sont pourtant pas les régions ayant le plus d’élevages qui sont la Nouvelle Aquitaine et l’Occitanie. Les autres régions ont des exploitations de moins de 100 UGB. La Bretagne est un cas particulier avec une exploitation de plus de 100 UGB.
Les exploitations bovines sont de taille plutôt uniforme sur toute la France (entre 50 et 200 UGB environ).
Les élevages caprins sont de tailles très hétérogènes selon les régions. Certaines régions n’ont que très peu d’exploitations caprines. Celles-ci sont principalement en Nouvelle-Aquitaines, Centre-Val de Loire et Pays-de-la-Loire. La Bretagne a peu d’exploitations caprines mais de tailles très élevées (une centaine d’UGB) alors que les autres régions ayant des élevages caprins ont des effectifs entre 20 et 50 UGB.
Aussi, en observant la colonne SAUTI, nous avons observé que le problème de “fausses classes” était le même que pour les UGB. Nous avons donc créé une colonne de SAUTI par classes.
Analyse préliminaire des OTEXE
On observe:
-Les charges totales fossiles sont fortement corrélées à SUT1GCFOU -LANMX : le loyer des animaux n’est corrélé à aucune autre variable. -Les achats d’aliments grossiers (ACHAG) ne sont pas corrélés à d’autres variables que celles qu’elles composent (Charges d’aliments grossiers car peu de variations de stock) -Toutes les variables du bilan comptable sont fortement corrélées entre elles (CHRTO,PBRTO,VAPBR,EBEXP,RESEX,CHSOX). L’intraconsommation de produits végétaux pour les animaux (Intraconsommation_animaux) est corrélée aux éléments comptables ((CHRTO,PBRTO,VAPBR,EBEXP,RESEX), car elle permet une réduction des charges tout en permettant une production de produits animaux. Elle est décorrélée des surface.
On trace la même ACP avec les variables bien représentées en noires et celles moins bien représentées sont en grise et seront retirées pour la prochaine ACP
Nous avons réalisé une nouvelle ACP sans le type d’exploitation “granivore” pour voir si elle influençait les résultats de l’ACP. Les exploitations qui se distinguent le plus sont les exploitations bovines ( viande et mixte) et petits ruminants possèdant plus de SFP par rapport aux autres exploitations. Les exploitations bovines laitières ne sont pas bien représentées. Les exploitations de poly élevage ont plus de charges alimentaires, ce qui est cohérent avec l’ACP précédente.
Etude des Achats alimentaires
On s’intéresse désormais aux achats alimentaires. Les charges alimentaires (charges d’aliments concentrés ‘CHRAC’ et grossiers ‘CHRAG’) sont liées à l’alimentation du troupeau. Nous avons créé une variable ‘ACHA_SOMME’ de la somme de ces deux charges.
Dans un premier temps, on peut comparer les Achats alimentaires pour différents types d’exploitation : en premier, la somme, puis les deux types d’aliments (concentrés et grossiers) :
En dehors du fait que l’échelle est écrasée par les valeurs très élevées de charges des élevages granivores,on observe que les charges dépendent fortement du type d’exploitation, et ce, en grande partie parce que les exploitations ont des nombre d’animaux à nourrir très différents. Pour réellement savoir si les exploitations reposent sur ces charges pour nourrir leurs animaux, nous créons des variables de charge par unité d’UGB (on divise ces charges par le nombre total d’UGB). En effet, on voit bien que les charges augmentent en fonction du nombre total d’UGB, ce qui est logique (et on a vu précédemment que les différents types d’exploitation avaient des nombre d’UGB distincts).
On trace les mêmes graphes de charges par type d’exploitation, mais désormais les charges sont divisées par le nombre d’UGB total pour chaque exploitation :
Les exploitations de granivores achètent davantage de concentrées par
rapport aux autres types d’exploitation. En bovin, ce sont les
exploitations spécialisées en lait qui achètent le plus de
concentré.
Nous avons des exploitations très hétérogènes certaines n’achètent aucun concentré alors que d’autres en consomment énormément.
Pour construire un modèle ayant comme variable à expliquer les Achats alimentaires, nous regardons dans un premier temps les relations (analyse bivariée) entre les Achats alimentaires et les autres variables :
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
Pour pouvoir proposer un modèle linéaire, nous avons tout d’abord réaliser des analyses bivariées entre notre variable d’intérêt “ACHA_SOMME” qui correspond à la somme des achats extérieurs d’aliment de l’exploitation en fonction de différentes variables : Charges totales ( hors charges sociales) (CHRGTO), l’EBE, l’intraconsommation des animaux (intra_conso_animaux), les charges foncières (LFERM+LMATE), le produit total brut (PBRTO), la production brute standard(PBUCE), le résultat de l’exploitation (RESEX), la SAU (SAUTI), la surface fourragère principale ( SFP), l’UGB total (UGBTO) et la valeur ajoutée (VAPBR).
Les variables UGBTO et SAUTI ne sont pas continues mais l’achat alimentaire augmente avec l’UGB de façon exponentielle.
Entre la SFP et l’achat alimentaire il semble y avoir une relation d’exponentielle décroissante nous aurions peut être besoin de modifier la variable SFP en exp(-SFP) pour proposer un modèle ultérieurement.
Il ne semble pas avoir de relation avec les charges foncières.
Pour les variables économiques nous avons différentes relations : - pour la charge totale, le produit total nous avons une relation linéaire. - pour la valeur ajoutée, l’EBE et le résultat la relation semble quadratique.
La relation entre intra consommation et l’achat alimentaire semble quadratique.
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
Les relations entre la variable d’intérêt ratio_achat et les autres variables n’ont pas de relations entre elles. Les valeurs des x sont très regroupées et il est diffcile d’établir des relations.
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
Idem que précédemment, on n’a pas de relation entre les variables
explicatives et la variable d’intérêt.
Pour mieux identifier les types d’exploitations, nous avons regrouper les OTEXE en 4 catégories : Polyculteurs, Petits ruminants, Bovins et Granivores.
Ce regroupement nous permet de mieux différencier les catérogies entre elles. Les granivores se distinguent toujours bien des autres exploitations herbivores et polyculteurs. Les exploitations de granivores se caractérisent par des achats d’aliments plus importants, une plus importante consommation d’électricité et des charges totales plus élevées.
Alors que les exploitations herbivores sont caractérisées par une SFP plus importantes et davantage de charges fossiles liées notamment à la surface. Plus il y a de surface, plus il y a des passages de tracteurs notamment.
Nous allons essayer d’analyser les grandes classes séparemment (la construction d’un modèle linéaire global ne respentent pas l’hypothèse d’homoscédasticité).
Modèle linéaire sur la grande classe Bovin :
##
## Call:
## lm(formula = ACHA_SOMME ~ OTEXE + REGIONS + SFP + SUT1GCU + VAPBR +
## EBEXP + CHRPV + Total_charges_fossiles + ELECT + PITOT +
## INTVP + Intra_conso_animaux + UGBTO_classe + SAUTI_classe,
## data = data_mod_Bovins)
##
## Residuals:
## Min 1Q Median 3Q Max
## -177625 -13510 -843 10906 147787
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.406e+04 3.652e+03 3.849 0.000123 ***
## OTEXE46 -1.590e+04 1.966e+03 -8.089 1.07e-15 ***
## OTEXE47 -1.059e+04 2.409e+03 -4.397 1.16e-05 ***
## OTEXE83 -1.961e+04 2.500e+03 -7.846 7.13e-15 ***
## REGIONSBourgogne-Franche-Comté -1.357e+03 2.348e+03 -0.578 0.563501
## REGIONSBretagne -1.181e+04 3.026e+03 -3.902 9.89e-05 ***
## REGIONSCentre-Val de Loire -6.375e+03 3.660e+03 -1.742 0.081670 .
## REGIONSGrand Est 4.124e+03 2.435e+03 1.694 0.090509 .
## REGIONSHauts-de-France 3.482e+03 2.948e+03 1.181 0.237632
## REGIONSIle-de-France 1.395e+03 1.130e+04 0.123 0.901805
## REGIONSNormandie -1.368e+03 2.787e+03 -0.491 0.623566
## REGIONSNouvelle-Aquitaine 4.238e+03 2.330e+03 1.819 0.069121 .
## REGIONSOccitanie 2.778e+03 2.777e+03 1.001 0.317188
## REGIONSPACA -9.355e+03 7.528e+03 -1.243 0.214109
## REGIONSPays de la Loire -4.903e+03 2.809e+03 -1.746 0.081016 .
## SFP -1.432e+02 2.002e+01 -7.153 1.21e-12 ***
## SUT1GCU -2.256e+02 4.573e+01 -4.933 8.81e-07 ***
## VAPBR -7.667e-03 2.879e-02 -0.266 0.790036
## EBEXP 9.845e-02 2.896e-02 3.400 0.000688 ***
## CHRPV 2.587e+00 1.610e-01 16.071 < 2e-16 ***
## Total_charges_fossiles 9.038e-01 9.751e-02 9.269 < 2e-16 ***
## ELECT 3.337e+00 2.292e-01 14.561 < 2e-16 ***
## PITOT 2.368e-01 2.131e-01 1.111 0.266696
## INTVP 4.829e+00 1.328e+01 0.364 0.716101
## Intra_conso_animaux -7.469e-01 2.232e-01 -3.347 0.000833 ***
## UGBTO_classe50-99 3.000e+03 2.609e+03 1.150 0.250447
## UGBTO_classe100-149 1.016e+04 3.063e+03 3.317 0.000927 ***
## UGBTO_classe150-199 1.914e+04 3.711e+03 5.157 2.77e-07 ***
## UGBTO_classe200-249 3.061e+04 4.426e+03 6.915 6.39e-12 ***
## UGBTO_classe250-299 3.088e+04 5.145e+03 6.002 2.33e-09 ***
## UGBTO_classe300-349 4.009e+04 6.587e+03 6.086 1.40e-09 ***
## UGBTO_classe350-399 7.832e+04 8.030e+03 9.753 < 2e-16 ***
## SAUTI_classe50-99 7.383e+03 3.008e+03 2.455 0.014185 *
## SAUTI_classe100-149 1.371e+04 4.083e+03 3.357 0.000805 ***
## SAUTI_classe150-199 2.143e+04 5.559e+03 3.855 0.000120 ***
## SAUTI_classe200-249 2.254e+04 7.233e+03 3.116 0.001862 **
## SAUTI_classe250-299 3.951e+04 9.322e+03 4.238 2.36e-05 ***
## SAUTI_classe300-349 3.845e+04 1.094e+04 3.514 0.000453 ***
## SAUTI_classe350-399 4.176e+04 1.415e+04 2.951 0.003212 **
## SAUTI_classe400-449 3.174e+04 1.773e+04 1.790 0.073596 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 26890 on 1890 degrees of freedom
## Multiple R-squared: 0.7586, Adjusted R-squared: 0.7537
## F-statistic: 152.3 on 39 and 1890 DF, p-value: < 2.2e-16## Start: AIC=39410.19
## ACHA_SOMME ~ OTEXE + REGIONS + SFP + SUT1GCU + VAPBR + EBEXP +
## CHRPV + Total_charges_fossiles + ELECT + PITOT + INTVP +
## Intra_conso_animaux + UGBTO_classe + SAUTI_classe## 2 3 4 5 6 19 33 36 43 70 111 120 139 159 180 183
## 2 3 4 5 6 19 33 36 43 70 111 120 139 159 180 183
## 193 199 201 207 234 255 267 270 271 275 284 289 300 312 315 342
## 193 199 201 207 234 255 267 270 271 275 284 289 300 312 315 342
## 363 373 380 401 409 411 415 432 441 454 457 463 470 484 485 531
## 363 373 380 401 409 411 415 432 441 454 457 463 470 484 485 531
## 534 547 556 557 588 618 626 627 636 675 676 687 688 697 711 716
## 534 547 556 557 588 618 626 627 636 675 676 687 688 697 711 716
## 723 734 766 772 773 781 790 792 800 804 806 813 818 830 835 879
## 723 734 766 772 773 781 790 792 800 804 806 813 818 830 835 879
## 881 886 891 897 903 905 908 914 919 922 929 943 956 968 978 979
## 881 886 891 897 903 905 908 914 919 922 929 943 956 968 978 979
## 986 1005 1006 1009 1019 1022 1028 1034 1040 1078 1097 1100 1103 1116 1121 1146
## 986 1005 1006 1009 1019 1022 1028 1034 1040 1078 1097 1100 1103 1116 1121 1146
## 1166 1177 1181 1187 1193 1213 1241 1250 1251 1271 1304 1349 1360 1372 1382 1411
## 1166 1177 1181 1187 1193 1213 1241 1250 1251 1271 1304 1349 1360 1372 1382 1411
## 1425 1430 1437 1440 1501 1515 1529 1548 1577 1583 1647 1683 1696 1700 1710 1715
## 1425 1430 1437 1440 1501 1515 1529 1548 1577 1583 1647 1683 1696 1700 1710 1715
## 1722 1734 1775 1782 1786 1803 1807 1844 1897 1921 1929
## 1722 1734 1775 1782 1786 1803 1807 1844 1897 1921 1929##
## Call:
## lm(formula = ACHA_SOMME ~ OTEXE + REGIONS + SFP + SUT1GCU + VAPBR +
## EBEXP + CHRPV + Total_charges_fossiles + ELECT + PITOT +
## INTVP + Intra_conso_animaux + UGBTO_classe + SAUTI_classe,
## data = data_mod_Bovins)
##
## Residuals:
## Min 1Q Median 3Q Max
## -177625 -13510 -843 10906 147787
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.406e+04 3.652e+03 3.849 0.000123 ***
## OTEXE46 -1.590e+04 1.966e+03 -8.089 1.07e-15 ***
## OTEXE47 -1.059e+04 2.409e+03 -4.397 1.16e-05 ***
## OTEXE83 -1.961e+04 2.500e+03 -7.846 7.13e-15 ***
## REGIONSBourgogne-Franche-Comté -1.357e+03 2.348e+03 -0.578 0.563501
## REGIONSBretagne -1.181e+04 3.026e+03 -3.902 9.89e-05 ***
## REGIONSCentre-Val de Loire -6.375e+03 3.660e+03 -1.742 0.081670 .
## REGIONSGrand Est 4.124e+03 2.435e+03 1.694 0.090509 .
## REGIONSHauts-de-France 3.482e+03 2.948e+03 1.181 0.237632
## REGIONSIle-de-France 1.395e+03 1.130e+04 0.123 0.901805
## REGIONSNormandie -1.368e+03 2.787e+03 -0.491 0.623566
## REGIONSNouvelle-Aquitaine 4.238e+03 2.330e+03 1.819 0.069121 .
## REGIONSOccitanie 2.778e+03 2.777e+03 1.001 0.317188
## REGIONSPACA -9.355e+03 7.528e+03 -1.243 0.214109
## REGIONSPays de la Loire -4.903e+03 2.809e+03 -1.746 0.081016 .
## SFP -1.432e+02 2.002e+01 -7.153 1.21e-12 ***
## SUT1GCU -2.256e+02 4.573e+01 -4.933 8.81e-07 ***
## VAPBR -7.667e-03 2.879e-02 -0.266 0.790036
## EBEXP 9.845e-02 2.896e-02 3.400 0.000688 ***
## CHRPV 2.587e+00 1.610e-01 16.071 < 2e-16 ***
## Total_charges_fossiles 9.038e-01 9.751e-02 9.269 < 2e-16 ***
## ELECT 3.337e+00 2.292e-01 14.561 < 2e-16 ***
## PITOT 2.368e-01 2.131e-01 1.111 0.266696
## INTVP 4.829e+00 1.328e+01 0.364 0.716101
## Intra_conso_animaux -7.469e-01 2.232e-01 -3.347 0.000833 ***
## UGBTO_classe50-99 3.000e+03 2.609e+03 1.150 0.250447
## UGBTO_classe100-149 1.016e+04 3.063e+03 3.317 0.000927 ***
## UGBTO_classe150-199 1.914e+04 3.711e+03 5.157 2.77e-07 ***
## UGBTO_classe200-249 3.061e+04 4.426e+03 6.915 6.39e-12 ***
## UGBTO_classe250-299 3.088e+04 5.145e+03 6.002 2.33e-09 ***
## UGBTO_classe300-349 4.009e+04 6.587e+03 6.086 1.40e-09 ***
## UGBTO_classe350-399 7.832e+04 8.030e+03 9.753 < 2e-16 ***
## SAUTI_classe50-99 7.383e+03 3.008e+03 2.455 0.014185 *
## SAUTI_classe100-149 1.371e+04 4.083e+03 3.357 0.000805 ***
## SAUTI_classe150-199 2.143e+04 5.559e+03 3.855 0.000120 ***
## SAUTI_classe200-249 2.254e+04 7.233e+03 3.116 0.001862 **
## SAUTI_classe250-299 3.951e+04 9.322e+03 4.238 2.36e-05 ***
## SAUTI_classe300-349 3.845e+04 1.094e+04 3.514 0.000453 ***
## SAUTI_classe350-399 4.176e+04 1.415e+04 2.951 0.003212 **
## SAUTI_classe400-449 3.174e+04 1.773e+04 1.790 0.073596 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 26890 on 1890 degrees of freedom
## Multiple R-squared: 0.7586, Adjusted R-squared: 0.7537
## F-statistic: 152.3 on 39 and 1890 DF, p-value: < 2.2e-16## Anova Table (Type II tests)
##
## Response: ACHA_SOMME
## Sum Sq Df F value Pr(>F)
## OTEXE 6.9980e+10 3 32.2522 < 2.2e-16 ***
## REGIONS 3.5185e+10 11 4.4225 1.344e-06 ***
## SFP 3.7008e+10 1 51.1685 1.207e-12 ***
## SUT1GCU 1.7600e+10 1 24.3347 8.806e-07 ***
## VAPBR 5.1291e+07 1 0.0709 0.7900361
## EBEXP 8.3614e+09 1 11.5606 0.0006876 ***
## CHRPV 1.8681e+11 1 258.2854 < 2.2e-16 ***
## Total_charges_fossiles 6.2132e+10 1 85.9059 < 2.2e-16 ***
## ELECT 1.5335e+11 1 212.0273 < 2.2e-16 ***
## PITOT 8.9279e+08 1 1.2344 0.2666960
## INTVP 9.5687e+07 1 0.1323 0.7161008
## Intra_conso_animaux 8.1016e+09 1 11.2015 0.0008333 ***
## UGBTO_classe 8.7551e+10 7 17.2929 < 2.2e-16 ***
## SAUTI_classe 2.1676e+10 8 3.7463 0.0002309 ***
## Residuals 1.3670e+12 1890
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Call:
## lm(formula = ACHA_SOMME ~ OTEXE + REGIONS + SFP + SUT1GCU + EBEXP +
## CHRPV + Total_charges_fossiles + ELECT + Intra_conso_animaux +
## UGBTO_classe + SAUTI_classe, data = data_mod_Bovins)
##
## Residuals:
## Min 1Q Median 3Q Max
## -178091 -13515 -814 10855 149401
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.431e+04 3.608e+03 3.966 7.58e-05 ***
## OTEXE46 -1.614e+04 1.951e+03 -8.274 2.42e-16 ***
## OTEXE47 -1.072e+04 2.402e+03 -4.461 8.63e-06 ***
## OTEXE83 -1.974e+04 2.489e+03 -7.932 3.65e-15 ***
## REGIONSBourgogne-Franche-Comté -1.401e+03 2.330e+03 -0.601 0.547672
## REGIONSBretagne -1.178e+04 2.964e+03 -3.976 7.27e-05 ***
## REGIONSCentre-Val de Loire -6.475e+03 3.633e+03 -1.782 0.074845 .
## REGIONSGrand Est 4.265e+03 2.390e+03 1.785 0.074458 .
## REGIONSHauts-de-France 3.238e+03 2.880e+03 1.124 0.260976
## REGIONSIle-de-France 1.005e+03 1.120e+04 0.090 0.928539
## REGIONSNormandie -1.200e+03 2.674e+03 -0.449 0.653708
## REGIONSNouvelle-Aquitaine 4.244e+03 2.321e+03 1.828 0.067702 .
## REGIONSOccitanie 2.773e+03 2.776e+03 0.999 0.317873
## REGIONSPACA -9.341e+03 7.518e+03 -1.242 0.214231
## REGIONSPays de la Loire -5.147e+03 2.754e+03 -1.869 0.061737 .
## SFP -1.426e+02 1.967e+01 -7.246 6.22e-13 ***
## SUT1GCU -2.251e+02 4.570e+01 -4.927 9.09e-07 ***
## EBEXP 9.157e-02 1.053e-02 8.697 < 2e-16 ***
## CHRPV 2.572e+00 1.603e-01 16.044 < 2e-16 ***
## Total_charges_fossiles 8.968e-01 9.728e-02 9.219 < 2e-16 ***
## ELECT 3.341e+00 2.247e-01 14.865 < 2e-16 ***
## Intra_conso_animaux -5.087e-01 6.374e-02 -7.980 2.51e-15 ***
## UGBTO_classe50-99 3.073e+03 2.603e+03 1.181 0.237865
## UGBTO_classe100-149 1.037e+04 3.047e+03 3.402 0.000682 ***
## UGBTO_classe150-199 1.945e+04 3.694e+03 5.264 1.57e-07 ***
## UGBTO_classe200-249 3.104e+04 4.405e+03 7.046 2.57e-12 ***
## UGBTO_classe250-299 3.148e+04 5.117e+03 6.153 9.25e-10 ***
## UGBTO_classe300-349 4.068e+04 6.561e+03 6.200 6.92e-10 ***
## UGBTO_classe350-399 7.941e+04 7.974e+03 9.959 < 2e-16 ***
## SAUTI_classe50-99 7.492e+03 3.005e+03 2.493 0.012750 *
## SAUTI_classe100-149 1.390e+04 4.078e+03 3.410 0.000664 ***
## SAUTI_classe150-199 2.168e+04 5.544e+03 3.911 9.52e-05 ***
## SAUTI_classe200-249 2.289e+04 7.220e+03 3.171 0.001544 **
## SAUTI_classe250-299 3.985e+04 9.306e+03 4.282 1.94e-05 ***
## SAUTI_classe300-349 3.885e+04 1.093e+04 3.555 0.000387 ***
## SAUTI_classe350-399 4.214e+04 1.411e+04 2.987 0.002858 **
## SAUTI_classe400-449 3.323e+04 1.768e+04 1.880 0.060287 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 26880 on 1893 degrees of freedom
## Multiple R-squared: 0.7585, Adjusted R-squared: 0.7539
## F-statistic: 165.1 on 36 and 1893 DF, p-value: < 2.2e-16Pour expliquer l’achat d’aliment (grossier+concentré), on regarde les variables suivantes : OTEXE, REGIONS, SFP, SUT1GCU, EBEXP, CHRPV, Total_charges_fossiles, ELECT, Intra_conso_animaux, UGBTO_classe et SAUTI_classe.
A l’aide d’un test de nullité des coefficients on observe que :
- Il y a une différence significative entre les exploitations spécialisées bovin lait (OTEXE 45, pris en référence) et toutes les autres exploitations bovines (viande, mixte, polycuture-élevage). Les exploitations laitières sont celles achetant le plus d’aliment extérieur, ce qui est cohérent avec des productions plus souvent en bâtiment (donc pas d’apport d’herbe au pâturage) et importantes (production laitière demandantes en nutriment, besoin d’ajout de concentrés). La différence entre les exploitations spécialisées en lait et celle de polyculture est du à un apport moins important de fourrage et concentré en intraconommation, qui sont cependant fréquent en polyculture.
## $lsmeans
## OTEXE lsmean SE df lower.CL upper.CL
## 45 81480 4120 1893 73392 89568
## 46 65338 4100 1893 57298 73378
## 47 70763 4260 1893 62407 79119
## 83 61740 4600 1893 52722 70757
##
## Results are averaged over the levels of: REGIONS, UGBTO_classe, SAUTI_classe
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## OTEXE45 - OTEXE46 16142 1950 1893 8.274 <.0001
## OTEXE45 - OTEXE47 10717 2400 1893 4.461 0.0001
## OTEXE45 - OTEXE83 19741 2490 1893 7.932 <.0001
## OTEXE46 - OTEXE47 -5425 2470 1893 -2.194 0.1251
## OTEXE46 - OTEXE83 3598 2620 1893 1.376 0.5148
## OTEXE47 - OTEXE83 9024 3140 1893 2.875 0.0213
##
## Results are averaged over the levels of: REGIONS, UGBTO_classe, SAUTI_classe
## P value adjustment: tukey method for comparing a family of 4 estimatesOn voit ensuite à l’aide d’une comparaison deux à deux sur la variable OTEXE (ci-dessus), qu’on trouve également une différence significative entre les exploitations bovin mixtes (lait-viande) et les exploitations en polyculture, avec significativement moins d’aliment exterieurs consommés par les exploitations de polyculture. On voit donc globalement que c’est la production laitière qui entraine le plus fréquemment le besoin d’achat d’aliment à l’extérieur, les exploitations viandes n’étant pas significativement différentes des exploitations en polyculture en termes d’achat d’aliment extérieur.
-La SFP et la SUT1GCU ont un impact significativement négatif sur l’achat d’aliment ce qui est également cohérent car on estime que la présence de surfaces fourragères et de surfaces culturales représentent un apport d’aliment interne à l’exploitation, et donc un besoin d’achat en aliment exterieur moins important. De même on voit un impact fortement négatif de l’intraconsommation d’aliment pour animaux sur l’achat d’aliment ce qui est cohérent.
On observe globalement que les charges énérgétiques (fossiles et éléctriques) ont un impact positif sur les achats d’aliment, ceci peut s’expliquer simplement par la correlation entre achat d’aliment et taille d’exploitation : en effet, plus l’exploitation est grande plus les charges énergétiques le seront mais plus les troupeaux le seront aussi. Ceci explique également l’effet positif des charges vétérinaires qui sont aussi liées aux tailles de cheptel. On voit cet effet des tailles de cheptel avec les effets significativement positifs de toutes les classes UGB supérieures à celle de référence qui est la plus faible.
On voit un effet positif important de l’EBEXP sur l’achat d’aliment, il semblerait donc que les exploitations achetant de l’aliment sont les exploitations les plus productives.
-Finalement, on observe que la seule région ayant un impact sur l’achat d’aliment est la Bretagne, on voit que moins d’aliment sont achetés en Bretagne en comparaison avec la région Auvergne-Rhône-Alpes. Ceci peut s’expliquer par la présence de plus de polyculture élevage en Bretagne ce qui augmente la part d’intraconsommation dans cette région, tandis qu’en Auvergne-Rhône-Alpes les terres sont moins propices à la polyculture (terres plus pauvres).
On va réaliser les mêmes modèles, mais pour l’autonomie : nous prenons comme variable à expliquer le ratio “Achats d’aliments / (Produit total - Subventions + Intraconsommation totale)”, c’est à dire que l’on divise les achats d’aliments par ce qui a été vendu (moins les subventions) et ce qui a été intraconsommé.
##
## Call:
## lm(formula = ratio_Achat ~ OTEXE + REGIONS + SFP + SUT1GCU +
## VAPBR + EBEXP + CHRPV + Total_charges_fossiles + ELECT +
## PITOT + INTVP + Intra_conso_animaux + UGBTO_classe + SAUTI_classe,
## data = data_mod_Bovins)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.24424 -0.05295 -0.00941 0.04087 0.82326
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.145e-01 1.175e-02 18.248 < 2e-16 ***
## OTEXE46 -4.179e-02 6.326e-03 -6.606 5.12e-11 ***
## OTEXE47 -9.695e-03 7.753e-03 -1.250 0.21130
## OTEXE83 -7.656e-02 8.046e-03 -9.515 < 2e-16 ***
## REGIONSBourgogne-Franche-Comté -2.051e-02 7.557e-03 -2.714 0.00671 **
## REGIONSBretagne -7.239e-02 9.740e-03 -7.432 1.61e-13 ***
## REGIONSCentre-Val de Loire -2.507e-02 1.178e-02 -2.128 0.03347 *
## REGIONSGrand Est -1.938e-02 7.838e-03 -2.473 0.01348 *
## REGIONSHauts-de-France -1.723e-02 9.487e-03 -1.817 0.06944 .
## REGIONSIle-de-France -1.106e-02 3.637e-02 -0.304 0.76111
## REGIONSNormandie -4.156e-02 8.970e-03 -4.633 3.85e-06 ***
## REGIONSNouvelle-Aquitaine -1.032e-02 7.500e-03 -1.376 0.16910
## REGIONSOccitanie 1.150e-02 8.938e-03 1.287 0.19827
## REGIONSPACA -1.544e-02 2.423e-02 -0.637 0.52409
## REGIONSPays de la Loire -4.067e-02 9.040e-03 -4.499 7.23e-06 ***
## SFP -3.285e-04 6.444e-05 -5.098 3.77e-07 ***
## SUT1GCU -7.289e-04 1.472e-04 -4.952 8.00e-07 ***
## VAPBR -2.855e-07 9.266e-08 -3.081 0.00209 **
## EBEXP 1.046e-07 9.319e-08 1.123 0.26176
## CHRPV 4.092e-06 5.181e-07 7.898 4.76e-15 ***
## Total_charges_fossiles 3.161e-07 3.139e-07 1.007 0.31398
## ELECT 1.428e-06 7.376e-07 1.936 0.05302 .
## PITOT -8.533e-07 6.859e-07 -1.244 0.21366
## INTVP 2.141e-05 4.273e-05 0.501 0.61645
## Intra_conso_animaux -5.677e-07 7.183e-07 -0.790 0.42943
## UGBTO_classe50-99 2.162e-02 8.399e-03 2.575 0.01011 *
## UGBTO_classe100-149 3.966e-02 9.860e-03 4.023 5.98e-05 ***
## UGBTO_classe150-199 5.177e-02 1.194e-02 4.334 1.54e-05 ***
## UGBTO_classe200-249 6.410e-02 1.425e-02 4.500 7.21e-06 ***
## UGBTO_classe250-299 5.408e-02 1.656e-02 3.266 0.00111 **
## UGBTO_classe300-349 6.329e-02 2.120e-02 2.985 0.00287 **
## UGBTO_classe350-399 9.069e-02 2.585e-02 3.509 0.00046 ***
## SAUTI_classe50-99 3.607e-03 9.680e-03 0.373 0.70945
## SAUTI_classe100-149 9.783e-03 1.314e-02 0.744 0.45668
## SAUTI_classe150-199 3.618e-02 1.789e-02 2.022 0.04329 *
## SAUTI_classe200-249 4.949e-02 2.328e-02 2.126 0.03366 *
## SAUTI_classe250-299 9.172e-02 3.000e-02 3.057 0.00226 **
## SAUTI_classe300-349 1.125e-01 3.522e-02 3.195 0.00142 **
## SAUTI_classe350-399 1.424e-01 4.555e-02 3.126 0.00180 **
## SAUTI_classe400-449 2.024e-01 5.707e-02 3.547 0.00040 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.08656 on 1890 degrees of freedom
## Multiple R-squared: 0.2739, Adjusted R-squared: 0.2589
## F-statistic: 18.28 on 39 and 1890 DF, p-value: < 2.2e-16## Start: AIC=-9405.68
## ratio_Achat ~ OTEXE + REGIONS + SFP + SUT1GCU + VAPBR + EBEXP +
## CHRPV + Total_charges_fossiles + ELECT + PITOT + INTVP +
## Intra_conso_animaux + UGBTO_classe + SAUTI_classe## 2 3 5 22 28 174 183 270 271 275 300 315 403 431 484 605
## 2 3 5 22 28 174 183 270 271 275 300 315 403 431 484 605
## 618 648 678 720 804 881 897 899 920 929 965 979 1029 1034 1078 1103
## 618 648 678 720 804 881 897 899 920 929 965 979 1029 1034 1078 1103
## 1116 1121 1196 1201 1304 1311 1317 1334 1438 1445 1487 1490 1495 1501 1527 1531
## 1116 1121 1196 1201 1304 1311 1317 1334 1438 1445 1487 1490 1495 1501 1527 1531
## 1537 1539 1561 1565 1582 1589 1628 1644 1673 1681 1719 1722 1789 1807 1813 1824
## 1537 1539 1561 1565 1582 1589 1628 1644 1673 1681 1719 1722 1789 1807 1813 1824
## 1845 1851 1868 1897 1907 1917 1918 1920 1921 1924 1929
## 1845 1851 1868 1897 1907 1917 1918 1920 1921 1924 1929##
## Call:
## lm(formula = ratio_Achat ~ OTEXE + REGIONS + SFP + SUT1GCU +
## VAPBR + EBEXP + CHRPV + Total_charges_fossiles + ELECT +
## PITOT + INTVP + Intra_conso_animaux + UGBTO_classe + SAUTI_classe,
## data = data_mod_Bovins)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.24424 -0.05295 -0.00941 0.04087 0.82326
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.145e-01 1.175e-02 18.248 < 2e-16 ***
## OTEXE46 -4.179e-02 6.326e-03 -6.606 5.12e-11 ***
## OTEXE47 -9.695e-03 7.753e-03 -1.250 0.21130
## OTEXE83 -7.656e-02 8.046e-03 -9.515 < 2e-16 ***
## REGIONSBourgogne-Franche-Comté -2.051e-02 7.557e-03 -2.714 0.00671 **
## REGIONSBretagne -7.239e-02 9.740e-03 -7.432 1.61e-13 ***
## REGIONSCentre-Val de Loire -2.507e-02 1.178e-02 -2.128 0.03347 *
## REGIONSGrand Est -1.938e-02 7.838e-03 -2.473 0.01348 *
## REGIONSHauts-de-France -1.723e-02 9.487e-03 -1.817 0.06944 .
## REGIONSIle-de-France -1.106e-02 3.637e-02 -0.304 0.76111
## REGIONSNormandie -4.156e-02 8.970e-03 -4.633 3.85e-06 ***
## REGIONSNouvelle-Aquitaine -1.032e-02 7.500e-03 -1.376 0.16910
## REGIONSOccitanie 1.150e-02 8.938e-03 1.287 0.19827
## REGIONSPACA -1.544e-02 2.423e-02 -0.637 0.52409
## REGIONSPays de la Loire -4.067e-02 9.040e-03 -4.499 7.23e-06 ***
## SFP -3.285e-04 6.444e-05 -5.098 3.77e-07 ***
## SUT1GCU -7.289e-04 1.472e-04 -4.952 8.00e-07 ***
## VAPBR -2.855e-07 9.266e-08 -3.081 0.00209 **
## EBEXP 1.046e-07 9.319e-08 1.123 0.26176
## CHRPV 4.092e-06 5.181e-07 7.898 4.76e-15 ***
## Total_charges_fossiles 3.161e-07 3.139e-07 1.007 0.31398
## ELECT 1.428e-06 7.376e-07 1.936 0.05302 .
## PITOT -8.533e-07 6.859e-07 -1.244 0.21366
## INTVP 2.141e-05 4.273e-05 0.501 0.61645
## Intra_conso_animaux -5.677e-07 7.183e-07 -0.790 0.42943
## UGBTO_classe50-99 2.162e-02 8.399e-03 2.575 0.01011 *
## UGBTO_classe100-149 3.966e-02 9.860e-03 4.023 5.98e-05 ***
## UGBTO_classe150-199 5.177e-02 1.194e-02 4.334 1.54e-05 ***
## UGBTO_classe200-249 6.410e-02 1.425e-02 4.500 7.21e-06 ***
## UGBTO_classe250-299 5.408e-02 1.656e-02 3.266 0.00111 **
## UGBTO_classe300-349 6.329e-02 2.120e-02 2.985 0.00287 **
## UGBTO_classe350-399 9.069e-02 2.585e-02 3.509 0.00046 ***
## SAUTI_classe50-99 3.607e-03 9.680e-03 0.373 0.70945
## SAUTI_classe100-149 9.783e-03 1.314e-02 0.744 0.45668
## SAUTI_classe150-199 3.618e-02 1.789e-02 2.022 0.04329 *
## SAUTI_classe200-249 4.949e-02 2.328e-02 2.126 0.03366 *
## SAUTI_classe250-299 9.172e-02 3.000e-02 3.057 0.00226 **
## SAUTI_classe300-349 1.125e-01 3.522e-02 3.195 0.00142 **
## SAUTI_classe350-399 1.424e-01 4.555e-02 3.126 0.00180 **
## SAUTI_classe400-449 2.024e-01 5.707e-02 3.547 0.00040 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.08656 on 1890 degrees of freedom
## Multiple R-squared: 0.2739, Adjusted R-squared: 0.2589
## F-statistic: 18.28 on 39 and 1890 DF, p-value: < 2.2e-16## Anova Table (Type II tests)
##
## Response: ratio_Achat
## Sum Sq Df F value Pr(>F)
## OTEXE 0.8044 3 35.7903 < 2.2e-16 ***
## REGIONS 0.6147 11 7.4594 1.174e-12 ***
## SFP 0.1947 1 25.9922 3.771e-07 ***
## SUT1GCU 0.1837 1 24.5232 7.995e-07 ***
## VAPBR 0.0711 1 9.4910 0.0020946 **
## EBEXP 0.0094 1 1.2602 0.2617634
## CHRPV 0.4674 1 62.3801 4.762e-15 ***
## Total_charges_fossiles 0.0076 1 1.0144 0.3139757
## ELECT 0.0281 1 3.7481 0.0530179 .
## PITOT 0.0116 1 1.5475 0.2136603
## INTVP 0.0019 1 0.2510 0.6164503
## Intra_conso_animaux 0.0047 1 0.6246 0.4294311
## UGBTO_classe 0.1918 7 3.6581 0.0006265 ***
## SAUTI_classe 0.1932 8 3.2232 0.0012043 **
## Residuals 14.1599 1890
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Call:
## lm(formula = ratio_Achat ~ OTEXE + REGIONS + SFP + SUT1GCU +
## EBEXP + CHRPV + Intra_conso_animaux + UGBTO_classe + SAUTI_classe,
## data = data_mod_Bovins)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.25701 -0.05408 -0.00865 0.04046 0.82224
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.248e-01 1.129e-02 19.910 < 2e-16 ***
## OTEXE46 -4.500e-02 5.700e-03 -7.895 4.85e-15 ***
## OTEXE47 -1.080e-02 7.629e-03 -1.416 0.156955
## OTEXE83 -7.827e-02 7.721e-03 -10.137 < 2e-16 ***
## REGIONSBourgogne-Franche-Comté -2.455e-02 7.490e-03 -3.277 0.001067 **
## REGIONSBretagne -8.040e-02 9.524e-03 -8.442 < 2e-16 ***
## REGIONSCentre-Val de Loire -3.061e-02 1.171e-02 -2.615 0.008984 **
## REGIONSGrand Est -2.521e-02 7.677e-03 -3.284 0.001044 **
## REGIONSHauts-de-France -2.436e-02 9.254e-03 -2.632 0.008562 **
## REGIONSIle-de-France -2.633e-02 3.615e-02 -0.728 0.466468
## REGIONSNormandie -5.045e-02 8.622e-03 -5.852 5.72e-09 ***
## REGIONSNouvelle-Aquitaine -1.214e-02 7.485e-03 -1.622 0.105054
## REGIONSOccitanie 1.210e-02 8.954e-03 1.351 0.176851
## REGIONSPACA -1.304e-02 2.426e-02 -0.537 0.591106
## REGIONSPays de la Loire -4.732e-02 8.840e-03 -5.353 9.72e-08 ***
## SFP -3.025e-04 6.315e-05 -4.791 1.79e-06 ***
## SUT1GCU -6.943e-04 1.433e-04 -4.844 1.37e-06 ***
## EBEXP -1.487e-07 3.268e-08 -4.552 5.66e-06 ***
## CHRPV 4.257e-06 5.054e-07 8.423 < 2e-16 ***
## Intra_conso_animaux -1.413e-06 2.050e-07 -6.889 7.62e-12 ***
## UGBTO_classe50-99 2.382e-02 8.391e-03 2.839 0.004572 **
## UGBTO_classe100-149 4.385e-02 9.760e-03 4.492 7.47e-06 ***
## UGBTO_classe150-199 5.751e-02 1.168e-02 4.923 9.26e-07 ***
## UGBTO_classe200-249 6.994e-02 1.384e-02 5.055 4.73e-07 ***
## UGBTO_classe250-299 6.014e-02 1.593e-02 3.774 0.000166 ***
## UGBTO_classe300-349 7.252e-02 2.062e-02 3.517 0.000447 ***
## UGBTO_classe350-399 9.971e-02 2.490e-02 4.003 6.48e-05 ***
## SAUTI_classe50-99 4.581e-03 9.680e-03 0.473 0.636072
## SAUTI_classe100-149 1.256e-02 1.310e-02 0.959 0.337554
## SAUTI_classe150-199 4.270e-02 1.780e-02 2.400 0.016502 *
## SAUTI_classe200-249 5.739e-02 2.313e-02 2.482 0.013167 *
## SAUTI_classe250-299 1.021e-01 2.983e-02 3.423 0.000632 ***
## SAUTI_classe300-349 1.214e-01 3.514e-02 3.456 0.000561 ***
## SAUTI_classe350-399 1.602e-01 4.534e-02 3.533 0.000421 ***
## SAUTI_classe400-449 2.136e-01 5.685e-02 3.757 0.000177 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.08677 on 1895 degrees of freedom
## Multiple R-squared: 0.2684, Adjusted R-squared: 0.2553
## F-statistic: 20.45 on 34 and 1895 DF, p-value: < 2.2e-16Pour expliquer le ratio d’achat (décrite par la relation : Achats d’aliments / (Produit total - Subventions + Intraconsommation totale)), on regarde les variables suivantes : OTEXE,REGIONS,SFP,SUT1GCU,EBEXP,CHRPV,Intra_conso_animaux,UGBTO_classe et SAUTI_classee.
A l’aide d’un test de nullité des coefficients on observe que :
- Il y a une différence significative entre les exploitations spécialisées bovin lait (OTEXE 45, pris en référence) et les exploitations ne produisant pas majoritairement de lait, soit viande et polycuture-élevage. Les exploitations laitières sont celles ayant la plus grande part d’aliment extérieur consommé par rapport à l’apport total d’aliment utilisé sur l’exploitation. Ceci est cohérent avec une forte utilisation de concentrés dans ce type d’exploitation. Regardons la comparaison deux à deux des OTEXE pour comparer les autres types d’exploitation entre elles :
## $lsmeans
## OTEXE lsmean SE df lower.CL upper.CL
## 45 0.247 0.0129 1895 0.221 0.272
## 46 0.202 0.0131 1895 0.176 0.227
## 47 0.236 0.0136 1895 0.209 0.262
## 83 0.168 0.0148 1895 0.139 0.197
##
## Results are averaged over the levels of: REGIONS, UGBTO_classe, SAUTI_classe
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## OTEXE45 - OTEXE46 0.0450 0.00570 1895 7.895 <.0001
## OTEXE45 - OTEXE47 0.0108 0.00763 1895 1.416 0.4895
## OTEXE45 - OTEXE83 0.0783 0.00772 1895 10.137 <.0001
## OTEXE46 - OTEXE47 -0.0342 0.00784 1895 -4.360 0.0001
## OTEXE46 - OTEXE83 0.0333 0.00842 1895 3.950 0.0005
## OTEXE47 - OTEXE83 0.0675 0.01010 1895 6.686 <.0001
##
## Results are averaged over the levels of: REGIONS, UGBTO_classe, SAUTI_classe
## P value adjustment: tukey method for comparing a family of 4 estimatesOn trouve également une différence significative entre les exploitations bovin mixtes (lait-viande) et les exploitations en polyculture, avec significativement moins d’aliment exterieurs consommés par les exploitations de polyculture, ce qui est cohérent avec la production d’aliment directement sur l’exploitation. On note également une plus grande part d’aliment achetée à l’exterieur en mixte qu’en viande ce qui tend à faire penser que c’est la production laitière qui a tendance à tirer la ratio vers des achats exterieurs importants (plus demandante en concentrés, qui sont plus rarement produits sur l’exploitation en particulier dans les régions tournées vers les systèmes herbagers comme l’Auvergne-Rhône-Alpes).
-La SFP et la SUT1GCU ont un impact significativement négatif sur le ratio d’achat, soit moins d’achat exterieur lorsque l’on a une production fourragère ou de grandes cultures sur l’exploitation. Ceci est également cohérent car on estime que la présence de surfaces fourragères et de surfaces culturales représentent un apport d’aliment interne à l’exploitation, et donc un besoin d’achat en aliment exterieur moins important. De même on voit un impact fortement négatif de l’intraconsommation d’aliment pour animaux sur l’achat d’aliment ce qui est cohérent.
- On voit un effet de la taille des cheptel avec des effets significativement positifs de toutes les classes UGB supérieures à celle de référence qui est la plus faible.
## $lsmeans
## UGBTO_classe lsmean SE df lower.CL upper.CL
## 0-49 0.160 0.0145 1895 0.131 0.188
## 50-99 0.183 0.0131 1895 0.158 0.209
## 100-149 0.204 0.0131 1895 0.178 0.229
## 150-199 0.217 0.0137 1895 0.190 0.244
## 200-249 0.230 0.0147 1895 0.201 0.258
## 250-299 0.220 0.0160 1895 0.188 0.251
## 300-349 0.232 0.0198 1895 0.193 0.271
## 350-399 0.259 0.0243 1895 0.212 0.307
##
## Results are averaged over the levels of: OTEXE, REGIONS, SAUTI_classe
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## (0-49) - (50-99) -0.02382 0.00839 1895 -2.839 0.0862
## (0-49) - (100-149) -0.04385 0.00976 1895 -4.492 0.0002
## (0-49) - (150-199) -0.05751 0.01170 1895 -4.923 <.0001
## (0-49) - (200-249) -0.06994 0.01380 1895 -5.055 <.0001
## (0-49) - (250-299) -0.06014 0.01590 1895 -3.774 0.0041
## (0-49) - (300-349) -0.07252 0.02060 1895 -3.517 0.0106
## (0-49) - (350-399) -0.09971 0.02490 1895 -4.003 0.0017
## (50-99) - (100-149) -0.02002 0.00586 1895 -3.420 0.0148
## (50-99) - (150-199) -0.03369 0.00817 1895 -4.126 0.0010
## (50-99) - (200-249) -0.04611 0.01070 1895 -4.293 0.0005
## (50-99) - (250-299) -0.03631 0.01320 1895 -2.744 0.1102
## (50-99) - (300-349) -0.04870 0.01840 1895 -2.640 0.1424
## (50-99) - (350-399) -0.07588 0.02310 1895 -3.279 0.0236
## (100-149) - (150-199) -0.01367 0.00686 1895 -1.992 0.4877
## (100-149) - (200-249) -0.02609 0.00936 1895 -2.789 0.0983
## (100-149) - (250-299) -0.01629 0.01190 1895 -1.368 0.8716
## (100-149) - (300-349) -0.02868 0.01730 1895 -1.659 0.7139
## (100-149) - (350-399) -0.05586 0.02210 1895 -2.523 0.1864
## (150-199) - (200-249) -0.01242 0.00900 1895 -1.380 0.8664
## (150-199) - (250-299) -0.00262 0.01130 1895 -0.232 1.0000
## (150-199) - (300-349) -0.01501 0.01660 1895 -0.906 0.9856
## (150-199) - (350-399) -0.04219 0.02150 1895 -1.966 0.5050
## (200-249) - (250-299) 0.00980 0.01180 1895 0.830 0.9914
## (200-249) - (300-349) -0.00259 0.01660 1895 -0.156 1.0000
## (200-249) - (350-399) -0.02977 0.02140 1895 -1.393 0.8608
## (250-299) - (300-349) -0.01239 0.01710 1895 -0.723 0.9963
## (250-299) - (350-399) -0.03957 0.02180 1895 -1.812 0.6118
## (300-349) - (350-399) -0.02718 0.02400 1895 -1.133 0.9497
##
## Results are averaged over the levels of: OTEXE, REGIONS, SAUTI_classe
## P value adjustment: tukey method for comparing a family of 8 estimatesIl semblerait que pour les classes faibles et intermédiairesd’UGB (jusqu’à 400 UGB) on ait un effet significatif du nombre d’UGB avec les classes les plus élevées qui utilisent une plus grande part d’aliments achetés à l’extérieur en comparaison avec les classes plus faibles. On perdrait donc en autonomie alimentaire en augmentant la taille des cheptels.
- On voit un effet très légérement négatif de l’EBEXP sur le ratio d’aliment alors que celui-ci avait un effet positif sur la quantité d’aliment acheté. Il semblerait donc que l’autonomie d’une exploitation est à mettre en lien avec sa rentabilité. Etre plus autonome alimentairement parlant favorise un bon EBE, même avec des achats importants dans l’absolu (par exemple pour les grandes exploitations).
Modèle linéaire sur la grande classe Petits ruminants :
On a proposé un modèle linéaire pour expliquer les achats alimentaires des élevages de petits ruminants. Il s’agit des exploitations mixtes et spécialisées.
##
## Call:
## lm(formula = ACHA_SOMME ~ PBUCE + OTE64F + REGIO + SFP + SUT1GCU +
## CHRGTO + PBRTO + VAPBR + EBEXP + UGBBO_classe + CHREN + CHRPH +
## CHRAC + CHRPV + Total_charges_fossiles + ELECT + LANMX +
## PITOT + ACHAC + ACHAG + Intra_conso_animaux + SAUTI_classe +
## ratio_Achat, data = data_mod_petits_ruminants)
##
## Residuals:
## Min 1Q Median 3Q Max
## -3.912e-10 -5.870e-12 7.700e-13 7.260e-12 1.582e-10
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -3.705e-11 2.259e-11 -1.640e+00 0.1019
## PBUCE -3.986e-17 4.619e-17 -8.630e-01 0.3888
## OTE64F4820 2.001e-11 8.723e-12 2.294e+00 0.0224 *
## OTE64F4830 -7.024e-12 7.271e-12 -9.660e-01 0.3347
## OTE64F4840 1.757e-11 1.319e-11 1.332e+00 0.1839
## REGIO23 -4.974e-11 3.734e-11 -1.332e+00 0.1838
## REGIO24 -2.345e-11 2.364e-11 -9.920e-01 0.3219
## REGIO25 1.878e-11 2.733e-11 6.870e-01 0.4923
## REGIO26 2.922e-11 2.407e-11 1.214e+00 0.2257
## REGIO41 1.117e-11 2.349e-11 4.750e-01 0.6349
## REGIO42 2.988e-11 3.794e-11 7.880e-01 0.4315
## REGIO43 6.860e-12 3.059e-11 2.240e-01 0.8227
## REGIO52 2.841e-11 2.579e-11 1.102e+00 0.2714
## REGIO53 2.424e-11 2.903e-11 8.350e-01 0.4043
## REGIO54 3.036e-11 2.275e-11 1.335e+00 0.1829
## REGIO72 2.008e-11 2.274e-11 8.830e-01 0.3780
## REGIO73 1.184e-11 2.280e-11 5.190e-01 0.6040
## REGIO74 1.741e-11 2.306e-11 7.550e-01 0.4508
## REGIO82 2.537e-11 2.223e-11 1.142e+00 0.2545
## REGIO83 1.136e-11 2.322e-11 4.890e-01 0.6249
## REGIO91 2.131e-11 2.302e-11 9.260e-01 0.3552
## REGIO93 1.184e-11 2.314e-11 5.120e-01 0.6092
## SFP -1.575e-14 3.466e-14 -4.550e-01 0.6497
## SUT1GCU 6.986e-14 8.050e-14 8.680e-01 0.3862
## CHRGTO 1.422e-16 9.385e-17 1.515e+00 0.1306
## PBRTO -5.413e-17 1.068e-16 -5.070e-01 0.6126
## VAPBR -1.328e-17 7.668e-17 -1.730e-01 0.8626
## EBEXP -3.064e-17 1.073e-16 -2.860e-01 0.7754
## UGBBO_classe50-99 -1.039e-11 1.068e-11 -9.740e-01 0.3309
## UGBBO_classe100-149 -9.718e-11 1.486e-11 -6.541e+00 2.34e-10 ***
## UGBBO_classe150-199 -3.681e-11 2.547e-11 -1.445e+00 0.1493
## UGBBO_classe200-249 -4.587e-11 2.520e-11 -1.820e+00 0.0697 .
## CHREN 4.960e-16 4.325e-16 1.147e+00 0.2523
## CHRPH 1.538e-15 7.241e-16 2.124e+00 0.0344 *
## CHRAC -1.030e-15 5.354e-16 -1.924e+00 0.0553 .
## CHRPV 6.758e-18 6.640e-16 1.000e-02 0.9919
## Total_charges_fossiles 6.078e-16 3.841e-16 1.582e+00 0.1145
## ELECT -1.569e-15 1.020e-15 -1.538e+00 0.1250
## LANMX 8.392e-16 1.065e-14 7.900e-02 0.9372
## PITOT -7.914e-16 1.260e-15 -6.280e-01 0.5303
## ACHAC 1.000e+00 5.212e-16 1.919e+15 < 2e-16 ***
## ACHAG 1.000e+00 2.080e-16 4.809e+15 < 2e-16 ***
## Intra_conso_animaux 8.104e-16 1.266e-15 6.400e-01 0.5227
## SAUTI_classe50-99 7.549e-12 5.169e-12 1.461e+00 0.1451
## SAUTI_classe100-149 8.235e-12 7.572e-12 1.088e+00 0.2776
## SAUTI_classe150-199 1.311e-11 1.004e-11 1.306e+00 0.1924
## SAUTI_classe200-249 1.300e-11 1.509e-11 8.610e-01 0.3897
## SAUTI_classe250-299 -7.197e-13 1.656e-11 -4.300e-02 0.9654
## SAUTI_classe300-349 1.145e-11 1.904e-11 6.010e-01 0.5480
## SAUTI_classe350-399 1.181e-11 2.349e-11 5.030e-01 0.6155
## SAUTI_classe400-449 1.970e-11 2.617e-11 7.530e-01 0.4522
## ratio_Achat 4.789e-13 5.091e-12 9.400e-02 0.9251
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 3.007e-11 on 329 degrees of freedom
## Multiple R-squared: 1, Adjusted R-squared: 1
## F-statistic: 1.34e+31 on 51 and 329 DF, p-value: < 2.2e-16On a pris un modèle complet puis nous avons opéré à une sélection de variable pour abourir à un modèle sans interaction avec les variables : Production brute standards en euros (PBUCE), Surface fourragère principale (SFP en ha), Charges totales (CHRTO) en euros, Produit brut total (PBRTO) en euros, Valeurs ajoutées en euros (VAPBR), EBE d’exploitation (EBEXP) en euros, Charges d’engrais en euros (CHREN), les charges en produits de défenses végétales en euros (CHRPH), les charges vétérinaires en euros (CHRPV), les charges en énergies fossiles (Total_charges_fossiles), l’intraconsommation en euros (PITOT), l’orientation technico économique (OTE64F), UGB bovin (UGBBO_classe) et les regions.
## Start: AIC=-18413.27
## ACHA_SOMME ~ PBUCE + OTE64F + REGIO + SFP + SUT1GCU + CHRGTO +
## PBRTO + VAPBR + EBEXP + UGBBO_classe + CHREN + CHRPH + CHRAC +
## CHRPV + Total_charges_fossiles + ELECT + LANMX + PITOT +
## ACHAC + ACHAG + Intra_conso_animaux + SAUTI_classe + ratio_Achat## 1 2 3 6 11 12 14 15 21 25 30 31 35 55 59 61 81 116 121 124
## 1 2 3 6 11 12 14 15 21 25 30 31 35 55 59 61 81 116 121 124
## 163 360 373
## 163 360 373##
## Call:
## lm(formula = ACHA_SOMME ~ PBUCE + OTE64F + SFP + CHRGTO + PBRTO +
## VAPBR + EBEXP + UGBBO_classe + CHREN + CHRPH + CHRPV + Total_charges_fossiles +
## PITOT + REGIONS, data = data_mod_petits_ruminants)
##
## Residuals:
## Min 1Q Median 3Q Max
## -64696 -7316 -421 6307 82073
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -7.264e+03 2.884e+03 -2.519 0.012224 *
## PBUCE 1.399e-01 2.210e-02 6.330 7.74e-10 ***
## OTE64F4820 -7.475e+03 4.657e+03 -1.605 0.109403
## OTE64F4830 7.647e+03 3.087e+03 2.478 0.013713 *
## OTE64F4840 -9.912e+03 6.711e+03 -1.477 0.140578
## SFP -1.240e+01 9.175e+00 -1.351 0.177554
## CHRGTO 1.382e-01 4.782e-02 2.891 0.004092 **
## PBRTO 1.221e-01 5.651e-02 2.161 0.031413 *
## VAPBR -9.934e-02 3.647e-02 -2.724 0.006792 **
## EBEXP -1.186e-01 5.548e-02 -2.138 0.033211 *
## UGBBO_classe50-99 5.910e+03 5.714e+03 1.034 0.301728
## UGBBO_classe100-149 2.773e+04 7.508e+03 3.694 0.000257 ***
## UGBBO_classe150-199 -1.405e+04 1.352e+04 -1.039 0.299661
## UGBBO_classe200-249 2.132e+04 1.306e+04 1.632 0.103640
## CHREN -5.946e-01 2.281e-01 -2.607 0.009537 **
## CHRPH -6.075e-01 3.621e-01 -1.678 0.094337 .
## CHRPV 1.795e+00 3.359e-01 5.343 1.68e-07 ***
## Total_charges_fossiles -1.309e+00 1.849e-01 -7.079 8.39e-12 ***
## PITOT -7.727e-01 1.109e-01 -6.970 1.65e-11 ***
## REGIONSBourgogne-Franche-Comté 5.472e+03 5.484e+03 0.998 0.319085
## REGIONSBretagne -1.091e+03 1.034e+04 -0.106 0.915984
## REGIONSCentre-Val de Loire 8.811e+03 6.068e+03 1.452 0.147423
## REGIONSGrand Est 1.334e+03 4.750e+03 0.281 0.778949
## REGIONSHauts-de-France -1.285e+03 1.190e+04 -0.108 0.914067
## REGIONSNormandie -3.823e+03 8.613e+03 -0.444 0.657452
## REGIONSNouvelle-Aquitaine 6.468e+03 2.881e+03 2.245 0.025412 *
## REGIONSOccitanie 8.708e+03 2.958e+03 2.944 0.003463 **
## REGIONSPACA -1.928e+03 4.192e+03 -0.460 0.645926
## REGIONSPays de la Loire 1.637e+03 7.439e+03 0.220 0.825933
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 16310 on 340 degrees of freedom
## Multiple R-squared: 0.852, Adjusted R-squared: 0.8398
## F-statistic: 69.91 on 28 and 340 DF, p-value: < 2.2e-16## Anova Table (Type II tests)
##
## Response: ACHA_SOMME
## Sum Sq Df F value Pr(>F)
## PBUCE 1.0662e+10 1 40.0694 7.744e-10 ***
## OTE64F 3.3709e+09 3 4.2227 0.005974 **
## SFP 4.8577e+08 1 1.8255 0.177554
## CHRGTO 2.2234e+09 1 8.3557 0.004092 **
## PBRTO 1.2423e+09 1 4.6688 0.031413 *
## VAPBR 1.9738e+09 1 7.4176 0.006792 **
## EBEXP 1.2166e+09 1 4.5720 0.033211 *
## UGBBO_classe 4.8821e+09 4 4.5868 0.001278 **
## CHREN 1.8085e+09 1 6.7963 0.009537 **
## CHRPH 7.4891e+08 1 2.8145 0.094337 .
## CHRPV 7.5974e+09 1 28.5513 1.675e-07 ***
## Total_charges_fossiles 1.3333e+10 1 50.1052 8.391e-12 ***
## PITOT 1.2927e+10 1 48.5795 1.655e-11 ***
## REGIONS 4.4457e+09 10 1.6707 0.086182 .
## Residuals 9.0472e+10 340
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## (Intercept) PBUCE
## -7.263997e+03 1.398978e-01
## OTE64F4820 OTE64F4830
## -7.474719e+03 7.647183e+03
## OTE64F4840 SFP
## -9.912151e+03 -1.239685e+01
## CHRGTO PBRTO
## 1.382374e-01 1.221041e-01
## VAPBR EBEXP
## -9.933590e-02 -1.186222e-01
## UGBBO_classe50-99 UGBBO_classe100-149
## 5.910464e+03 2.773260e+04
## UGBBO_classe150-199 UGBBO_classe200-249
## -1.404780e+04 2.131801e+04
## CHREN CHRPH
## -5.946129e-01 -6.075295e-01
## CHRPV Total_charges_fossiles
## 1.794622e+00 -1.308792e+00
## PITOT REGIONSBourgogne-Franche-Comté
## -7.726587e-01 5.472257e+03
## REGIONSBretagne REGIONSCentre-Val de Loire
## -1.091265e+03 8.810991e+03
## REGIONSGrand Est REGIONSHauts-de-France
## 1.334377e+03 -1.285003e+03
## REGIONSNormandie REGIONSNouvelle-Aquitaine
## -3.822513e+03 6.468095e+03
## REGIONSOccitanie REGIONSPACA
## 8.707719e+03 -1.927821e+03
## REGIONSPays de la Loire
## 1.637269e+03Toutes nos variables explicatives finales ont un effet significatif sur la variable d’intérêt achat-somme.
Le modèle retenu : ACHA_SOMME ~ PBUCE + OTE64F + SFP + CHRGTO + PBRTO + VAPBR + EBEXP + UGBBO_classe + CHREN + CHRPH + CHRPV + Total_charges_fossiles + PITOT+ REGIONS
En valeur absolue, pour les variables économiques en euros, les variables des charges vétérinaires (coefficient positifs) et de produits de défense végétaux (coefficient négatif) sont ceux qui ont le plus d’effet sur la variable d’intérêt. En effet, plus l’agriculteur possède de culture, plus il va utiliser de produits phytosaniatires et donc moins acheter à l’extérieur. Et plus il y a de frais véto plus il y a d’animaux et donc plus il y a des achats alimentaires.
## $lsmeans
## UGBBO_classe lsmean SE df lower.CL upper.CL
## 0-49 32653 2870 340 27013 38293
## 50-99 38563 5330 340 28071 49056
## 100-149 60385 7320 340 45978 74793
## 150-199 18605 13000 340 -7033 44243
## 200-249 53971 13000 340 28481 79460
##
## Results are averaged over the levels of: OTE64F, REGIONS
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## (0-49) - (50-99) -5910 5710 340 -1.034 0.8393
## (0-49) - (100-149) -27733 7510 340 -3.694 0.0024
## (0-49) - (150-199) 14048 13500 340 1.039 0.8372
## (0-49) - (200-249) -21318 13100 340 -1.632 0.4781
## (50-99) - (100-149) -21822 7940 340 -2.747 0.0494
## (50-99) - (150-199) 19958 13600 340 1.466 0.5854
## (50-99) - (200-249) -15408 13500 340 -1.138 0.7863
## (100-149) - (150-199) 41780 14100 340 2.956 0.0274
## (100-149) - (200-249) 6415 14100 340 0.456 0.9910
## (150-199) - (200-249) -35366 17900 340 -1.979 0.2783
##
## Results are averaged over the levels of: OTE64F, REGIONS
## P value adjustment: tukey method for comparing a family of 5 estimates## $lsmeans
## OTE64F lsmean SE df lower.CL upper.CL
## 4810 43270 5270 340 32900 53641
## 4820 35796 5540 340 24890 46702
## 4830 50918 5680 340 39751 62084
## 4840 33358 6940 340 19701 47015
##
## Results are averaged over the levels of: UGBBO_classe, REGIONS
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## OTE64F4810 - OTE64F4820 7475 4660 340 1.605 0.3771
## OTE64F4810 - OTE64F4830 -7647 3090 340 -2.478 0.0653
## OTE64F4810 - OTE64F4840 9912 6710 340 1.477 0.4525
## OTE64F4820 - OTE64F4830 -15122 5470 340 -2.762 0.0307
## OTE64F4820 - OTE64F4840 2437 7360 340 0.331 0.9875
## OTE64F4830 - OTE64F4840 17559 6620 340 2.653 0.0414
##
## Results are averaged over the levels of: UGBBO_classe, REGIONS
## P value adjustment: tukey method for comparing a family of 4 estimatesPour les variables qualitatives
UGBBO_classe :
Les classes avec des différences significatives sont : 1) 0-49 et 100-149 UGB bovin 2)50-99 et 200-249 UGB bovin 3) 100-149 et 150-199 UGB bovin
Et pour les OTE64F :
On a des différences significatives entre les OTE : 1) 4830 et 4840 Exploitations caprines spécialisées et Exploitations d’herbivores 2) 4830 et 4820 Exploitations caprines spécialisée et Exploitations avec ovins et bovins combinés
Les hypothèses sont respectées une fois que les points aberrants ont été
retirés.
On s’intéresse ensuite à l’autonomie alimentaire à l’aide de ratio_achat
##
## Call:
## lm(formula = ratio_Achat ~ PBUCE + OTE64F + REGIO + SFP + SUT1GCU +
## CHRGTO + VAPBR + EBEXP + UGBBO_classe + CHREN + CHRPH + CHRAC +
## CHRPV + Total_charges_fossiles + ELECT + LANMX + Intra_conso_animaux +
## SAUTI_classe + CHRAC, data = data_mod_petits_ruminants)
##
## Residuals:
## Min 1Q Median 3Q Max
## -5.1590 -0.0779 -0.0083 0.0750 1.4108
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.371e-01 2.487e-01 0.953 0.3413
## PBUCE 9.323e-07 4.969e-07 1.876 0.0615 .
## OTE64F4820 -4.778e-03 9.838e-02 -0.049 0.9613
## OTE64F4830 -3.617e-02 7.873e-02 -0.459 0.6463
## OTE64F4840 -2.712e-02 1.434e-01 -0.189 0.8500
## REGIO24 6.294e-02 2.699e-01 0.233 0.8158
## REGIO25 6.532e-03 3.017e-01 0.022 0.9827
## REGIO26 1.279e-01 2.672e-01 0.479 0.6323
## REGIO41 8.369e-02 2.591e-01 0.323 0.7469
## REGIO42 8.122e-02 4.166e-01 0.195 0.8455
## REGIO43 3.648e-01 3.369e-01 1.083 0.2797
## REGIO52 -9.827e-02 2.835e-01 -0.347 0.7291
## REGIO53 -9.453e-02 3.189e-01 -0.296 0.7671
## REGIO54 -1.298e-01 2.510e-01 -0.517 0.6056
## REGIO72 6.616e-02 2.510e-01 0.264 0.7923
## REGIO73 4.494e-02 2.520e-01 0.178 0.8585
## REGIO74 1.989e-01 2.545e-01 0.781 0.4352
## REGIO82 8.178e-02 2.451e-01 0.334 0.7389
## REGIO83 9.549e-02 2.568e-01 0.372 0.7102
## REGIO91 8.977e-02 2.542e-01 0.353 0.7242
## REGIO93 7.541e-02 2.552e-01 0.295 0.7678
## SFP -9.766e-05 3.802e-04 -0.257 0.7975
## SUT1GCU -7.796e-04 8.793e-04 -0.887 0.3759
## CHRGTO -1.269e-07 4.749e-07 -0.267 0.7896
## VAPBR -5.926e-07 8.134e-07 -0.729 0.4668
## EBEXP -2.760e-07 7.892e-07 -0.350 0.7268
## UGBBO_classe50-99 -1.310e-01 1.184e-01 -1.106 0.2694
## UGBBO_classe100-149 -8.401e-02 1.613e-01 -0.521 0.6028
## UGBBO_classe150-199 -1.844e-01 2.856e-01 -0.646 0.5190
## UGBBO_classe200-249 1.290e-02 2.773e-01 0.047 0.9629
## CHREN -2.599e-06 4.785e-06 -0.543 0.5875
## CHRPH 2.489e-07 7.983e-06 0.031 0.9751
## CHRAC 2.195e-06 1.125e-06 1.951 0.0519 .
## CHRPV -5.694e-06 7.405e-06 -0.769 0.4425
## Total_charges_fossiles -4.788e-06 4.224e-06 -1.134 0.2578
## ELECT -5.226e-06 1.120e-05 -0.467 0.6410
## LANMX 8.873e-06 1.178e-04 0.075 0.9400
## Intra_conso_animaux -3.797e-06 2.459e-06 -1.544 0.1236
## SAUTI_classe50-99 -3.840e-02 5.812e-02 -0.661 0.5093
## SAUTI_classe100-149 4.381e-03 8.433e-02 0.052 0.9586
## SAUTI_classe150-199 2.803e-02 1.119e-01 0.251 0.8023
## SAUTI_classe200-249 7.561e-02 1.719e-01 0.440 0.6603
## SAUTI_classe250-299 7.685e-02 1.894e-01 0.406 0.6852
## SAUTI_classe300-349 1.458e-01 2.096e-01 0.696 0.4870
## SAUTI_classe350-399 4.325e-02 2.585e-01 0.167 0.8673
## SAUTI_classe400-449 1.282e-01 2.904e-01 0.441 0.6592
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3319 on 323 degrees of freedom
## Multiple R-squared: 0.1264, Adjusted R-squared: 0.004674
## F-statistic: 1.038 on 45 and 323 DF, p-value: 0.4111La p value du test globale est supérieur à 0.05 donc notre modèle peut être saturé, nous allons procéder à une sélection de variable.
## Start: AIC=-771.17
## ratio_Achat ~ PBUCE + OTE64F + REGIO + SFP + SUT1GCU + CHRGTO +
## VAPBR + EBEXP + UGBBO_classe + CHREN + CHRPH + CHRAC + CHRPV +
## Total_charges_fossiles + ELECT + LANMX + Intra_conso_animaux +
## SAUTI_classe + CHRAC`` Le modèle proposé contient encore beaucoup de variables donc allons donc procéder à une sélection manuelle avec les p value pour avoir un modèle explicatif.
##
## Call:
## lm(formula = ratio_Achat ~ PBUCE + VAPBR + CHREN + Intra_conso_animaux,
## data = data_mod_petits_ruminants)
##
## Residuals:
## Min 1Q Median 3Q Max
## -5.3275 -0.0859 -0.0071 0.0776 1.5670
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.646e-01 2.772e-02 9.544 < 2e-16 ***
## PBUCE 6.599e-07 2.520e-07 2.619 0.00919 **
## VAPBR -8.274e-07 3.647e-07 -2.268 0.02389 *
## CHREN -5.137e-06 3.127e-06 -1.643 0.10125
## Intra_conso_animaux -4.928e-06 1.756e-06 -2.806 0.00529 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3268 on 364 degrees of freedom
## Multiple R-squared: 0.04527, Adjusted R-squared: 0.03478
## F-statistic: 4.315 on 4 and 364 DF, p-value: 0.002012Le modèle retenu : ratio_Achat ~ PBUCE + VAPBR + CHREN + Intra_conso_animaux
## 6 117 135 147 235
## 3 112 130 142 228## Anova Table (Type II tests)
##
## Response: ratio_Achat
## Sum Sq Df F value Pr(>F)
## PBUCE 0.732 1 6.8584 0.009191 **
## VAPBR 0.550 1 5.1459 0.023887 *
## CHREN 0.288 1 2.6994 0.101249
## Intra_conso_animaux 0.841 1 7.8734 0.005286 **
## Residuals 38.875 364
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## (Intercept) PBUCE VAPBR CHREN
## 2.645746e-01 6.598604e-07 -8.273773e-07 -5.137315e-06
## Intra_conso_animaux
## -4.927822e-06
Toutes nos variables explicatives sont en euros nous pouvons donc
comparer leurs coefficients entre eux. En valeur absolue, il s’agit de
la variable charge d’engrais (CHREN) qui semble la plus influente sur le
ratio d’achat. Plus l’agriculteur a des charges d’engrais plus le ratio
diminue en effet s’il a des charges d’engrais importants il a souvent
beaucoup de terres et donc moins besoin d’acheter à l’extérieur. PBUCE
est la variable qui lorsque le production augmente le ratio aussi. Nous
retrouvons les résultats que l’intra consommation diminue lorsque le
ratio achat augmente.
Les hypothèses sont également vérifiées. #####MANQUE AUTONOMIE#####
Modèle linéaire sur la classe Granivores :
Comme il y a avait une grande classe granivores, on va tester le modèle linéaire sur les élevages spécialisées porcins (engraissement, naisseur et les deux), ayant plus de 20UGBPO.
Les hypothèses sont globalement vérifiées, donc nous allons affiner le
modèle à l’aide d’un step() (critère utilisé : AIC car on souhaite un
modèle plutôt explicatif). On obtient le modèle suivant (à la suite de
test des facteurs pour affiner davantage à la main) :
## Start: AIC=4216.58
## ACHAC ~ OTE64F + UGBTO_classe + REGIONS + SFP + SUT1GCU + CHRGTO +
## PBRTO + VAPBR + EBEXP + CHREN + CHRPH + CHRPV + Total_charges_fossiles +
## ELECT + PITOT + INTVP + Intra_conso_animaux + PBUCE + UTANS +
## UTATO## Analysis of Variance Table
##
## Response: ACHAC
## Df Sum Sq Mean Sq F value Pr(>F)
## OTE64F 2 1.8087e+12 9.0433e+11 45.804 < 2.2e-16 ***
## REGIONS 10 2.6559e+12 2.6559e+11 13.452 < 2.2e-16 ***
## SUT1GCU 1 6.1463e+11 6.1463e+11 31.131 1.000e-07 ***
## CHRGTO 1 2.7746e+13 2.7746e+13 1405.331 < 2.2e-16 ***
## PBRTO 1 2.5218e+11 2.5218e+11 12.773 0.0004644 ***
## Total_charges_fossiles 1 3.9847e+11 3.9847e+11 20.182 1.338e-05 ***
## Residuals 161 3.1787e+12 1.9743e+10
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Call:
## lm(formula = ACHAC ~ OTE64F + REGIONS + SUT1GCU + CHRGTO + PBRTO +
## Total_charges_fossiles, data = data_mod_porcins)
##
## Residuals:
## Min 1Q Median 3Q Max
## -589368 -49445 9312 67971 368001
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.165e+05 6.406e+04 -1.819 0.070831 .
## OTE64F5120 8.750e+04 5.705e+04 1.534 0.127055
## OTE64F5130 1.062e+05 5.653e+04 1.878 0.062203 .
## REGIONSBourgogne-Franche-Comté 8.610e+04 7.414e+04 1.161 0.247238
## REGIONSBretagne 9.042e+04 4.865e+04 1.859 0.064888 .
## REGIONSCentre-Val de Loire -1.447e+03 7.317e+04 -0.020 0.984242
## REGIONSGrand Est -1.595e+04 7.496e+04 -0.213 0.831808
## REGIONSHauts-de-France 7.471e+04 6.953e+04 1.075 0.284200
## REGIONSNormandie 6.398e+04 6.312e+04 1.014 0.312307
## REGIONSNouvelle-Aquitaine 1.287e+04 5.556e+04 0.232 0.817072
## REGIONSOccitanie 2.010e+04 7.475e+04 0.269 0.788346
## REGIONSPACA 7.118e+04 1.484e+05 0.480 0.632045
## REGIONSPays de la Loire -1.864e+04 5.753e+04 -0.324 0.746395
## SUT1GCU -6.728e+02 2.578e+02 -2.610 0.009917 **
## CHRGTO 1.008e+00 9.732e-02 10.362 < 2e-16 ***
## PBRTO -2.917e-01 8.289e-02 -3.519 0.000563 ***
## Total_charges_fossiles -5.963e+00 1.327e+00 -4.492 1.34e-05 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 140500 on 161 degrees of freedom
## Multiple R-squared: 0.9133, Adjusted R-squared: 0.9047
## F-statistic: 106 on 16 and 161 DF, p-value: < 2.2e-16## Warning: les observations ayant un 'leverage = 0' ne sont pas représentées sur le graphique :
## 178
On observe que toutes les variables présentes
(OTE64F+REGIONS+SUT1GCU+CHRGTO+PBRTO+Total_charges_fossiles) ont un
impact significatif sur les achats d’aliments concentrés. Concernant
l’impact des différentes variables sur les achats d’aliments concentrés,
on observe que :
* Pour les types d’exploitations, les exploitations 5130 (combinant
l’élevage et l’engraissement de porcins) ont des achats plus élevés que
les 5110 (pris pour référence, qui sint spécialisées porcins d’élevage),
et que les 5120 (spécialisées porcins engraissement, qui n’ont elles pas
des achats différents des 5110).
* Pour les régions, seule la Bretagne ressort significativement
différente de la région AURA prise pour référence, (avec, comme on
pouvait s’y attendre, un effet très positif sur les achats d’aliments).
* Pour les variables quantitatives, on observe que les surfaces de
culture (SUT1GCU) ont un impact négatif sur les achats de concentrés :
c’est tout à fait logique ; plus on a de surfaces en culture, moins on a
besoin d’en acheter à l’extérieur. Le produit total (incluant les
subventions, PBRTO) a un impact très faible (mais significatif, et
négatif) sur les achats extérieurs, tout comme le total des charges
fossiles.
On va réaliser les mêmes modèles, mais pour l’autonomie : on prend comme variable à expliquer le ratio “Achats d’aliments / (Produit total - Subventions + Intraconsommation totale)”, c’est à dire que l’on divise les achats d’aliments par ce qui a été vendu (moins les subventions) et ce qui a été intraconsommé.
Les hypothèses sont globalement vérifiées, avec un léger problème
d’homoscedasticité. On l’a gérer en enlevant des outliers après
l’affinage du modèle. Nous allons affiner le modèle à l’aide d’un step()
(critère utilisé : AIC car on souhaite un modèle plutôt explicatif). On
obtient le modèle suivant (à la suite de test des facteurs pour affiner
davantage à la main) :
## Start: AIC=-636.5
## ratio_Achat ~ OTE64F + UGBTO_classe + REGIONS + SFP + SUT1GCU +
## CHRGTO + PBRTO + VAPBR + EBEXP + CHREN + CHRPH + CHRPV +
## Total_charges_fossiles + ELECT + PITOT + INTVP + Intra_conso_animaux +
## PBUCE + UTANS + UTATO## Analysis of Variance Table
##
## Response: ratio_Achat
## Df Sum Sq Mean Sq F value Pr(>F)
## UGBTO_classe 10 0.6029 0.06029 2.4049 0.0108705 *
## SUT1GCU 1 1.0422 1.04218 41.5684 1.25e-09 ***
## CHRGTO 1 0.2964 0.29639 11.8217 0.0007443 ***
## PBRTO 1 0.3054 0.30541 12.1815 0.0006219 ***
## VAPBR 1 0.4336 0.43357 17.2935 5.18e-05 ***
## EBEXP 1 0.3058 0.30577 12.1960 0.0006174 ***
## Residuals 162 4.0616 0.02507
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Anova Table (Type II tests)
##
## Response: ratio_Achat
## Sum Sq Df F value Pr(>F)
## UGBTO_classe 0.4956 10 1.9765 0.0389756 *
## SUT1GCU 0.6839 1 27.2778 5.349e-07 ***
## CHRGTO 0.3701 1 14.7606 0.0001749 ***
## PBRTO 0.1880 1 7.4968 0.0068726 **
## VAPBR 0.7231 1 28.8429 2.685e-07 ***
## EBEXP 0.3058 1 12.1960 0.0006174 ***
## Residuals 4.0616 162
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Call:
## lm(formula = ratio_Achat ~ UGBTO_classe + SUT1GCU + CHRGTO +
## PBRTO + VAPBR + EBEXP, data = data_mod_porcins)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.48148 -0.06265 0.02140 0.09693 0.35122
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.290e-01 9.154e-02 2.502 0.013348 *
## UGBTO_classe50-99 5.532e-02 1.293e-01 0.428 0.669338
## UGBTO_classe100-149 3.712e-01 1.211e-01 3.067 0.002538 **
## UGBTO_classe150-199 2.224e-01 1.033e-01 2.154 0.032755 *
## UGBTO_classe200-249 2.640e-01 1.074e-01 2.459 0.015002 *
## UGBTO_classe250-299 2.122e-01 1.034e-01 2.053 0.041681 *
## UGBTO_classe300-349 3.159e-01 1.008e-01 3.135 0.002042 **
## UGBTO_classe350-399 2.833e-01 1.038e-01 2.729 0.007064 **
## UGBTO_classe400-449 1.968e-01 1.064e-01 1.850 0.066204 .
## UGBTO_classe450-499 1.949e-01 1.101e-01 1.770 0.078556 .
## UGBTO_classe500+ 2.584e-01 9.814e-02 2.633 0.009284 **
## SUT1GCU -1.176e-03 2.251e-04 -5.223 5.35e-07 ***
## CHRGTO 5.930e-07 1.544e-07 3.842 0.000175 ***
## PBRTO -4.150e-07 1.516e-07 -2.738 0.006873 **
## VAPBR -9.798e-07 1.824e-07 -5.371 2.69e-07 ***
## EBEXP 8.433e-07 2.415e-07 3.492 0.000617 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.1583 on 162 degrees of freedom
## Multiple R-squared: 0.4237, Adjusted R-squared: 0.3704
## F-statistic: 7.941 on 15 and 162 DF, p-value: 3.136e-13
Les variables restantes ont des impacts significatifs sur le ratio
d’autonomie alimentaire calculé. détailler un peu plus.
##Etude des charges énergétiques fossiles##
Le total des charges fossiles correspond à la somme des charges des combustibles de l’exploitation.
