title: Multilevel bootstrap analysis with assumptions violated.
abstract: BACKGROUND
  Likelihood-based methods can work poorly when the residuals are not normally distributed and the variances across clusters are heterogeneous.
  
  
  METHODS
  The performance of two estimation methods, the non-parametric residual bootstrap (RB) and the restricted maximum likelihood (REML) for fitting multilevel models are compared through simulation studies in terms of bias, coverage, and precision.
  
  
  RESULTS
  We find that (a) both methods produce unbiased estimates of the fixed parameters, but biased estimates of the random parameters, although the REML was more prone to give biased estimates for the variance components; (b) the RB method yields substantial reductions in the difference between nominal and actual confidence interval coverage, compared with the REML method; and (c) for the square root of the mean squared error (RMSE) of the fixed effects, the RB method performed slightly better than the REML method. For the variance components, however, the RB method did not offer a systematic improvement over the REML method in terms of RMSE.
  
  
  CONCLUSIONS
  It can be stated that the RB method is, in general, superior to the REML method with violated assumptions.
authors.affiliation: Universidad de Oviedo.
authors.forename: Guillermo
authors.initials: G
authors.lastname: Vallejo Seco
authors.forename: Manuel
authors.initials: M
authors.lastname: Ato García
authors.forename: María Paula
authors.initials: MP
authors.lastname: Fernández García
authors.forename: Pablo Esteban
authors.initials: PE
authors.lastname: Livacic Rojas
mesh_terms: D000704:Analysis of Variance; D015982:Bias; D016013:Likelihood Functions; D018709:Statistics, Nonparametric