RE = "all" all random effects are reported RE = "RE" just group effects are reported RE = "OLRE" just observation-level effects are reported
GetOLRE(sigma, model, X, data, RE = c("all", "RE", "OLRE"))