I am using R and I have some issues with replicating the output of the following:
mod1 <- glm(TVAR ~ .,data = df, family = "binomial") y <- predic