Shapiro–Wilk testleri, ölçülen tüm değişkenlerin dağılımlarının her iki sezonda da normallikten önemli ölçüde sapmadığını göstermiştir (p > 0.05). Varyans analizi, kompost çayı düzeyi (Ct) ile yaprak biostimülan uygulaması (F) arasındaki etkileşimin, her iki yetiştirme sezonunda da vejetatif büyüme özelliklerini, şemsiye sayısını ve yağ ile ilgili özellikleri (yağ yüzdesi, bitki başına yağ verimi ve hektar başına yağ verimi) önemli ölçüde etkilediğini ortaya koymuştur.
Kontrol uygulamasına göre verim yanıtı
İlk sezonda (Şekil 2), hektar başına yağ verimi (Oy ha⁻1), artan kompost çayı seviyesiyle birlikte kademeli olarak artmıştır. Tüm Ct3 (60 L ha⁻1) uygulamaları, kontrol uygulamasından (T1: yaprak püskürtmesiz Ct1 × F1, 360 L ha⁻1) daha yüksek yağ verimi sağlamıştır. En yüksek yağ verimi (5,58 kg ha⁻1), 2 cm3 L⁻1 dozunda Spirulina ile kombine edilen Ct3 (T13) uygulamasından elde edilmiş olup; bu değer, LSD değerini (2,93 kg ha⁻1) aşarak kontrolü ve tüm Ct1 uygulamalarını anlamlı şekilde geride bırakmıştır. Kontrole kıyasla, bu uygulama hektar başına yağ veriminde belirgin bir artış göstermiştir (Şekil 3).

Şekil 3Birinci ve ikinci bileşenler kullanılarak; dört vejetatif özellik, bitki başına şemsiye sayısı ve üç yağla ilgili özellik temelinde, etkileşimli on beş uygulama için biplot temel bileşen analizi; her bir nokta bir uygulamayı temsil etmektedir. Grafik içindeki konumlar, giriş faktörlerine (8 büyüme ve yağ özelliği) göre işlemler arasındaki benzerlik veya farklılık derecesini göstermektedir; noktalar birbirine ne kadar yakınsa o kadar benzer, ne kadar uzaksa o kadar belirgin şekilde farklıdırlar. PH: bitki boyu (cm), NB: bitki başına dal sayısı, HFW: toplam taze bitki ağırlığı (g), HDW: toplam kuru bitki ağırlığı (g), O%: uçucu yağ içeriği, Oy/p: bitki başına yağ verimi (cm3), Oy/ha: hektar başına yağ verimi (L) ve NU: bitki başına şemsiye sayısı 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), T1 (Ct3 × F1), T12 (Ct3 × F2), T13 (Ct3 × F3), T14 (Ct3 × F4) ve T15 (Ct3 × F5). Lütfen bu şeklin daha büyük bir versiyonunu görüntülemek için buraya tıklayın.
Bitki başına yağ verimi ve şemsiye sayısı için de benzer eğilimler gözlemlenmiş olup, hektar başına daha yüksek yağ verimi sağlayan uygulamaların daha fazla üreme çıktısı ve biyokütle birikimi ile ilişkili olduğu belirlenmiştir. Bu ilişkiler, nedensel mekanizmalar çıkarılmaksızın, ölçülen özellikler arasındaki istatistiksel ilişkiler olarak sunulmuştur.
İkinci sezonda, Ct3 × F3 (T13) hektar başına yeniden en yüksek yağ verimini (58,85 kg ha⁻1) vermiş olup, kontrol grubunu ve çoğu Ct2 işlemini önemli ölçüde aşmıştır (LSD = 3,17 kg ha⁻1). İşlemlerin kontrol grubuna göre sıralaması mevsimler boyunca büyük ölçüde tutarlı seyretmiş olup, bu durum test edilen çevresel koşullar altında tekrarlanabilir işlem etkilerine işaret etmektedir.
Temel bileşenler analizi (PCA)
15 uygulama genelinde ölçülen sekiz özellik arasındaki çok değişkenli ilişkileri özetlemek için keşifsel bir araç olarak temel bileşenler analizi uygulanmıştır. İlk iki temel bileşen, toplam varyansın %98,4'ünü açıklamış olup, PC₁ %96,5'ini ve PC₂ %1,9'unu oluşturmuştur (Tablo 2; Şekil 4). PC₁'in baskınlığı, veri setindeki farklı temel biyolojik süreçlerden ziyade, büyüme ve verimle ilgili özellikler arasındaki güçlü kovaryansı yansıtmaktadır.
| Parametre | PC1 | PC2 | PC3 | PC4 |
| 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 |
| 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 |
| Özdeğer | 7.7201 | 0.1518 | 0.0721 | 0.0346 |
| Varyans Yüzdesi | 96.5 | 1.9 | 0.9 | 0.43 |
| Kümülatif (%) | 96.5 | 98.4 | 99.3 | 99.73 |
| PH: bitki boyu (cm), NB: bitki başına dal sayısı, HFW: toplam taze bitki ağırlığı (g), HDW: toplam kuru bitki ağırlığı (g), O%: uçucu yağ içeriği, Oy/p: bitki başına yağ verimi (cm³), Oy/ha: hektar başına yağ verimi (L) ve NU: bitki başına şemsiye sayısı |
Tablo 2: 15 işlem kombinasyonu genelinde sekiz büyüme ve yağla ilgili özelliğe dayalı ilk dört temel bileşenin özdeğerleri, özvektörleri, varyans yüzdesi ve kümülatif varyansı. Özellik kısaltmaları şöyledir: PH (bitki boyu, cm), NB (bitki başına dal sayısı), HFW (bitki taze ağırlığı, g), HDW (bitki kuru ağırlığı, g), O% (uçucu yağ içeriği), Oy/p (bitki başına yağ verimi, cm3), Oy/ha (hektar başına yağ verimi, L) ve NU (bitki başına şemsiye sayısı).

