Have: several treatments measured across time, many parameters recorded
dat <- expand.grid(
time = 1:3,
trt = letters[1:2],
param = LETTERS[1:3])
dat$value <- c(
rep(1.0, 6), # A: constant, should drop
rep(c(1.0, 1.5, 2.0), 2), # B: trts identical within days, should drop