proc genmod moel y/n = x

a <- c(8, 2, 0, 18, 1, 1, 28, 1, 1

     , 10, 2, 0, 20, 3, 2, 32, 1, 0

     , 12, 3, 0, 22, 2, 1, 34, 1, 1

     , 14, 3, 0, 24, 1, 0, 38, 3, 2

     , 16, 3, 0, 26, 1, 1)

b <- matrix(a, ncol = 3, byrow=T)

head(b)

c <- data.frame(b)

colnames(c) = c('li', 'n', 'y')

    

c <- transform(c,

               b = n -y,

               t = y/n,

               li_sq = li ^ 2)

head(c)

x <- c$li

outcome <- cbind(c$y, c$b)

model.c <- glm(outcome ~ c$li, family = binomial(link = "logit"))

summary(model.c)

model.c1 <- glm(outcome ~ c$li + c$li_sq, family = binomial(link = "logit"))

summary(model.c1)