SLIDwomen1 <- read.table("c:/temp/SLID-women.txt",sep = "\t",header=TRUE)

SLIDwomen1.glm1 <- glm(working ~ as.factor(region) +I(kids0004)
                  +I(kids0509) +I(kids1014) +education +familyIncome
                  +I(kids0004)*familyIncome +I(kids0509)*familyIncome 
                  +I(kids1014)*familyIncome    +I(kids0004)*education 
                  +I(kids0509)*education 
                  +I(kids1014)*education, data=SLIDwomen1,family=binomial())

summary(SLIDwomen1.glm1)

SLIDwomen1.glm2 <- glm(working ~ as.factor(region) +I(kids0004)
              +I(kids0509) +education +familyIncome, data=SLIDwomen1,family=binomial())

summary(SLIDwomen1.glm2)