The code works as I want, but the points are too big. How to decrease them? (Where to insert: size = 0.8?) p1 <- p + geom_point(aes(size = Stage), alpha = 1/3) + xlab ("X") + ylab("Y") + geom_smooth()