Şekil 415 etkileşimli uygulamaya tabi tutulan 8 incelenen özellik için özvektör yükleme grafiği: PC₁ çubukları sol ekseni; PC₂, PC₃, PC₄ ise sağ eksen üzerinde farklı nokta şekillerini ve renklerini kullanır. PH: bitki boyu (cm), NB: bitki başına dal sayısı, HFW: toplam otsu taze ağırlık (g), HDW: toplam otsu kuru ağırlık (g), O%: uçucu yağ içeriği, Oy/p: bitki başına yağ verimi (cm3), Oy/ha: hektar başına yağ verimi (L) ve NU: bitki başına şemsiye sayısıdır. Hata çubukları ortalamanın standart hatalarını (SEM) temsil etmektedir.. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
Şemsiye sayısı, yağ yüzdesi, bitki başına yağ verimi, hektar başına yağ verimi, bitki boyu, dal sayısı, bitki taze ağırlığı ve bitki kuru ağırlığı dahil olmak üzere tüm özellikler PC₁ üzerinde güçlü ve benzer pozitif yüklemeler göstermiştir (Tablo 2; Şekil 5). Buna göre PC₁, işlemler arasındaki biyokütle ve yağ verimi özelliklerindeki varyasyonu özetleyen genel bir üretkenlik gradyanı olarak yorumlanabilir. PC₂ varyansın yalnızca küçük bir kısmını açıklamış ve işlemler arasındaki ek ayrışmaya sınırlı katkıda bulunmuştur.

Şekil 5Ölçülen özelliklere dayalı olarak 15 uygulama kombinasyonunun dağılımını gösteren PCA skor grafiği. Uygulamalar şu şekilde tanımlanmıştır: 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), T1 (Ct3 × F1), T12 (Ct3 × F2), T13 (Ct3 × F3), T14 (Ct3 × F4) ve T15 (Ct3 × F5). Hata çubukları ortalamanın standart hatalarını (SEM) temsil etmektedir. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
PCA skorları ve kontrolle karşılaştırma
Yüksek pozitif PC₁ skorlarına sahip uygulamalar (özellikle T13, T12, T15 ve T14), çoğu özellik genelinde kontrole (T1) kıyasla en yüksek yüzde artışlarını göstermiştir (Şekil 5 ve Şekil 6). Örneğin T13, kontrole kıyasla şemsiye sayısı bakımından %94,4 ve hektar başına yağ verimi bakımından %98,18 artış göstermiştir.

