mean_set_pl <- mean(iris$Petal.Length[iris$Species == "setosa"])
mean_set_pw <- mean(iris$Petal.Width[iris$Species == "setosa"])
mean_ver_pl <- mean(iris$Petal.Length[iris$Species == "versicolor"])
mean_ver_pw <- mean(iris$Petal.Width[iris$Species == "versicolor"])
mean_vir_pl <- mean(iris$Petal.Length[iris$Species == "virginica"])
mean_vir_pw <- mean(iris$Petal.Width[iris$Species == "virginica"])
scatter_iris <- ggplot(iris,
aes(x=Petal.Length, y=Petal.Width, color=Species)) +
geom_point(alpha=1) +
geom_point(aes(mean_vir_pl, mean_vir_pw),
alpha=.05, size= 7, color="blue") +
geom_point(aes(mean_set_pl, mean_set_pw),
alpha=.05, size= 7, color="red") +
geom_point(aes(mean_ver_pl, mean_ver_pw),
alpha=.05, size= 7, color="green") +
labs(title = "Iris petal length vs width by species")
scatter_iris