Ugh, even shorter, 367

install.packages("modelsummary")
library(modelsummary)
z<-read.csv("entrepreneurs.csv")
summary(z)
aggregate(.~z[,5],z[(3:10)],mean)
s<-function(i){j<-z[z[,5]!=i,]
j$b<-ifelse(j[,5]=="control",0,1)
return(j)}
q<-s("true winners")
g<-function(y){modelsummary(list(lm(y[,10]~b,y),lm(y[,11]~b,y)),
stars=TRUE,gof_map=c("nobs","adj.r.squared"))}
g(q)
e<-s("treated")
g(e)