14
$\begingroup$

Possible Duplicate:
Under what conditions should one use multilevel/hierarchical analysis?

I have been reading various papers dealing with multilevel analysis, and to be honest, I am still confused re: choosing simple OLS regression or multilevel (random or mixed) models when data is clustered, but the main focus of the research is at individual-level (while group-level effects are treated as a nuisance)? What would be the most "proper" way to proceed in this case?

$\endgroup$
2
  • $\begingroup$ This question is very broad. can you be more specific? $\endgroup$ Oct 1, 2012 at 21:24
  • $\begingroup$ The trouble with OLS is that it assumes the errors are independent, and, if data are clustered, that assumption is violated. $\endgroup$
    – Peter Flom
    Oct 1, 2012 at 21:25

1 Answer 1

27
$\begingroup$

You have two major options:

  1. multilevel analysis that you must have been reading about;
  2. OLS with clustered standard errors (Peter Flom made a comment that OLS assumes that the errors are independent, but that assumption is easy to circumvent with the right choice of the covariance matrix estimator)

Multilevel analysis surely is fancy and hot. That's also the reason it is misused a lot, because everybody seems to want to do something multilevel, no matter whether their data are suitable for it or not. My reaction to about 2/3 of the questions with this tag on this site is that the goals of the study (except for being published in a highly ranked journal, which is often THE goal of many studies) are better addressed by other methods. In multilevel analysis, you have to make strong assumptions: (i) that your random effects are normal (or, if you have random slopes as long as random intercepts, that the joint distribution is multivariate normal), (ii) that your model contains all relevant variables, so that you are safe assuming that errors and regressors are uncorrelated at all levels, (iii) you have enough observations at each level to really utilize the asymptotic theory results concerning the likelihood ratio test statistics and inverse of the information matrix as the estimator of the variances of the parameter estimates. These assumptions are swept under the carpet, most of the time, and rarely if ever checked. The methods that deal with them do exist, but they would require a Ph.D. in statistics to read them. There are also alternative Bayesian solutions which too require a solid stat sequence in Bayesian computing before you even dare to open these papers.

OLS with clustered errors makes fewer assumptions: something like (ii) above, i.e., to be able to convince yourself that the regressors and errors are uncorrelated, and something like (iii), that you have enough clusters so that the variance-covariance estimate is obtained as a sum over sufficiently many independent terms. Note that you don't need to have asymptotics in terms of the number of observations per cluster, unlike multilevel models. An unpleasant side effect concerning OLS with clustered standard errors is that you may run out of degrees of freedom if you have a model with 40 variables and only 30 clusters. (Well if you have 30 clusters, you're screwed anyway.)

An interesting feature of multilevel models is that they can address interactions between levels (e.g., how does the education and experience of a teacher affect the student gains?) It is messier, but possible, to address in OLS as well by explicitly constructing the interactions and using them as explanatory variables in your regression.

With enough data, you can run both analyses and construct a Hausman specification test on the difference between the efficient estimator (multilevel model) and a less efficient and more robust estimator (OLS with clustered standard errors) for the parameters that both models estimate. Most of the time, I would trust the OLS with clustered standard errors more than I would multilevel analysis, frankly.

$\endgroup$
1
  • 7
    $\begingroup$ A little niggle: I don't think you need a phd in statistics to deal with random effects correlated with regressors by adding averages of individual level covariates as predictors for group level intercepts. This fixes assumption ii) across levels, and if you don't like it within level then you ought to think twice about regression models altogether. $\endgroup$ Nov 18, 2013 at 21:16

Not the answer you're looking for? Browse other questions tagged or ask your own question.