Shapiro-Wilk-toetsen gaven aan dat de verdelingen van alle gemeten variabelen in geen van beide seizoenen significant afweken van de normaliteit (p > 0.05). Variantieanalyse wees uit dat de interactie tussen het niveau van compostthee (Ct) en de behandeling met bladerbiostimulantia (F) de vegetatieve groeikenmerken, het aantal schermbloemen en de olie-gerelateerde kenmerken (oliepercentage, olieopbrengst per plant en olieopbrengst per hectare) in beide groeiseizoenen significant beïnvloedde.
Opbrengstrespons ten opzichte van de controlebehandeling
In het eerste seizoen (Figuur 2) nam de olieopbrengst per hectare (Oy ha⁻1) progressief toe bij een toenemend niveau van compostthee. Alle Ct3-behandelingen (600 L ha⁻1) produceerden hogere olieopbrengsten dan de controlebehandeling (T1: Ct1 × F1, 360 L ha⁻1 zonder bladbespuiting). De hoogste olieopbrengst (55,58 kg ha⁻1) werd verkregen uit Ct3 gecombineerd met Spirulina bij 2 cm3 L⁻1 (T13), wat significant hoger was dan de controle en alle Ct1-behandelingen, waarbij de LSD-waarde (2,93 kg ha⁻1) werd overschreden. Vergeleken met de controle vertoonde deze behandeling een duidelijke toename in de olieopbrengst per hectare (Figuur 3).

Figuur 3: Biplot hoofdcomponentanalyse voor de vijftien interactieve behandelingen, gebaseerd op de vier vegetatieve kenmerken, het aantal schermen per plant en de drie olie-gerelateerde kenmerken, gebruikmakend van de eerste en tweede componenten; een behandeling wordt vertegenwoordigd door elk punt. De locatie in de plot geeft de mate van overeenkomst of verschil tussen de behandelingen weer volgens de inputfactoren (8 groeikenmerken en oliekarakteristieken); hoe dichter bij elkaar, des te gelijker, en hoe verder uit elkaar, des te duidelijker het onderscheid. PH: plantlengte (cm), NB: aantal takken/plant, HFW: totaal versgewicht van het kruid (g), HDW: totaal drooggewicht van het kruid (g), O%: essentieel oliegehalte, Oy/p: olieopbrengst per plant (cm3), Oy/ha: olieopbrengst per hectare (L) en NU: aantal schermen/plant 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), and T15 (Ct3 × F5). Klik hier om een grotere versie van deze figuur te bekijken.
Vergelijkbare trends werden waargenomen voor de olieopbrengst per plant en het aantal schermbloemen, wat aangeeft dat behandelingen die een hogere olieopbrengst per hectare opleverden, ook geassocieerd waren met een grotere reproductieve output en biomassa-accumulatie. Deze relaties worden gepresenteerd als statistische associaties tussen gemeten kenmerken, zonder causale mechanismen af te leiden.
In het tweede seizoen produceerde Ct3 × F3 (T13) opnieuw de hoogste olieopbrengst per hectare (58,85 kg ha⁻1), wat significant hoger was dan de controle en de meeste Ct2-behandelingen (LSD = 3,17 kg ha⁻1). De rangschikking van de behandelingen ten opzichte van de controle was grotendeels consistent over de seizoenen heen, wat wijst op reproduceerbare effecten van de behandelingen onder de geteste omgevingsomstandigheden.
Principecomponentenanalyse (PCA)
Principal component analyse werd toegepast als exploratief instrument om de multivariate relaties tussen de acht gemeten kenmerken over de 15 behandelingen samen te vatten. De eerste twee hoofdcomponenten verklaarden 98,4% van de totale variantie, waarbij PC₁ verantwoordelijk was voor 96,5% en PC₂ voor 1,9% (Tabel 2; Figuur 4). De dominantie van PC₁ weerspiegelt een sterke covariantie tussen groeigerelateerde en opbrengstgerelateerde kenmerken in de dataset, in plaats van onderscheidende onderliggende biologische processen.
| Parameter | PC1 | PC2 | PC3 | PC4 |
| NU | 0.3516 | 0.3064 | -0.3732 | 0.772 |
| 0% | 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 |
| NB | 0.3495 | -0.26 | 0.7504 | 0.3693 |
| HFW | 0.3551 | 0.3138 | -0.1821 | -0.1535 |
| HDW | 0.358 | -0.0575 | -0.2675 | -0.2847 |
| Eigenwaarde | 7.7201 | 0.1518 | 0.0721 | 0.0346 |
| % van de variantie | 96.5 | 1.9 | 0.9 | 0.43 |
| Cumulatief (%) | 96.5 | 98.4 | 99.3 | 99.73 |
| PH: plantelhoogte (cm), NB: aantal takken/plant, HFW: totaal vers gewicht van de kruiden (g), HDW: totaal drooggewicht van de kruiden (g), O%: gehalte aan essentiële olie, Oy/p: olieopbrengst per plant (cm³), Oy/ha: olieopbrengst per hectare (L) en NU: aantal schermkransen/plant |
Tabel 2: Eigenwaarden, eigenvectoren, percentage variantie en cumulatieve variantie van de eerste vier hoofdcomponenten op basis van acht groeigerelateerde en olie-gerelateerde kenmerken over 15 behandelcombinaties.Afkortingen voor kenmerken zijn als volgt: PH (planthoogte, cm), NB (aantal takken per plant), HFW (vers gewicht van het kruid, g), HDW (droog gewicht van het kruid, g), O% (essentiële oliegehalte), Oy/p (olieopbrengst per plant, cm3), Oy/ha (olieopbrengst per hectare, L) en NU (aantal schermen per plant).

