SLIDwomen1 <- read.table("c:/temp/SLID-women.txt",sep = "\t",header=TRUE) SLIDwomen1.glm1 <- glm(as.numeric(working) ~ as.factor(region) +as.numeric(kids0004) +as.numeric(kids0509) +as.numeric(kids1014) +education +familyIncome +I(as.numeric(kids0004)*familyIncome) +I(as.numeric(kids0509)*familyIncome) +I(as.numeric(kids1014)*familyIncome) +I(as.numeric(kids0004)*education) +I(as.numeric(kids0509)*education) +I(as.numeric(kids1014)*education), data=SLIDwomen1,family=binomial()) summary(SLIDwomen1.glm1) SLIDwomen1.glm2 <- glm(as.numeric(working) ~ as.factor(region) +as.numeric(kids0004) +as.numeric(kids0509) +education +familyIncome, data=SLIDwomen1,family=binomial()) summary(SLIDwomen1.glm2)