Ajustei diferentes modelos considerando a variável de resposta (massaseca) em função do (tempo) para cada nível de tratamento (teor) usando o pacote ggplot2 combinada com a função stat_poly_eq.
Porém, como pode ser visto no gráfico a seguir, as legendas das retas estimadas estão sobrepostas. Eu gostaria que eles fossem empilhados no canto esquerdo separadamente. Ao usar a função stat_regline_equation (label.y = 380, label.x = 1000) é possível mover a legenda, no entanto, eles ainda estão sobrepostos.
dados: https://drive.google.com/file/d/1Y-GsNNcYINqtO-hcJfNRgaj545JZXZIS/view?usp=sharing
library(ggplot2)
library(ggpubr)
library(ggpmisc)
my.formula <- y ~ x
ggplot(dadosnew, aes(x = Tempo, y = massaseca, group = interaction(Fator,Trat),
color=interaction(Fator,Trat))) +
stat_summary(geom = "point", fun = mean) +
stat_smooth(method = "lm", se=FALSE, formula=y ~ poly(x, 1, raw=TRUE)) +
stat_poly_eq(formula = my.formula,eq.with.lhs = "As-italic(hat(y))~`=`~",
aes(label = paste(..eq.label.., ..rr.label.., sep = "*plain(\",\")~")),
parse = TRUE, size = 5, label.y = 35)+
labs(title = "",
x = "Time (Minutes)",
y = "Weight (mg)") + theme_bw() +
theme(axis.title = element_text(size = 23,color="black"),
axis.text = element_text(size = 18,color="black"),
text = element_text(size = 20,color="black")) + facet_wrap(~Fator)