Penso em dois tipos de gráficos para a sua situação, mas todos ao redor de pacotes diferentes do que você usuou, (base R).
Mas segue como eu faria, pode ser que lhe ajude!
Primeiro gero alguns dados que parecem ter a mesma estrutura que os seus:
library(tidyverse)
n <-100
set.seed(42)
theme_set(theme_minimal())
tb <-
tibble(
coleta = map_chr(1:5 , ~paste('A', . , sep = '')) %>%
rep(20),
tratamento = rbinom(size = 1, n = 100, prob = 0.5),
peso = rnorm(mean = tratamento*0.03, 100),
tratamento_char = tratamento %>% paste('tratamento', .)
)
A primeira maneira, é usando um facet
para cada tratamento, que eu não gosto muito no caso de tratamento em si, pois você não vê o efeito de uma mudança diretamente:
tb %>%
group_by(coleta, tratamento_char) %>%
summarise(avg = mean(peso)) %>%
ggplot(data = . , aes(x = coleta, y = avg)) +
geom_point() +
facet_wrap(~ tratamento_char)
O segundo método é mais minimalista, mas passa bem a mensagem, pois facilita a comparação dos tratamentos em uma coleta:
tb %>%
group_by(coleta, tratamento_char) %>%
summarise(avg = mean(peso)) %>%
ggplot(data = . , aes(x = coleta, y = avg)) +
geom_point(aes(color = tratamento_char,
pch = tratamento_char), size = 2)
O último método, é o meu preferido, pois facilita a comparação entre os grupos e mostra a tendência do todo no experimento.
tb %>%
group_by(coleta, tratamento_char) %>%
summarise(avg = mean(peso)) %>%
ggplot(data = . , aes(x = tratamento_char, y = avg)) +
geom_point(aes(color = coleta)) +
geom_line(aes(group = coleta, color = coleta),
alpha = 0.2,
size = 2.5)