R code: logistic regression