n <- 3 strata <- rep(1:4, each=n) y <- rnorm(n =12) x <- 1:12 category <- rep(c("A", "B", "C"), times = 4) df <- cbind.data.frame(……