I do the following graph:
> ddd
UV.NF TRIS volvol
2 145.1923 31 500 μl / 625 μl
3 116.3462 50 500 μl / 625 μl
4 127.1635 60 500 μl / 625 μl
5 125.9615 69 500 μl / 625 μl
6 162.0192 30 1 ml / 625 μl
7 166.8269 50 1 ml / 625 μl
8 176.4423 60 1 ml / 625 μl
9 171.6346 70 1 ml / 625 μl
19 292.3077 31 500 μl / 2500 μl
20 321.1538 50 500 μl / 2500 μl
21 225.0000 60 500 μl / 2500 μl
22 263.4615 69 500 μl / 2500 μl
23 301.9231 30 1 ml / 2500 μl
24 350.0000 50 1 ml / 2500 μl
25 282.6923 60 1 ml / 2500 μl
26 282.6923 70 1 ml / 2500 μl
35 133.6207 31 500 μl / 625 μl
ggplot() +
geom_point(aes(y = log(UV.NF), x = TRIS, colour=ddd[,"volvol"], shape=ddd[,"volvol"]),
data=ddd) +
labs(colour = "volvol", shape="volvol") + xlab("TRIS (mM)") +
guides(colour = guide_legend(title="Vol. lyo. / Vol. reconst."),
shape=guide_legend(title="Vol. lyo. / Vol. reconst.")) +
scale_shape_manual(values = c(19,19,3,3)) + scale_colour_manual(values = c(2,4,2,4))
I want to add the regression line lm(y~x)
for each of the four groups appearing in the legend. I have done many attempts with geom_smooth()
but without success.
See Question&Answers more detail:
os 与恶龙缠斗过久,自身亦成为恶龙;凝视深渊过久,深渊将回以凝视…