Binary Logistic Regression

The required format is an annotated R script (.R) or a Markdown document (to .pdf or .html).
Use a model with no more than two categorical and two continuous predictors.
To get a 10, make sure your project has the following points correct:
1) the research question is formulated and the choice of predictors is argued;
2) explicit model equation is specified (e.g. log[p(Married)/p(1-Married)] = b0 + b1*Female + b2*Age_centered + b3*Female*Age_centered);
3) scale types and descriptive statistics are provided and analyzed prior to the analysis;
4) model summary and model comparison with interpretations;
5) odds ratios with confidence intervals, predicted probabilities, and average marginal effects;
6) interpretation of all coefficients in the best model;
7) classification accuracy with interpretation;
8) overall goodness of fit measure(s), e.g. pseudo-R-squared measures or the Hosmer-Lemeshow test;
9) diagnostics (linearity and multicollinearity checks); and
10) a one-paragraph conclusion with the results.