library(Zelig)

dta <- read.csv("Rwanda.csv")

rwpove <- dta[,1]
treatm <- dta[,2]
nhouse <- dta[,3]
sexhhd <- dta[,4]

llogit <- lm(rwpove ~ treatm + nhouse + sexhhd)
print(summary.lm(llogit))