Les tests de Shapiro–Wilk ont indiqué que les distributions de toutes les variables mesurées ne s'écartaient pas de manière significative de la normalité au cours des deux saisons (p > 0,05). Une analyse de la variance a révélé que l'interaction entre le niveau d'infusion de compost (Ct) et le traitement foliaire au biostimulant (F) affectait de manière significative les caractéristiques de croissance végétative, le nombre d'ombelles et les caractéristiques liées à l'huile (pourcentage d'huile, rendement en huile par plante et rendement en huile par hectare) au cours des deux saisons de croissance.
Réponse en rendement par rapport au traitement témoin
Lors de la première saison (Figure 2), le rendement en huile par hectare (Oy ha⁻1) a augmenté progressivement avec le niveau croissant de thé de compost. Tous les traitements Ct3 (600 L ha⁻1) ont produit des rendements en huile supérieurs à celui du témoin (T1 : Ct1 × F1, 360 L ha⁻1 sans pulvérisation foliaire). Le rendement en huile le plus élevé (55,58 kg ha⁻1) a été obtenu avec le traitement combinant Ct3 et Spirulina à 2 cm3 L⁻1 (T13), dépassant significativement le témoin et tous les traitements Ct1, avec une différence supérieure à la valeur LSD (2,93 kg ha⁻1). Par rapport au témoin, ce traitement a montré une augmentation marquée du rendement en huile par hectare (Figure 3).

Figure 3 : Analyse en composantes principales sous forme de biplot pour les quinze traitements en interaction, basée sur les quatre caractères végétatifs, le nombre d'ombelles par plante et les trois caractères liés à l'huile, en utilisant les deux premières composantes ; chaque point représente un traitement. La position à l'intérieur du graphique indique le degré de similitude ou de différence entre les traitements selon les facteurs d'entrée (8 caractéristiques de croissance et d'huile) ; plus les points sont proches, plus les traitements sont similaires, et plus ils sont éloignés, plus ils sont distincts. PH : hauteur de la plante (cm), NB : nombre de branches/plante, HFW : masse fraîche totale de la plante (g), HDW : masse sèche totale de la plante (g), O% : teneur en huile essentielle, Oy/p : rendement en huile par plante (cm3), Oy/ha : rendement en huile par hectare (L) et NU : nombre d'ombelles/plante T1 (Ct1 × F1), T2 (Ct1 × F2), T3 (Ct1 × F3), T4 (Ct1 × F4), T5 (Ct1 × F5), T6 (Ct2 × F1), T7 (Ct2 × F2), T8 (Ct2 × F3), T9 (Ct2 × F4), T10 (Ct2 × F5), T11 (Ct3 × F1), T12 (Ct3 × F2), T13 (Ct3 × F3), T14 (Ct3 × F4), et T15 (Ct3 × F5). Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
Des tendances similaires ont été observées pour le rendement en huile par plante et le nombre d'ombelles, indiquant que les traitements produisant un rendement en huile plus élevé par hectare étaient également associés à une production reproductive plus importante et à une accumulation accrue de biomasse. Ces relations sont présentées comme des associations statistiques entre les caractères mesurés, sans inférer de mécanismes causaux.
Au cours de la deuxième saison, le traitement Ct3 × F3 (T13) a de nouveau produit le rendement en huile le plus élevé par hectare (58,85 kg ha⁻1), dépassant significativement le témoin et la plupart des traitements Ct2 (LSD = 3,17 kg ha⁻1). Le classement des traitements par rapport au témoin était largement cohérent d'une saison à l'autre, indiquant des effets reproductibles des traitements dans les conditions environnementales testées.
Analyse en composantes principales (PCA)
L'analyse en composantes principales a été appliquée comme outil exploratoire afin de résumer les relations multivariées entre les huit caractères mesurés dans les 15 traitements. Les deux premières composantes principales expliquent 98,4 % de la variance totale, la CP₁ représentant 96,5 % et la CP₂ 1,9 % (Tableau 2 ; Figure 4). La prédominance de la CP₁ reflète une forte covariance entre les caractères liés à la croissance et au rendement dans le jeu de données, plutôt que des processus biologiques sous-jacents distincts.
| Paramètre | CP1 | CP2 | CP3 | CP4 |
| NU | 0.3516 | 0.3064 | -0.3732 | 0.772 |
| O% | 0.349 | -0.6068 | -0.1873 | -0.0634 |
| Oy/p | 0.3591 | -0.11 | -0.0028 | -0.0974 |
| Oy/ha | 0.3583 | -0.1579 | -0.1046 | -0.2021 |
| PH | 0.3476 | 0.5757 | 0.3834 | -0.3296 |
| BN | 0.3495 | -0.26 | 0.7504 | 0.3693 |
| PMF | 0.3551 | 0.3138 | -0.1821 | -0.1535 |
| PMS | 0.358 | -0.0575 | -0.2675 | -0.2847 |
| Valeur propre | 7.7201 | 0.1518 | 0.0721 | 0.0346 |
| % de variance | 96,5 | 1,9 | 0,9 | 0,43 |
| Cumulatif (%) | 96,5 | 98,4 | 99,3 | 99,73 |
| PH : hauteur de la plante (cm), BN : nombre de branches/plante, PMF : masse fraîche totale de la plante (g), PMS : masse sèche totale de la plante (g), O% : teneur en huile essentielle, Oy/p : rendement en huile par plante (cm³), Oy/ha : rendement en huile par hectare (L) et NU : nombre d'ombelles/plante |
Tableau 2 : Valeurs propres, vecteurs propres, pourcentage de variance et variance cumulée des quatre premières composantes principales basés sur huit caractères liés à la croissance et à la production d'huile, pour 15 combinaisons de traitements. Les abréviations des caractères sont les suivantes : PH (hauteur de la plante, cm), NB (nombre de branches par plante), HFW (masse fraîche de la plante, g), HDW (masse sèche de la plante, g), O% (teneur en huile essentielle), Oy/p (rendement en huile par plante, cm3), Oy/ha (rendement en huile par hectare, L) et NU (nombre d'ombelles par plante).

