# Linearly regress readk on the interaction between stark and ethnicity
read_int <- lm(readk ~ stark*ethnicity, data=STAR)
read_int
# Bootstrapped standard errors
#read_int_boot <- replicate(500, est.coefs.lm(read_int, resample.cases(STAR)))

#merge(read_int_boot)
#merge(x, y, by="row.names")
#Reduce(function(x, y) merge(x, y, by="row.names"), read_int_boot)
#Reduce(function(x, y) {merge(x, y, by="row.names")}, read_int_boot)
#Reduce(function(...) merge(..., all=T), read_int_boot)

#dat <- as.data.frame(read_int_boot[[1]])
#for(i in 2:length(read_int_boot)){
#  dat <- merge(dat, as.data.frame(read_int_boot[[i]]), by="row.names")
#}

#library(reshape)
#merge_all(read_int_boot)


#rn <- rownames(df1)
#l <- list(df1, df2, df3, df4)
#dat <- l[[1]]
#for(i in 2:length(l)) {
#  dat <- merge(dat, l[[i]],  by= "row.names", all.x= F, all.y= F) [,-1]
#  rownames(dat) <- rn
#}




#read_int_se <- bootstrap.ses(read_int_boot[[500]])

#library(plyr)
#yo <- join_all(as.data.frame(read_int_boot), by="rn", type="inner")