Y <- c(1,0,1,0,0,0,1,0,0,0,1,1,0,1,1,0,1,1,1,0)
x1 <- c(1,1,0,1,0,0,0,1,1,0,0,1,0,1,0,1,1,1,0,0)
x2 <- c(1,0,1,1,1,0,1,0,0,0,1,1,0,1,0,0,0,0,1,1)
data <- data.frame(Y,x1,x2)
 reglog<-glm(Y~x1+x2,family=binomial,data=data)
summary(reglog)
## 
## Call:
## glm(formula = Y ~ x1 + x2, family = binomial, data = data)
## 
## Coefficients:
##             Estimate Std. Error z value Pr(>|z|)  
## (Intercept)  -1.0986     0.9522  -1.154   0.2486  
## x1            0.4055     1.0206   0.397   0.6912  
## x2            1.7918     1.0206   1.756   0.0792 .
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 27.726  on 19  degrees of freedom
## Residual deviance: 24.274  on 17  degrees of freedom
## AIC: 30.274
## 
## Number of Fisher Scoring iterations: 4