Residual errors not normally distributed
The residuals of the linear mixed-effects models in the three studies violated the assumption of normality. Even though linear mixed-effects models tend to be quite robust to deviations from normality (Knief & Forstmeier, 2021; Schielzeth et al., 2020), we sought to verify our results. To this end, we attempted to run robust models using two methods, neither of which worked. The methods are nonetheless described below.
Method A: robustlmm model
The first method drew on the R package ‘robustlmm’ v2.4-4 (Koller, 2016). To calculate the \(p\) values, we followed the procedure of Sleegers et al. (2021), but used the Kenward-Roger method instead of Satterthwaite (see Luke, 2017).
Method B: Inverse Gaussian model with identity link function
In the second approach, we followed a method proposed by Lo and Andrews (2015), based on generalized linear mixed-effects models (GLMM) implementing an identity link function. According to Lo and Andrews (2015), the link function helps avoid directly transforming the dependent variable, which can hinder the interpretability of the results (also see Knief & Forstmeier, 2021).
GLMMs require the use of families of distributions. Lo and Andrews (2015) tested the Gaussian, Gamma and Inverse Gaussian families, with either an identity or an inverse link function. The authors found that the Inverse Gaussian family with an identity link yielded the most normal residuals. The Inverse Gaussian and the Gamma families only accept positive values in the outcome variable (see Table 15.2 in Fox, 2016). Due to this restriction, the dependent variable in the present model is raw RT, unlike the standardised RT that was used in the main analysis.
\(P\) values were to be calculated through parametric bootstrapping, which is the most robust method for GLMMs, as the Kenward-Roger and Satterthwaite methods are not available for these models (Luke, 2017; Singmann et al., 2021).
Neither Method A nor Method B could finally be used, as the code produced errors. These errors are shown in the corresponding scripts inside the ‘model_diagnostics’ folder in each study.
The residuals of the final models are shown in the corresponding sections below.
References
Thesis: https://doi.org/10.17635/lancaster/thesis/1795.
Online book created using the R package bookdown.