I test di Shapiro-Wilk hanno indicato che le distribuzioni di tutte le variabili misurate non si discostavano in modo significativo dalla normalità in nessuna delle due stagioni (p > 0,05). L'analisi della varianza ha rivelato che l'interazione tra il livello di tè di compost (Ct) e il trattamento con biostimolante fogliare (F) ha influenzato in modo significativo le caratteristiche della crescita vegetativa, il numero di ombrelle e le caratteristiche relative all'olio (percentuale di olio, resa di olio per pianta e resa di olio per ettaro) in entrambe le stagioni di crescita.
Risposta in termini di resa rispetto al trattamento di controllo
Nella prima stagione (Figura 2), la resa di olio per ettaro (Oy ha⁻1) è aumentata progressivamente all’aumentare del livello di tè di compost. Tutti i trattamenti Ct3 (600 L ha⁻1) hanno prodotto rese di olio superiori rispetto al trattamento di controllo (T1: Ct1 × F1, 360 L ha⁻1 senza applicazione fogliare). La resa di olio più elevata (55,58 kg ha⁻1) è stata ottenuta con Ct3 combinato con Spirulina alla concentrazione di 2 cm3 L⁻1 (T13), superando significativamente il controllo e tutti i trattamenti Ct1, con un valore superiore a quello del LSD (2,93 kg ha⁻1). Rispetto al controllo, questo trattamento ha mostrato un marcato incremento della resa di olio per ettaro (Figura 3).

Figura 3: Analisi delle componenti principali (biplot) per i quindici trattamenti interagenti, basata sui quattro caratteri vegetativi, sul numero di ombrelli per pianta e sui tre caratteri legati all'olio, utilizzando la prima e la seconda componente; ogni punto rappresenta un trattamento. La posizione all'interno del grafico indica il grado di somiglianza o differenza tra i trattamenti in base ai fattori in ingresso (8 caratteristiche di crescita e dell'olio); più vicini sono tra loro, più simili sono i trattamenti, più distanti sono, più diversi risultano. PH: altezza della pianta (cm), NB: numero di rami/pianta, HFW: peso fresco totale dell'erbetta (g), HDW: peso secco totale dell'erbetta (g), O%: contenuto di olio essenziale, Oy/p: resa di olio per pianta (cm3), Oy/ha: resa di olio per ettaro (L) e NU: numero di ombrelli/pianta 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), e T15 (Ct3 × F5). Cliccare qui per visualizzare una versione ingrandita di questa figura.
Tendenze simili sono state osservate per la resa di olio per pianta e per il numero di ombrelle, indicando che i trattamenti che producono una resa maggiore di olio per ettaro sono associati anche a un maggiore output riproduttivo e all'accumulo di biomassa. Queste relazioni sono presentate come associazioni statistiche tra caratteri misurati, senza inferire meccanismi causali.
Nella seconda stagione, Ct3 × F3 (T13) ha nuovamente prodotto il rendimento oleaginoso più elevato per ettaro (58,85 kg ha⁻1), superando significativamente il controllo e la maggior parte dei trattamenti Ct2 (LSD = 3,17 kg ha⁻1). La classifica dei trattamenti rispetto al controllo è stata in gran parte coerente tra le stagioni, indicando effetti dei trattamenti riproducibili nelle condizioni ambientali testate.
Analisi delle componenti principali (PCA)
L'analisi delle componenti principali è stata applicata come strumento esplorativo per riassumere le relazioni multivariate tra gli otto caratteri misurati nei 15 trattamenti. Le prime due componenti principali spiegano il 98,4% della varianza totale, con la PC₁ che contribuisce per il 96,5% e la PC₂ per l'1,9% (Tabella 2; Figura 4). Il predominio della PC₁ riflette una forte covarianza tra i caratteri legati alla crescita e alla produzione nel dataset, piuttosto che distinti processi biologici sottostanti.
| Parametro | 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 |
| Autovalore | 7.7201 | 0.1518 | 0.0721 | 0.0346 |
| % di Varianza | 96,5 | 1,9 | 0,9 | 0,43 |
| Cumulativo (%) | 96,5 | 98,4 | 99,3 | 99,73 |
| PH: altezza della pianta (cm), NB: numero di rami/pianta, HFW: peso fresco totale dell’erba (g), HDW: peso secco totale dell’erba (g), O%: contenuto di olio essenziale, Oy/p: resa di olio per pianta (cm³), Oy/ha: resa di olio per ettaro (L) e NU: numero di ombrelli/pianta |
Tabella 2: Autovalori, autovettori, varianza percentuale e varianza cumulativa dei primi quattro componenti principali basati su otto caratteristiche legate alla crescita e all'olio da 15 combinazioni di trattamento. Le abbreviazioni dei caratteri sono le seguenti: PH (altezza della pianta, cm), NB (numero di rami per pianta), HFW (peso fresco dell'erbetta, g), HDW (peso secco dell'erbetta, g), O% (contenuto di olio essenziale), Oy/p (resa di olio per pianta, cm3), Oy/ha (resa di olio per ettaro, L) e NU (numero di ombrelli per pianta).

Figura 4: Grafico dei carichi del vettore Eigen per 8 caratteri studiati sottoposti a 15 trattamenti interagenti: le barre di PC₁ utilizzano l'asse sinistro; PC₂, PC₃, PC₄ utilizzano forme e colori distinti dei punti sull'asse destro. PH: altezza della pianta (cm), NB: numero di rami/pianta, HFW: peso fresco totale dell'erbetta (g), HDW: peso secco totale dell'erbetta (g), O%: contenuto di olio essenziale, Oy/p: resa di olio per pianta (cm3), Oy/ha: resa di olio per ettaro (L) e NU: numero di ombrelli/pianta. Le barre di errore rappresentano gli errori standard della media (SEM). Cliccare qui per visualizzare una versione ingrandita di questa figura.
Tutti i caratteri hanno mostrato un'elevata e simile correlazione positiva con PC₁, inclusi il numero di ombrelle, la percentuale di olio, la resa di olio per pianta, la resa di olio per ettaro, l'altezza della pianta, il numero di rami, il peso fresco dell'erbetta e il peso secco dell'erbetta (Tabella 2; Figura 5). Di conseguenza, PC₁ può essere interpretato come un gradiente generale di produttività che riassume la variabilità nei caratteri legati alla biomassa e alla resa di olio tra i diversi trattamenti. PC₂ ha spiegato solo una piccola parte della varianza e ha contribuito in misura limitata alla discriminazione aggiuntiva tra i trattamenti.

Figura 5: Grafico dei punteggi PCA che illustra la distribuzione di 15 combinazioni di trattamento in base ai caratteri misurati. I trattamenti sono definiti come segue: 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) e T15 (Ct3 × F5). Le barre di errore rappresentano gli errori standard della media (SEM). Cliccare qui per visualizzare una versione ingrandita di questa figura.
Punteggi PCA e confronto con il controllo
I trattamenti con elevati punteggi positivi di PC₁ (in particolare T13, T12, T15 e T14) hanno mostrato anche gli aumenti percentuali maggiori rispetto al controllo (T1) nella maggior parte dei caratteri (Figura 5 e Figura 6). Ad esempio, T13 ha mostrato un aumento del 94,44% nel numero di ombrelle e del 98,18% nella resa di olio per ettaro rispetto al controllo.

