Блог

#r

#r

Вопрос:

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

Комментарии:

1. Для меня пример работает нормально (используется R 4.0.3). Можете ли вы сделать это в новом сеансе R

2. Я также использую 4.0.3 и попробовал его в новом сеансе R. Я обнаружил, что добавление кавычек вокруг биномиала заставило его работать.

3. Возможно, у вас есть объект с именем binomial в вашей среде. Попробуйте выполнить rm(binomial) и запустить код снова.