Os testes de Shapiro–Wilk indicaram que as distribuições de todas as variáveis medidas não diferiram significativamente da normalidade em ambas as estações (p > 0,05). A análise de variância revelou que a interação entre o nível de chá de compostagem (Ct) e o tratamento com bioestimulante foliar (F) afetou significativamente as características de crescimento vegetativo, o número de umbelas e as características relacionadas ao óleo (percentagem de óleo, produtividade de óleo por planta e produtividade de óleo por hectare) em ambas as estações de cultivo.
Resposta de produtividade em relação ao tratamento controle
Na primeira estação (Figura 2), o rendimento de óleo por hectare (Oy ha⁻1) aumentou progressivamente com o aumento do nível de chá de composto. Todos os tratamentos Ct3 (600 L ha⁻1) produziram maiores rendimentos de óleo do que o tratamento controle (T1: Ct1 × F1, 360 L ha⁻1 sem pulverização foliar). O maior rendimento de óleo (55,58 kg ha⁻1) foi obtido com Ct3 combinado com Spirulina a 2 cm3 L⁻1 (T13), que superou significativamente o controle e todos os tratamentos Ct1, ultrapassando o valor de LSD (2,93 kg ha⁻1). Em relação ao controle, este tratamento apresentou um aumento acentuado no rendimento de óleo por hectare (Figura 3).

Figura 3: Análise de componentes principais em biplot para os quinze tratamentos interativos, com base nas quatro características vegetativas, número de umbelas por planta e nas três características relacionadas ao óleo, utilizando os dois primeiros componentes; cada ponto representa um tratamento. A posição no gráfico indica o grau de semelhança ou diferença entre os tratamentos de acordo com os fatores de entrada (8 características de crescimento e do óleo); quanto mais próximos, mais semelhantes, e quanto mais distantes, mais distintos. PH: altura da planta (cm), NB: número de ramos/planta, HFW: peso fresco total da erva (g), HDW: peso seco total da erva (g), O%: teor de óleo essencial, Oy/p: produtividade de óleo por planta (cm3), Oy/ha: produtividade de óleo por hectare (L) e NU: número de umbelas/planta 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). Clique aqui para visualizar uma versão maior desta figura.
Tendências semelhantes foram observadas para a produtividade de óleo por planta e o número de umbelas, indicando que os tratamentos que resultaram em maior produtividade de óleo por hectare também estiveram associados a uma maior produção reprodutiva e acúmulo de biomassa. Essas relações são apresentadas como associações estatísticas entre características medidas, sem inferir mecanismos causais.
Na segunda estação, Ct3 × F3 (T13) produziu novamente o maior rendimento de óleo por hectare (58,85 kg ha⁻1), excedendo significativamente a testemunha e a maioria dos tratamentos Ct2 (LSD = 3,17 kg ha⁻1). A classificação dos tratamentos em relação à testemunha foi amplamente consistente entre as estações, indicando efeitos reprodutíveis dos tratamentos nas condições ambientais testadas.
Análise de componentes principais (ACP)
A análise de componentes principais foi aplicada como uma ferramenta exploratória para resumir as relações multivariadas entre as oito características medidas ao longo dos 15 tratamentos. Os dois primeiros componentes principais explicaram 98,4% da variância total, sendo que o CP₁ representou 96,5% e o CP₂, 1,9% (Tabela 2; Figura 4). A dominância do CP₁ reflete uma forte covariância entre características relacionadas ao crescimento e à produtividade no conjunto de dados, em vez de processos biológicos subjacentes distintos.
| Parâmetro | 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 |
| Autovalor | 7.7201 | 0.1518 | 0.0721 | 0.0346 |
| % de Variância | 96.5 | 1.9 | 0.9 | 0.43 |
| Acumulado (%) | 96.5 | 98.4 | 99.3 | 99.73 |
| PH: altura da planta (cm), NB: número de ramos/planta, HFW: peso fresco total da erva (g), HDW: peso seco total da erva (g), O%: teor de óleo essencial, Oy/p: rendimento de óleo por planta (cm³), Oy/ha: rendimento de óleo por hectare (L) e NU: número de umbelas/planta |
Tabela 2: Autovalores, autovetores, variância percentual e variância acumulada dos quatro primeiros componentes principais com base em oito características de crescimento e relacionadas ao óleo em 15 combinações de tratamento. As abreviações das características são as seguintes: PH (altura da planta, cm), NB (número de ramos por planta), HFW (peso fresco da erva, g), HDW (peso seco da erva, g), O% (teor de óleo essencial), Oy/p (rendimento de óleo por planta, cm3), Oy/ha (rendimento de óleo por hectare, L) e NU (número de umbelas por planta).

