Fitting a Binomial GLM on Probabilities: A Deep Dive into Logistic Regression for Regression with the Quasibinomial Family Function in R