In edxu96/TidyDynamics: Tidy Multivariate (Non)Linear Dynamic Systems

Example 6.4

# Clear variables
rm(list=ls())

data_6.4 <- read.table("seed.data_63")
y <- data_6.4[,1]  # success
n <- data_6.4[,2]  # number of planted seeds

## See Table 6.4
seed <- c(rep(0, 11), rep(1, 10))
extract <- c(rep(0, 5), rep(1, 6), rep(0, 5), rep(1, 5))

head(cbind(y, n, seed, extract))

glm(cbind(y, n - y) ~ factor(seed) * factor(extract), family = binomial)