Figura 4: Gráfico de cargas do vetor próprio para 8 características estudadas submetidas aos 15 tratamentos interativos: as barras de PC₁ utilizam o eixo esquerdo; PC₂, PC₃, PC₄ utilizam formas e cores distintas de pontos no eixo direito. PH: altura da planta (cm), NB: número de ramos/planta, HFW: peso fresco total da erva (g), HDW: peso seco total da erva (g), O%: teor de óleo essencial, Oy/p: rendimento de óleo por planta (cm3), Oy/ha: rendimento de óleo por hectare (L) e NU: número de umbelas/planta. Barras de erro representam os erros padrão da média (SEM). Clique aqui para visualizar uma versão maior desta figura.
Todos os caracteres apresentaram cargas fortes e semelhantes positivas no PC₁, incluindo número de umbelas, porcentagem de óleo, produtividade de óleo por planta, produtividade de óleo por hectare, altura da planta, número de ramos, peso fresco da erva e peso seco da erva (Tabela 2; Figura 5). Assim, o PC₁ pode ser interpretado como um gradiente geral de produtividade que resume a variação nos caracteres de biomassa e produtividade de óleo entre os tratamentos. O PC₂ explicou apenas uma pequena proporção da variância e contribuiu com discriminação adicional limitada entre os tratamentos.

Figura 5: Gráfico de escores da ACP ilustrando a distribuição de 15 combinações de tratamentos com base nas características medidas. Os tratamentos são definidos da seguinte forma: 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). As barras de erro representam os erros padrão da média (SEM). Clique aqui para visualizar uma versão maior desta figura.
Escores da ACP e comparação com o controle
Os tratamentos com altos escores positivos no PC₁ (notadamente T13, T12, T15 e T14) também apresentaram os maiores aumentos percentuais em relação ao controle (T1) na maioria das características (Figura 5 e Figura 6). Por exemplo, T13 mostrou aumentos de 94,44% no número de umbelas e de 98,18% no rendimento de óleo por hectare em comparação com o controle.

Figura 6: Aumentos percentuais relativos das características medidas em comparação com o tratamento controle (T1: Ct1 × F1). Os tratamentos são definidos da seguinte forma: 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). Clique aqui para visualizar uma versão maior desta figura.
Tratamentos com escores intermediários de PC₁ (por exemplo, T8 e T10) apresentaram melhorias moderadas em relação ao controle, enquanto tratamentos com escores negativos de PC₁, incluindo o controle e outras combinações Ct1, exibiram mudanças relativas mínimas. Esses padrões indicam que os escores de PC₁ estavam fortemente alinhados com as respostas aos tratamentos ajustadas à linha de base, reforçando os resultados da ANOVA univariada. No entanto, a ACP foi utilizada estritamente para redução de dados e visualização de padrões multivariados e não substitui testes estatísticos inferenciais.
Associações de características no biplot da ACP
O biplot da ACP (Figura 4) ilustra associações estreitas entre rendimento de óleo por planta, rendimento de óleo por hectare e peso seco da erva, refletindo fortes correlações entre essas variáveis. A altura da planta, o peso fresco da erva e o número de umbelas também exibiram orientações vetoriais semelhantes, indicando variação coordenada entre essas características dentro do conjunto de dados. O número de ramos apresentou separação parcial dos demais componentes de produtividade, sugerindo associações mais fracas sob as condições experimentais atuais. Essas relações descrevem apenas correlações observadas e não implicam ligações mecanicistas.
Os tratamentos (T12–T15) que receberam chá de compostagem na dose de 600 L ha⁻1 combinada com bioestimulantes foliares apresentaram desempenho multivariado superior nas características de crescimento e de óleo. Em contraste, os tratamentos (T7–T11) exibiram valores comparativamente mais baixos na maioria dos parâmetros avaliados. Tratamentos posicionados em direção oposta à maioria dos vetores, como Ct1 (360 L ha⁻1) sem aplicação foliar ou combinado com Nostoc, estiveram associados a valores mais baixos na maioria das características mensuradas.
No geral, a ACP forneceu um resumo multivariado conciso que complementa a ANOVA, ilustrando como os tratamentos diferiram do controle em múltiplas características simultaneamente. A consistência entre os padrões da ACP e os resultados univariados ajustados pela linha de base apoia a robustez dos efeitos dos tratamentos, sem ampliar a interpretação para além das variáveis medidas.
DISPONIBILIDADE DE DADOS:
Os dados brutos para este manuscrito estão disponíveis nos arquivos suplementares.
Tabela Suplementar 1: Dados experimentais brutos dos resultados. Clique aqui para baixar este arquivo.