Several sources recommend reporting regression coefficients in a table for every mixed-effects model. For continuous predictors that's fine because I only get one coefficient for that predictor. But what if I have categorical predictors with more than two levels? For example
mod <- lmer(angle~temp*recipe + (1|replicate), data=cake)
summary(mod)
Fixed effects:
Estimate Std. Error df t value
(Intercept) 2.379365 6.199942 262.454162 0.384
temp 0.153714 0.029819 250.000012 5.155
recipeB -3.649206 8.464773 250.000006 -0.431
recipeC -1.941270 8.464773 250.000006 -0.229
temp:recipeB 0.010857 0.042170 250.000006 0.257
temp:recipeC 0.002095 0.042170 250.000006 0.050
How should I report this information? Should I report the coefficients for temp
, recipeB
and recipeC
and their corresponding std. errors? Recipe
has 3 levels (A
,B
,C
) and mod
used A
as reference. Should I also change the reference category in recipe
(eg to B
) and re-run the model?
What about the interactions? Does it make sense to include all of the different combinations in the table?
EDIT: Model with two categorical variables, three levels each
cake2 <- cake[which(cake$temp < 200),]
mod <- lmer(angle~temperature*recipe + (1|replicate), data=cake2)
summary(mod)
Fixed effects:
Estimate Std. Error df t value
(Intercept) 29.33333 1.67210 17.27956 17.543
temperature.L 3.44125 1.12498 112.00000 3.059
temperature.Q -0.08165 1.12498 112.00000 -0.073
recipeA 1.15556 0.91854 112.00000 1.258
recipeC 0.20000 0.91854 112.00000 0.218
temperature.L:recipeA -2.26274 1.59096 112.00000 -1.422
temperature.Q:recipeA -1.19753 1.59096 112.00000 -0.753
temperature.L:recipeC -0.75425 1.59096 112.00000 -0.474
temperature.Q:recipeC 0.81650 1.59096 112.00000 0.513
Here I'm not sure what is irrelevant and can be be left out of the table. For example the way the model is set up, the reference is always temperature=175
and recipe=B
. I think I should also report the interaction effects using other references right? Or will readers still be able to calculate the values of the other effects only using the values from the table above?