As an example I use Pinheiro, J. C. & Bates, D. M. 2000. Mixed-effects models in S and S-PLUS. Springer, New York. page 225. Rats whose body mass has been measured are fed by 3 different diets over time.
Response: Body mass, fixed effects Time*Diet, random effect ~Time|Rat. The main question with this test was if the interaction term is significant (i.e. difference in growth rate between diets). However, my question is could I also look at the p-values of the main effects to say that body mass increased significantly with time for Diet1 (which is the "dummy variable")?
From Pinheiro, J. C. & Bates, D. M. (2000)
Fixed effects: weight ~Time * Diet Value St.error DF t-value p-value Intercept 251.60 13.068 157 19.254 <.0001 Time 0.36 0.088 13 4.084 0.0001 Diet2 200.78 22.657 13 8.862 <.0001 Diet3 252.17 22.662 157 11.127 <.0001 Time*Diet2 0.60 0.155 157 3.871 0.0002 Time*Diet3 0.30 0.156 157 1.893 0.0602 As stated by Pinheiro & Bates, the growth rate of diet 2 (TimeDiet2) differs significantly from diet 1. Although could I state like this for the effect of time on Diet1: f(x) = 251.60 (+/-13.068) + 0.36 x (+/- 0.088), t = 4.084, p = 0.0001? I have seen that people have claimed that it is wrong to interpret p-values for the main effects when the interaction is significant. And is it more proper to split the data and run the test (weight ~Time) for each diet separately, when only looking at the effect of time on body mass?
Best Answer
If there is interaction in the model, the interpretations for main effects would change. For example, in your model, if there is no interaction between Time and Diet, Diet2 means the difference between Diet2 and Diet1 regardless the value of Time; however, if you add the interaction Time*Diet, Diet2 means the difference between Diet2 and Diet1 when Time equals 0, i.e. the difference of intercepts. This depends on model formula, but not significance.
You can say that "body mass increased significant with Time for Diet1" just based on the significance of Time; but for Diet2 and Diet3, to test the slopes, you may need to test the linear combinations of parameters, say for Diet2, to test the significance of Time+Time*Diet2.
- for
Diet1:weight = 251.60 + 0.36*Time; - for
Diet2:weight = (251.60 + 200.78) + (0.36 + 0.60)*Time; - for
Diet3:weight = (251.60 + 252.17) + (0.36 + 0.30)*Time.
You can split the data and run the test separately, but it is formal to integrate the three regressions into one by using the interaction.
By the way, if you need to test the significance of Time*Diet, you may use anova() since it is actually a factor.
Similar Posts:
- Solved – Mean center the dependent variable for within-subject comparison
- Solved – Can interaction be statistically significant while only one of two main effect are significant
- Solved – Insignificant interaction makes variable of interest significant
- Solved – Insignificant interaction makes variable of interest significant
- Solved – repeated measures ANOVA, crossover trial, R