Figura 6: Aumenti percentuali relativi dei caratteri misurati rispetto al trattamento di controllo (T1: Ct1 × F1). I trattamenti sono definiti come segue: 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) e T15 (Ct3 × F5). Cliccare qui per visualizzare una versione ingrandita di questa figura.
I trattamenti con punteggi PC₁ intermedi (ad esempio T8 e T10) hanno mostrato miglioramenti moderati rispetto al controllo, mentre i trattamenti con punteggi PC₁ negativi, inclusi il controllo e altre combinazioni Ct1, hanno mostrato variazioni relative minime. Questi schemi indicano che i punteggi PC₁ erano strettamente correlati alle risposte ai trattamenti aggiustate rispetto al valore iniziale, confermando i risultati dell'ANOVA univariata. Tuttavia, l'ACP è stata utilizzata esclusivamente per la riduzione dei dati e la visualizzazione degli schemi multivariati e non sostituisce i test statistici inferenziali.
Associazioni tra caratteri nel biplot PCA
Il biplot PCA (Figura 4) illustra associazioni strette tra resa di olio per pianta, resa di olio per ettaro e peso secco dell'erbetta, riflettendo forti correlazioni tra queste variabili. Altezza della pianta, peso fresco dell'erbetta e numero di ombrelli hanno mostrato anche orientamenti vettoriali simili, indicando una variazione coordinata tra questi caratteri all'interno del dataset. Il numero di rami ha mostrato una parziale separazione dagli altri componenti della resa, suggerendo associazioni più deboli nelle attuali condizioni sperimentali. Queste relazioni descrivono solo correlazioni osservate e non implicano collegamenti meccanicistici.
I trattamenti (T12–T15) che hanno ricevuto tè di compost a 600 L ha⁻1 in combinazione con biostimolanti fogliari hanno mostrato prestazioni multivariate superiori per quanto riguarda le caratteristiche di crescita e di olio. Al contrario, i trattamenti (T7–T11) hanno mostrato valori relativamente più bassi nella maggior parte dei parametri valutati. I trattamenti posizionati nella direzione opposta alla maggior parte dei vettori, come Ct1 (360 L ha⁻1) senza applicazione fogliare o combinato con Nostoc, erano associati a valori più bassi nella maggior parte delle caratteristiche misurate.
Nell'insieme, l'ACP ha fornito un riassunto multivariato conciso che integra l'ANOVA illustrando in che modo i trattamenti differivano dal controllo attraverso più caratteristiche simultaneamente. La coerenza tra i modelli dell'ACP e i risultati univariati aggiustati rispetto al valore iniziale sostiene la robustezza degli effetti dei trattamenti senza estendere l'interpretazione oltre le variabili misurate.
DISPONIBILITÀ DEI DATI:
I dati grezzi di questo manoscritto sono disponibili nei file supplementari.
Tabella supplementare 1: Dati sperimentali grezzi dei risultati. Cliccare qui per scaricare il file.