“…In our simulation scheme, we fix n i = 4 measurements for each subject i and generate data according to the model of (2) where = (0.0, 0.1, 0.2, 0.4, 0.8) , b ijk ∼ N(0, 0.6) and ijk ∼ N(0, 1) . To generate multicollinear data, using the "EnvStats" package in R, each predictor variable is assumed to arise from N (5,1), and the correlation between predictor variables is taken from the set {0.0, 0.2, 0.5, 0.7, 0.9} . Initial values of the variance components are set to be the estimates obtained from fitting a mixed model with no ridge component.…”