Figuur 4: Loading plot van de eigenvectoren voor 8 bestudeerde kenmerken behandeld met de 15 interactieve behandelingen: PC₁ balken gebruiken de linkeras; PC₂, PC₃, PC₄ gebruiken verschillende puntvormen en kleuren op de rechteras. PH: planthoogte (cm), NB: aantal takken/plant, HFW: totaal vers gewicht kruid (g), HDW: totaal drooggewicht kruid (g), O%: gehalte aan essentiële olie, Oy/p: olieopbrengst per plant (cm3), Oy/ha: olieopbrengst per hectare (L) en NU: aantal schermen/plant. Foutenbalken vertegenwoordigen de standaardfouten van het gemiddelde (SEM). Klik hier om een grotere versie van deze figuur te bekijken.
Alle kenmerken vertoonden sterke en vergelijkbare positieve ladingen op PC₁, waaronder het aantal schermtoppen, het oliepercentage, de olieopbrengst per plant, de olieopbrengst per hectare, de plantlengte, het aantal takken, het natte gewicht van de kruiden en het droge gewicht van de kruiden (Tabel 2; Figuur 5). Bijgevolg kan PC₁ worden geïnterpreteerd als een algemene productiviteitsgradiënt die de variatie in biomassa- en olieopbrengstkenmerken over de behandelingen samenvat. PC₂ verklaarde slechts een klein deel van de variantie en droeg beperkt bij aan aanvullende discriminatie tussen de behandelingen.

Figuur 5PCA-scoreplot die de distributie van 15 behandelingscombinaties illustreert op basis van de gemeten kenmerken. De behandelingen zijn als volgt gedefinieerd: 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) en T15 (Ct3 × F5). Foutenbalken representeren de standaardfouten van het gemiddelde (SEM). Klik hier om een grotere versie van deze figuur te bekijken.
PCA-scores en vergelijking met de controle
Behandelingen met hoge positieve PC₁-scores (met name T13, T12, T15 en T14) vertoonden ook de grootste procentuele toenames ten opzichte van de controle (T1) voor de meeste kenmerken (Figuur 5 en Figuur 6). Zo vertoonde T13 toenames van 94,44% in het aantal schermbloemen en 98,18% in de olieopbrengst per hectare vergeleken met de controle.

Figuur 6: Relatieve percentage-toenames van de gemeten kenmerken vergeleken met de controlebehandeling (T1: Ct1 × F1). Behandelingen zijn als volgt gedefinieerd: 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) en T15 (Ct3 × F5). Klik hier om een grotere versie van deze figuur te bekijken.
Behandelingen met gemiddelde PC₁-scores (bijv. T8 en T10) vertoonden matige verbeteringen ten opzichte van de controle, terwijl behandelingen met negatieve PC₁-scores, inclusief de controle en andere Ct1-combinaties, minimale relatieve veranderingen vertoonden. Deze patronen geven aan dat de PC₁-scores nauw aansloten bij de op de baseline gecorrigeerde behandelingsresponsen, wat de resultaten van de univariate ANOVA bevestigde. PCA werd echter strikt gebruikt voor datareductie en visualisatie van multivariate patronen en vervangt geen inferentiële statistische toetsing.
Kenmerkassociaties in de PCA-biplot
De PCA-biplot (Figuur 4) illustreert nauwe associaties tussen de olieopbrengst per plant, de olieopbrengst per hectare en het drooggewicht van het kruid, wat wijst op sterke correlaties tussen deze variabelen. De plantlengte, het versgewicht van het kruid en het aantal schermvijzels vertoonden ook vergelijkbare vectororiëntaties, wat duidt op een gecoördineerde variatie tussen deze kenmerken binnen de dataset. Het aantal takken vertoonde een gedeeltelijke scheiding van andere opbrengstcomponenten, wat suggereert dat de associaties onder de huidige experimentele omstandigheden zwakker zijn. Deze relaties beschrijven uitsluitend de waargenomen correlaties en impliceren geen mechanische verbanden.
Behandelingen (T12–T15) die compostthee ontvingen bij 600 L ha⁻1 in combinatie met bladbiostimulanten vertoonden een superieure multivariate prestatie voor zowel groeikenmerken als olie-eigenschappen. In contrast hiermee vertoonden behandelingen (T7–T11) relatief lagere waarden voor de meeste geëvalueerde parameters. Behandelingen die in de tegenovergestelde richting van de meeste vectoren waren gepositioneerd, zoals Ct1 (360 L ha⁻1) zonder bladbespuiting of in combinatie met Nostoc, waren geassocieerd met lagere waarden voor de meeste gemeten kenmerken.
Over het geheel genomen bood de PCA een beknopte multivariate samenvatting die de ANOVA aanvult door te illustreren hoe de behandelingen gelijktijdig over meerdere kenmerken verschilden van de controle. De consistentie tussen de PCA-patronen en de op de baseline gecorrigeerde univariate resultaten ondersteunt de robuustheid van de behandelingseffecten, zonder de interpretatie uit te breiden voorbij de gemeten variabelen.
DATABESCHIKBAARHEID:
De ruwe gegevens voor dit manuscript zijn beschikbaar in de aanvullende bestanden.
Aanvullende tabel 1: Ruwe experimentele gegevens voor de resultaten. Klik hier om dit bestand te downloaden.