Binary GLMs when response is not a factor

```r
heart <- data.frame(
  ck = 0:11*40 + 20,
  ha = c(2, 13, 30, 30, 21, 19, 18, 13, 19, 15, 7, 8),
  ok = c(88, 26, 8, 5, 0, 1, 1, 1, 1, 0, 0, 0)
)
mod.0 <- glm(cbind(ha,ok) ~ ck, family = binomial(link = "logit"),
             data = heart)
sure::resids(mod.0)
```