Pseudo Analysis.R

# Load the geepack library

library(geepack)

# Add an ID to each observation, since it's required by the geese function

n <- length(netlixx_pseudo$time)

netlixx_pseudo <- data.frame(netlixx_pseudo, id = 1:n)

# Fit a regression model

summary(fit <- geese(pseudos ~ female + age + coupon,
  data = netlixx_pseudo, id = id, jack = TRUE, family = gaussian,
  corstr = "independence", scale.fix = FALSE))