# \donttest{
data(sw_c)
dat <- sw_c[sw_c$cluster %in% sort(unique(sw_c$cluster))[1:6], ]
fit <- mrstdlcrt_fit(
data = dat,
formula = y ~ trt + factor(period) + x1 + x2 + (1 | cluster),
cluster_id = "cluster", period = "period", trt = "trt",
method = "lmer", family = "gaussian"
)
plot(fit)
# }
Run the code above in your browser using DataLab