Logistic regression with pre-made output
require(rms) mod1b <- lrm(Day0_Active ~ MultyGamer, data = mydata) print(mod1b)