edu = read.csv("edu.csv")
Education = edu$Education
Income = edu$Income
Millenial= edu$Millenial

#Plot in 2j
plot(Education, Income, type = "n", main="Scatter Plot", ylim = c(4, 40))
points(jitter(Education[Millenial==F],1), Income[Millenial == F], pch = 1, cex = 1)
points(jitter(Education[Millenial == T],1), Income[Millenial == T], pch = 17, cex = 1)
legend("topright", c("Non-Millenial", "Millenial"), pch = c(1, 17), cex = 1, bty = "n")