I want to conduct a mediational analyses with three variables:
- Predictor: it is the result of a memory test with range -1 to 1.
- Mediator: it is the absolute error made by the participant when estimating a magnitude that goes from -1 to 1. Thus, it is a variable that can only take positive values, because it is the absolute deviation from a fixed value. Also, most errors are small so the distribution is skewed (lots of small values).
- Dependent Variable: it is another absolute error made by the participant when estimating a different magnitude. Again, this is a variable with a skewed distribution, all positive values.
My prediction is that the Predictor will affect the Dependent variable through the Mediator. However, I don't think the usual OSL regression with normal functions that most software produces will work here, because of the skewed distributions.
Which approach should I use in this case? I haven't tried but I guess there are a couple things to do: transformations, non-nomral functions (Gamma?)... Still, I prefer receiving input from experts rather than trying things alone.
References are welcome.
EDIT: Thanks to one commenter, I know that I have to look at the residuals. Thus, I have examined the q-q plots for the normal model and the Gamma model (which I think would make sense given my variable).
These are the q-q plots from the Gaussian model – Identity link:
Predictor + Mediator on dependent var.

If I try a Gamma function instead:
Predictor + Mediator on dependent var.

So, I am no expert on q-q plots but these don’t look OK to me. What would you conclude from these plots?
Also: does it make sense to use one distribution for modelling one of the steps of the mediational analysis, and a different one for the others? (I mean, it seems that the residuals look better for Gamma when predicting the mediator, but Gaussian is better when predicting the outcome variable). Prehaps what I am asking makes no sense, sorry.

