Testy Shapiro–Wilka wykazały, że rozkłady wszystkich mierzonych zmiennych nie odbiegały istotnie od rozkładu normalnego w żadnym z sezonów (p > 0,05). Analiza wariancji ujawniła, że interakcja między poziomem herbatki kompostowej (Ct) a traktowaniem biostymulatorem liściowym (F) istotnie wpływała na cechy wzrostu wegetatywnego, liczbę baldachów oraz cechy związane z olejkiem (procentową zawartość olejku, wydajność olejku z rośliny oraz wydajność olejku z hektara) w obu sezonach wegetacyjnych.
Wydajność w stosunku do grupy kontrolnej
W pierwszym sezonie (Rysunek 2) plon oleju z hektara (Oy ha⁻1) wzrastał stopniowo wraz ze zwiększeniem poziomu herbatki kompostowej. Wszystkie warianty Ct3 (600 L ha⁻1) pozwoliły uzyskać wyższe plony oleju niż w wariancie kontrolnym (T1: Ct1 × F1, 360 L ha⁻1 bez oprysku dolistnego). Najwyższy plon oleju (55,58 kg ha⁻1) uzyskano dla Ct3 w połączeniu ze Spirulina w stężeniu 2 cm3 L⁻1 (T13), co znacząco przekroczyło wartość kontrolną oraz wszystkie warianty Ct1, wykraczając poza wartość LSD (2,93 kg ha⁻1). W porównaniu z kontrolą, wariant ten wykazał wyraźny wzrost plonu oleju z hektara (Rysunek 3).

Rysunek 3Wykres biplot analizy głównych składowych dla piętnastu wchodzących w interakcję zabiegów, opracowany na podstawie czterech cech wegetatywnych, liczby baldachów na roślinę oraz trzech cech związanych z olejem, z wykorzystaniem pierwszej i drugiej składowej; każdy punkt reprezentuje pojedynczy zabieg. Położenie wewnątrz wykresu wskazuje stopień podobieństwa lub różnicy między wariantami w zależności od czynników wejściowych (8 cech wzrostu i właściwości oleju); im bliższe są one względem siebie, tym bardziej podobne, a im dalej od siebie, tym bardziej odmienne. PH: wysokość rośliny (cm), NB: liczba pędów/roślinę, HFW: całkowita świeża masa nadziemna (g), HDW: całkowita sucha masa nadziemna (g), O%: zawartość olejku eterycznego, Oy/p: wydajność olejku z rośliny (cm)3), Oy/ha: plon oleju z hektara (L) oraz NU: liczba baldachów na roślinę 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) oraz T15 (Ct3 × F5). Aby wyświetlić powiększoną wersję tej figury, kliknij tutaj.
Podobne trendy zaobserwowano w przypadku wydajności oleju z rośliny oraz liczby baldachów, co wskazuje na to, że zabiegi zapewniające wyższą wydajność oleju z hektara wiązały się również z większym potencjałem rozrodczym i akumulacją biomasy. Relacje te przedstawiono jako powiązania statystyczne między mierzonymi cechami, bez wyciągania wniosków na temat mechanizmów przyczynowo-skutkowych.
W drugim sezonie wariant Ct3 × F3 (T13) ponownie przyniósł najwyższą wydajność oleju z hektara (58,85 kg ha⁻1), znacząco przewyższając kontrolę oraz większość zabiegów Ct2 (LSD = 3,17 kg ha⁻1). Ranking zabiegów w stosunku do kontroli był w dużej mierze spójny w obu sezonach, co wskazuje na powtarzalność efektów zabiegów w badanych warunkach środowiskowych.
Analiza głównych składowych (PCA)
Analizę głównych składowych zastosowano jako narzędzie eksploracyjne do podsumowania zależności wielowymiarowych pomiędzy ośmioma pomierzonymi cechami w 15 wariantach doświadczalnych. Dwie pierwsze główne składowe wyjaśniły 98,4% całkowitej wariancji, przy czym PC₁ odpowiadała za 96,5%, a PC₂ za 1,9% (Tabela 2; Rysunek 4). Dominacja PC₁ odzwierciedla silną kowariancję cech związanych ze wzrostem i plonowaniem w zbiorze danych, a nie odrębne, podstawowe procesy biologiczne.
| Parametr | Komputer osobisty1 | PC2 | PC3 | Komputer osobisty4 |
| Jednostka Nukleotydowa | 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 |
| 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 |
| Wartość własna | 7.7201 | 0.1518 | 0.0721 | 0.0346 |
| % wariancji | 96.5 | 1.9 | 0.9 | 0.43 |
| Skumulowane (%) | 96.5 | 98.4 | 99.3 | 99.73 |
| PH: wysokość rośliny (cm), NB: liczba pędów na roślinę, HFW: całkowita świeża masa zielna (g), HDW: całkowita sucha masa zielna (g), O%: zawartość olejku eterycznego, Oy/p: wydajność olejku na roślinę (cm³), Oy/ha: wydajność olejku z hektara (L) oraz NU: liczba baldachów na roślinę |
Tabela 2: Wartości własne, wektory własne, procent wariancji i wariancja skumulowana pierwszych czterech głównych składowych na podstawie ośmiu cech związanych z wzrostem i zawartością oleju w 15 kombinacjach zabiegów. Skróty cech są następujące: PH (wysokość rośliny, cm), NB (liczba pędów na roślinę), HFW (świeża masa części nadziemnej, g), HDW (sucha masa części nadziemnej, g), O% (zawartość olejku eterycznego), Oy/p (wydajność oleju z rośliny, cm3), Oy/ha (wydajność oleju z hektara, L) oraz NU (liczba baldachów na roślinę).

