Script_barplot.R

# I used this online tool to extract the data from the scatter plot:

# https://apps.automeris.io/wpd/

# need Hmisc:
install.packages("Hmisc")

# read the data from my mac into R:
ssgac <- read.csv("ssgac.csv", header=FALSE)


PRS<- ssgac$V1
EduYears <- ssgac$V2

college_complete <- ssgac$V2 > 15

quintiles_PRS <- Hmisc::cut2(x = PRS,g=5)

table<- table(college_complete, quintiles_PRS)

barplot(table[2,]/(table[1,]+table[2,]),names.arg = 1:5,ylim =c(0,0.6))