Figure 4 : Graphique des chargements du vecteur propre pour les 8 caractères étudiés soumis aux 15 traitements combinés : les barres de CP₁ utilisent l'axe gauche ; CP₂, CP₃, CP₄ utilisent des formes et des couleurs de points distinctes sur l'axe droit. PH : hauteur de la plante (cm), NB : nombre de branches/plante, HFW : poids frais total de la plante (g), HDW : poids sec total de la plante (g), O% : teneur en huile essentielle, Oy/p : rendement en huile par plante (cm3), Oy/ha : rendement en huile par hectare (L) et NU : nombre d'ombelles/plante. Les barres d'erreur représentent les erreurs types de la moyenne (ETM). Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
Tous les caractères présentaient de fortes charges positives similaires sur la CP₁, notamment le nombre d'ombelles, le pourcentage d'huile, le rendement en huile par plante, le rendement en huile par hectare, la hauteur de la plante, le nombre de branches, le poids frais de la plante et le poids sec de la plante (Tableau 2 ; Figure 5). Par conséquent, la CP₁ peut être interprétée comme un gradient général de productivité résumant la variation des caractères liés à la biomasse et au rendement en huile entre les traitements. La CP₂ expliquait seulement une faible proportion de la variance et apportait une discrimination supplémentaire limitée entre les traitements.

