[R-sig-ME] Interpreting the coefficients of a main effect
Sverre Stausland
johnsen at fas.harvard.edu
Tue Mar 8 00:00:54 CET 2011
Hi mixed-models users,
I am trying to interpret the coefficients I get from lmer models. In
model A, I have these variables:
Dependent variable:
Correct response (0=incorrect, 1=correct)
Independent variables:
Trial (continuous)
ReactionTime (continuous)
Stimulus (binary: 'same' or 'different')
Here is the output:
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.756272 1.344013 2.795 0.00519 **
Trial 0.019962 0.007494 2.664 0.00773 **
ReactionTime -2.354231 1.236785 -1.904 0.05697 .
Stimulussame -5.364448 1.203831 -4.456 8.34e-06 ***
Trial:ReactionTime -0.014469 0.006982 -2.072 0.03822 *
Trial:Stimulussame -0.010312 0.002525 -4.085 4.41e-05 ***
ReactionTime:Stimulussame 4.838763 1.105341 4.378 1.20e-05 ***
I am curious about the overall effect of ReactionTime. If I understand
it correctly, I should add the all the coefficients with ReactionTime
in it to see that, i.e. (-2.25)+(-0.02)+4.84 = 2.57. In other words,
slower responses give more accurate responses.
Question 1:
Is that the correct way to see the overall effect of ReactionTime?
In model B, I do the same as above, but I leave out interactions with
ReactionTime. Now I get this:
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 1.992515 0.598379 3.330 0.000869 ***
Trial 0.006496 0.002061 3.152 0.001622 **
ReactionTime -0.739240 0.452896 -1.632 0.102626
Stimulussame -0.331084 0.353614 -0.936 0.349126
Trial:Stimulussame -0.011652 0.002440 -4.776 1.79e-06 ***
Question 2:
How do I interpret the fact that the overall effect of ReactionTime in
model A with interactions is positive, whereas the overall effect in
model B without interactions is negative? Is the coefficient in model
B simply not trustworthy as a result of leaving out significant
interactions?
Question 3:
Normally I do model comparisons with anova to determine the effect of
a variable. In model A, where I am curious about the significance of
ReactionTime in predicting the dependent variable, should I compare
model A (= the superset model) with a subset model where I remove the
main effect ReactionTime _as well as_ all the interactions it is a
part of?
Thank you
Sverre
More information about the R-sig-mixed-models
mailing list