I am calculating a two-level linear multilevel analysis. A look at the random intercept random slope model showed me a significant decrease in my model deviance if I include two dummy variables. Those are from a set of 6 dummyvariables in my multilevel model (which was originally a categorical variable with 7 cagegories). The other dummies arent significant (I use the WALD statistics to judge that).

Can I put just those two dummy variables into the random part of my mixed-model or do I have to usw all six dummy variables even though the other four arent significant? Or do I run into any problem with i.e. the reference group of the six dummies or maybe the interpretation of the two dummies?

I googled, I looked into several books, but I do not find an answer to this.

My entities are countries. My dependent variable is a Likert-Scale. My independent variables on the individual level are mostly dummies, i.e. education is split into six dummy categories (with one category as the reference group). The metric variables are grand_mean centered. On the macro level I have four dummy variables (and a refernce group) and two metric grand_mean centered variables. I calculated the empty model with just my dep.V. Then I put everything into the fixed part and after that I tested with the WALD statistics and the lr-test which variables I should let differ across the countries. My first two education categories and one metric variable were significant and I am unsure if I can let just two categories be random.

Does this help? I certainly can show some code but I do not know which part would be helpful.


Yes, you can. A dummy variable is no different, mathematically, from any other fixed effect you might choose to include. Presumably, if some of the categories all have the same impact on the response, it would make sense to zero them out and push their effect into the intercept term.

That said, using Wald statistics to cull variables is risky. You may get the right set of variables, but this isn't necessarily the case.

When you say "random slope", are you talking about the coefficient of the categorical variable? If so, I would do some model checking. Look at the estimated random effects and see if they are trying to cover small but real differences in your categories.

To clarify that last point: suppose I have 4 categories: A,B,C and D. I decide to omit the dummies associated with C and D. The intercept in the model now corresponds to the case where categories C or D occur. It's like I'm recoding to A, B and Other. But let's suppose that C and D really are real, but just fairly small.

Now fit the random effects model. You will get random intercepts for individuals coded "Other" ... but if you plot these against the true categories (C and D), you might find that the C effects are large and the D effects are small (say), or vice versa. When you add random effects, you are giving a bunch of extra parameters to your model, which it could use to cover up for defects of the model itself.

| cite | improve this answer | |
  • 1
    $\begingroup$ Thanks! I'll also look at the likelihood ratio test to detect the best random slopes and all variables in general. With random slope I mean that the variable can differ across my entities (countries in my case). By zero them out, do you mean I should put them into one category and make them the reference category? Why would that be helpful? But I am very unsure what you mean by the last part. Could you explain that in other words? I am not very confident with this analysis yet. $\endgroup$ – user58673 Oct 15 '14 at 18:59
  • $\begingroup$ @Emma. Welcome to the site. If you find the answer helpful, please upvote (click the grey up arrow) $\endgroup$ – Placidia Oct 15 '14 at 20:08
  • $\begingroup$ By "zero out", I simply meant that the parameters associated with them would be zero. I have edited the last paragraph to make it more explicit. $\endgroup$ – Placidia Oct 15 '14 at 20:09
  • $\begingroup$ Unfortunately I cannot upvote yet, I do not have enough reputations (according to the error message). $\endgroup$ – user58673 Oct 15 '14 at 20:10
  • $\begingroup$ I think I understand. But since I have only one categorie chosen to be the reference group there should be no problem. All dummies are still gonna be in the fixed part, but only two will differ across my entities. $\endgroup$ – user58673 Oct 15 '14 at 20:33

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy