library(readxl)
library(ggplot2)
library(dplyr)Gráficos
Visualização do progresso da severidade da doença
Nesta seção, serão gerados gráficos para visualizar a evolução da severidade da ferrugem asiática ao longo do tempo, considerando os diferentes tratamentos. Serão apresentadas duas abordagens: um gráfico único com todas as curvas e gráficos individuais facetados por tratamento.
Gráfico único
Para dar início as análises, vamos carregar os pacotes necessários para essas análises:
Primeiramente, os dados experimentais são importados do arquivo Excel contendo as avaliações de severidade ao longo dos dias para cada unidade experimental.
dados <- read_excel("Trabalho final Emerson.xlsx", sheet = "Planilha1")Os dados são então resumidos para obter a média da severidade por dia para cada tratamento. Esse resumo facilita a visualização da tendência geral do progresso da doença em cada grupo avaliado.
resumo <- dados %>%
group_by(Tratamento, Dia) %>%
summarise(Severidade_media = mean(Severidade, na.rm = TRUE)) %>%
ungroup()Gráfico único com todos os tratamentos
O gráfico a seguir mostra a curva média de severidade da doença ao longo do tempo para todos os tratamentos simultaneamente. Isso permite uma comparação direta entre os tratamentos em termos de velocidade e intensidade da progressão da ferrugem.
ggplot(resumo, aes(x = Dia, y = Severidade_media, color = Tratamento)) +
geom_line(size = 1) +
geom_point(size = 2) +
labs(
title = "Severidade da Doença ao Longo do Tempo",
x = "Dias após Inoculação",
y = "Severidade Média (%)",
color = "Tratamento"
) +
theme_minimal() +
theme(
plot.title = element_text(hjust = 0.5, size = 14, face = "bold"),
axis.title = element_text(size = 12),
legend.position = "right"
)
Criando gráficos individuais
Para facilitar a análise visual de cada tratamento separadamente, os dados também foram representados em gráficos individuais, usando a função facet_wrap. Essa abordagem destaca o comportamento isolado de cada tratamento ao longo do tempo.
ggplot(resumo, aes(x = Dia, y = Severidade_media)) +
geom_line(color = "steelblue", size = 1) +
geom_point(color = "steelblue", size = 2) +
facet_wrap(~ Tratamento) +
labs(
title = "Progresso da Severidade da Doença por Tratamento",
x = "Dias após Inoculação",
y = "Severidade Média (%)"
) +
theme_minimal() +
theme(
strip.text = element_text(face = "bold", size = 12),
plot.title = element_text(hjust = 0.5, size = 14, face = "bold"),
axis.title = element_text(size = 12)
)