Figure 5 : Graphique des scores de l'ACP illustrant la répartition des 15 combinaisons de traitements en fonction des caractères mesurés. Les traitements sont définis comme suit : T1 (Ct1 × F1), T2 (Ct1 × F2), T3 (Ct1 × F3), T4 (Ct1 × F4), T5 (Ct1 × F5), T6 (Ct2 × F1), T7 (Ct2 × F2), T8 (Ct2 × F3), T9 (Ct2 × F4), T10 (Ct2 × F5), T11 (Ct3 × F1), T12 (Ct3 × F2), T13 (Ct3 × F3), T14 (Ct3 × F4) et T15 (Ct3 × F5). Les barres d'erreur représentent les erreurs standards de la moyenne (SEM). Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
Scores de l'ACP et comparaison avec le témoin
Les traitements présentant des scores PC₁ fortement positifs (notamment T13, T12, T15 et T14) ont également montré les augmentations en pourcentage les plus élevées par rapport au témoin (T1) pour la plupart des caractères (Figure 5 et Figure 6). Par exemple, T13 a présenté une augmentation de 94,44 % du nombre d'ombelles et de 98,18 % du rendement en huile par hectare par rapport au témoin.

Figure 6 : Augmentations relatives en pourcentage des caractères mesurés par rapport au traitement témoin (T1 : Ct1 × F1). Les traitements sont définis comme suit : T1 (Ct1 × F1), T2 (Ct1 × F2), T3 (Ct1 × F3), T4 (Ct1 × F4), T5 (Ct1 × F5), T6 (Ct2 × F1), T7 (Ct2 × F2), T8 (Ct2 × F3), T9 (Ct2 × F4), T10 (Ct2 × F5), T11 (Ct3 × F1), T12 (Ct3 × F2), T13 (Ct3 × F3), T14 (Ct3 × F4) et T15 (Ct3 × F5). Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
Les traitements ayant des scores PC₁ intermédiaires (par exemple, T8 et T10) ont montré des améliorations modérées par rapport au témoin, tandis que les traitements ayant des scores PC₁ négatifs, y compris le témoin et d'autres combinaisons Ct1, ont présenté des changements relatifs minimes. Ces profils indiquent que les scores PC₁ étaient étroitement associés aux réponses aux traitements ajustées par rapport à la ligne de base, confirmant ainsi les résultats de l'ANOVA univariée. Toutefois, l'ACP a été utilisée uniquement pour la réduction des données et la visualisation des profils multivariés, et ne remplace pas les tests statistiques inférentiels.
Associations entre caractères dans le biplot ACP
Le biplot de l'ACP (Figure 4) illustre des associations étroites entre le rendement en huile par plante, le rendement en huile par hectare et le poids sec des parties aériennes, reflétant de fortes corrélations entre ces variables. La hauteur de la plante, le poids frais des parties aériennes et le nombre d'ombelles présentaient également des orientations vectorielles similaires, indiquant une variation coordonnée entre ces caractères au sein du jeu de données. Le nombre de branches montrait une séparation partielle par rapport aux autres composantes du rendement, suggérant des associations plus faibles dans les conditions expérimentales actuelles. Ces relations décrivent uniquement des corrélations observées et n'impliquent pas de liens mécanistiques.
Les traitements (T12–T15) ayant reçu du thé de compost à raison de 600 L ha⁻1 combiné à des biostimulants foliaires ont montré une performance multivariée supérieure pour les caractéristiques de croissance et de teneur en huile. En revanche, les traitements (T7–T11) ont présenté des valeurs relativement plus faibles pour la plupart des paramètres évalués. Les traitements situés en direction opposée à la plupart des vecteurs, comme Ct1 (360 L ha⁻1) sans pulvérisation foliaire ou combiné à Nostoc, étaient associés à des valeurs plus faibles pour la plupart des caractères mesurés.
Dans l'ensemble, l'ACP a fourni un résumé multivarié concis qui complète l'ANOVA en illustrant en quoi les traitements différaient du témoin selon plusieurs caractères simultanément. La cohérence entre les profils de l'ACP et les résultats univariés ajustés au niveau de base soutient la robustesse des effets des traitements, sans étendre l'interprétation au-delà des variables mesurées.
DISPONIBILITÉ DES DONNÉES :
Les données brutes de cet article sont disponibles dans les fichiers supplémentaires.
Tableau supplémentaire 1 : Données expérimentales brutes des résultats. Veuillez cliquer ici pour télécharger ce fichier.