Rysunek 4: Wykres obciążeń wektora własnego dla 8 badanych cech poddanych 15 interakcyjnym zabiegom: słupki PC₁ odnoszą się do lewej osi; PC₂, PC₃, PC₄ są przedstawione za pomocą różnych kształtów punktów i kolorów na prawej osi. PH: wysokość rośliny (cm), NB: liczba pędów/roślina, HFW: całkowita świeża masa zielna (g), HDW: całkowita sucha masa zielna (g), O%: zawartość olejku eterycznego, Oy/p: wydajność olejku na roślinę (cm3), Oy/ha: wydajność olejku z hektara (L) oraz NU: liczba baldachów/roślina. Słupki błędów reprezentują standardowe błędy średniej (SEM). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Wszystkie cechy wykazały silne i podobne dodatnie ładunki na PC₁, w tym liczba baldachów, zawartość oleju, wydajność oleju z rośliny, wydajność oleju z hektara, wysokość rośliny, liczba gałęzi, świeża masa zielnej części rośliny oraz sucha masa zielnej części rośliny (Tabela 2; Rysunek 5). W związku z tym PC₁ można interpretować jako ogólny gradient produktywności podsumowujący zmienność cech biomasy i wydajności oleju w poszczególnych wariantach. PC₂ wyjaśnił jedynie niewielką część wariancji i wniósł ograniczoną dodatkową dyskryminację między wariantami.

Rysunek 5: Wykres wyników PCA ilustrujący rozkład 15 kombinacji zabiegowych na podstawie mierzonych cech. Zabiegi zdefiniowano następująco: 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) oraz T15 (Ct3 × F5). Słupki błędów reprezentują błędy standardowe średniej (SEM). Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.
Wyniki PCA i porównanie z grupą kontrolną
Warianty z wysokimi dodatnimi wartościami PC₁ (w szczególności T13, T12, T15 i T14) wykazały również największe procentowe wzrosty w stosunku do kontroli (T1) w większości badanych cech (Rysunek 5 i Rysunek 6). Na przykład w przypadku T13 odnotowano wzrost o 94,44% w liczbie baldachów oraz o 98,18% w plonie oleju z hektara w porównaniu z kontrolą.

Rycina 6: Względne wzrosty procentowe mierzonych cech w porównaniu do zabiegu kontrolnego (T1: Ct1 × F1). Zabiegi zdefiniowano następująco: 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) oraz T15 (Ct3 × F5). Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Zabiegi z pośrednimi wynikami PC₁ (np. T8 i T10) wykazały umiarkowaną poprawę w stosunku do grupy kontrolnej, natomiast zabiegi z ujemnymi wynikami PC₁, w tym grupa kontrolna i inne kombinacje Ct1, wykazały minimalne zmiany względne. Wzorce te wskazują, że wyniki PC₁ były ściśle powiązane z odpowiedziami na zabiegi skorygowanymi o wartość bazową, co potwierdza wyniki jednoczynnikowej analizy wariancji (ANOVA). Niemniej jednak PCA została wykorzystana wyłącznie do redukcji danych i wizualizacji wzorców wielowymiarowych i nie zastępuje inferencyjnych testów statystycznych.
Powiązania cech na biwykresie PCA
Wykres biplot PCA (Rysunek 4) obrazuje bliskie powiązania między wydajnością oleju z rośliny, wydajnością oleju z hektara a suchą masą ziół, co odzwierciedla silne korelacje między tymi zmiennymi. Wysokość roślin, świeża masa ziół oraz liczba baldachów również wykazały podobne orientacje wektorów, co wskazuje na skoordynowaną zmienność tych cech w zbiorze danych. Liczba gałęzi wykazała częściową separację od pozostałych komponentów plonu, co sugeruje słabsze powiązania w obecnych warunkach eksperymentalnych. Relacje te opisują jedynie zaobserwowane korelacje i nie implikują powiązań mechanistycznych.
Warianty (T12–T15), w których zastosowano herbatkę kompostową w dawce 600 L ha⁻1 w połączeniu z biostymulatorami dolistnymi, wykazały wyższą wydajność w analizie wielowymiarowej w zakresie cech wzrostu i zawartości oleju. W przeciwieństwie do nich, warianty (T7–T11) wykazały porównawczo niższe wartości dla większości ocenianych parametrów. Warianty usytuowane w kierunku przeciwnym do większości wektorów, takie jak Ct1 (360 L ha⁻1) bez oprysku dolistnego lub w połączeniu z Nostoc, wiązały się z niższymi wartościami dla większości mierzonych cech.
Podsumowując, analiza PCA dostarczyła zwięzłego wielowymiarowego zestawienia, które uzupełnia analizę ANOVA, ilustrując różnice między grupami traktowanymi a grupą kontrolną w obrębie wielu cech jednocześnie. Spójność między wzorcami PCA a jednowymiarowymi wynikami skorygowanymi o wartość bazową potwierdza trwałość efektów traktowania bez rozszerzania interpretacji poza zmierzone zmienne.
DOSTĘPNOŚĆ DANYCH:
Surowe dane do niniejszego manuskryptu znajdują się w plikach uzupełniających.
Tabela uzupełniająca 1: Surowe dane eksperymentalne dla uzyskanych wyników. Kliknij tutaj, aby pobrać ten plik.