[R-sig-ME] Interpreting the coefficients of a main effect
Charles E. (Ted) Wright
cewright at uci.edu
Tue Mar 8 00:27:06 CET 2011
It seems dicey to try to describe an "overall effect of reaction time" based
on this analysis, given that p-value, which may be somewhat suspect, is only
marginally significant, and that there are two important interactions with
reaction time. I would certainly want to explore these interactions
carefully -- are they ordinal or disordinal for example -- before trying to
arrive at any summary.
Ted Wright
On Mon, 7 Mar 2011, Sverre Stausland wrote:
> 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
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
>
More information about the R-sig-mixed-models
mailing list