Şekil 6Ölçülen özelliklerin kontrol uygulamasına (T1: Ct1 × F1) kıyasla göreceli yüzde artışlarıUygulamalar şu şekilde tanımlanmıştır: 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), T1 (Ct3 × F1), T12 (Ct3 × F2), T13 (Ct3 × F3), T14 (Ct3 × F4) ve T15 (Ct3 × F5). Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
Orta düzeyde PC₁ skorlarına sahip tedaviler (örneğin, T8 ve T10), kontrole göre orta dereceli iyileşmeler gösterirken; kontrol ve diğer Ct1 kombinasyonları dahil olmak üzere negatif PC₁ skorlarına sahip tedaviler, nispeten minimal değişiklikler sergilemiştir. Bu örüntüler, PC₁ skorlarının bazal düzeye göre ayarlanmış tedavi yanıtlarıyla yakından ilişkili olduğunu göstermekte ve tek değişkenli ANOVA sonuçlarını desteklemektedir. Bununla birlikte, PCA kesinlikle veri indirgeme ve çok değişkenli örüntülerin görselleştirilmesi amacıyla kullanılmış olup çıkarımsal istatistiksel testlerin yerini tutmamaktadır.
PCA biplotundaki özellik ilişkileri
PCA bi-plotu (Şekil 4); bitki başına yağ verimi, hektar başına yağ verimi ve bitki kuru ağırlığı arasındaki yakın ilişkileri göstermekte olup, bu değişkenler arasındaki güçlü korelasyonları yansıtmaktadır. Bitki boyu, bitki taze ağırlığı ve şemsiye sayısı da benzer vektör yönelimleri sergilemiş ve bu durum, veri seti içindeki bu özellikler arasında koordineli bir varyasyon olduğunu göstermiştir. Dal sayısı, diğer verim bileşenlerinden kısmen ayrışmış olup, mevcut deneysel koşullar altında daha zayıf ilişkiler olduğunu ortaya koymuştur. Bu ilişkiler yalnızca gözlemlenen korelasyonları tanımlamakta olup, mekanistik bağlantılar anlamına gelmemektedir.
60 L ha⁻1 kompost çayı ile yaprak gübresi biostimülanların kombinasyonunun uygulandığı uygulamalar (T12–T15), büyüme ve yağ özellikleri genelinde üstün çok değişkenli performans göstermiştir. Buna karşın, uygulamalar (T7–T11), değerlendirilen çoğu parametre için karşılaştırmalı olarak daha düşük değerler sergilemiştir. Yaprak spreyi uygulanmayan veya Nostoc ile kombine edilen Ct1 (360 L ha⁻1) gibi çoğu vektörün zıt yönünde konumlanan uygulamalar, ölçülen çoğu özellik için daha düşük değerlerle ilişkilendirilmiştir.
Genel olarak PCA, tedavilerin kontrol grubundan birden fazla özellik üzerinden aynı anda nasıl farklılaştığını göstererek ANOVA'yı tamamlayan özlü bir çok değişkenli özet sunmuştur. PCA örüntüleri ile başlangıç değerine göre düzeltilmiş tek değişkenli sonuçlar arasındaki tutarlılık, yorumlamayı ölçülen değişkenlerin ötesine taşımadan tedavi etkilerinin sağlamlığını desteklemektedir.
VERİ ERİŞİLEBİLİRLİĞİ:
Bu makalenin ham verileri ek dosyalar içerisinde mevcuttur.
Ek Tablo 1: Sonuçlar için ham deneysel veriler. Lütfen bu dosyayı indirmek için buraya tıklayın.