library(AER) data("STAR") STAR<- STAR[complete.cases(STAR), ] STAR<-STAR %>%select(stark,star1,readk,read1,read2) student_ID <- rownames(STAR)