Error: $ operator is invalid for atomic vectors for code without $?

Question

Here is my code:

BCG <- c(1,0)
cases <- c(101, 159)
controls <- c(46028, 34594)
BCG_model <- glm(cbind(cases, controls)~BCG, family = binomial)