Meu gráfico está sendo gerado com a legenda triplicada, ou seja, uma para cada escala. Estou tentando plotar os modelos de estimação para cada linha de tendência na legenda.
Abaixo são apresentados os dados,o código utilizado e a saída com a legenda triplicada:
Os modelos de estimação:
eqs<-list(D1 = "italic(hat(Y))== 2.816 ^\"***\" - 0.094 *x^\"***\" *\",\"~~~italic(R)^{2}== 0.4632",
D2 = "italic(hat(Y))== 3.3635 ^\"***\" - 0.1202 *x^\"***\" *\",\"~~~italic(R)^{2}== 0.622",
D3 = "italic(hat(Y))== 3.7997 ^\"***\" - 0.152 *x^\"***\" *\",\"~~~italic(R)^{2}== 0.7494",
D4 = "italic(hat(Y))== 4.0633 ^\"***\" - 0.167 *x^\"***\" *\",\"~~~italic(R)^{2}== 0.7954")
Exemplo de código do gráfico:
ggplot(data.frame(trat=c(1:4),resp=c(1:4),grupo=c(LETTERS[1:4])),aes(x=trat,y=resp))+
geom_line(color="black",aes(linetype=grupo,group=grupo))+
geom_point(aes(fill=grupo,shape=grupo,group=grupo),size=3,color="black")+
scale_shape_manual("Modelos",labels=parse(text=eqs),values=c(21,22,8,24))+
scale_fill_manual("Modelos",labels=parse(text=eqs),values=viridis(direction = -1, 4))+
scale_linetype_manual("Modelos",labels=parse(text=eqs),values=c(1:4))