Multilevel model
Part of a series on 
Regression analysis 

Models 
Estimation 
Background 

Multilevel models (also known as hierarchical linear models, linear mixedeffect model, mixed models, nested data models, random coefficient, randomeffects models, random parameter models, or splitplot designs) are statistical models of parameters that vary at more than one level.^{[1]} An example could be a model of student performance that contains measures for individual students as well as measures for classrooms within which the students are grouped. These models can be seen as generalizations of linear models (in particular, linear regression), although they can also extend to nonlinear models. These models became much more popular after sufficient computing power and software became available.^{[1]}
Multilevel models are particularly appropriate for research designs where data for participants are organized at more than one level (i.e., nested data).^{[2]} The units of analysis are usually individuals (at a lower level) who are nested within contextual/aggregate units (at a higher level).^{[3]} While the lowest level of data in multilevel models is usually an individual, repeated measurements of individuals may also be examined.^{[2]}^{[4]} As such, multilevel models provide an alternative type of analysis for univariate or multivariate analysis of repeated measures. Individual differences in growth curves may be examined.^{[2]} Furthermore, multilevel models can be used as an alternative to ANCOVA, where scores on the dependent variable are adjusted for covariates (e.g. individual differences) before testing treatment differences.^{[5]} Multilevel models are able to analyze these experiments without the assumptions of homogeneityofregression slopes that is required by ANCOVA.^{[2]}
Multilevel models can be used on data with many levels, although 2level models are the most common and the rest of this article deals only with these. The dependent variable must be examined at the lowest level of analysis.^{[1]}
Level 1 regression equation
When there is a single level 1 independent variable, the level 1 model is:
[math]\displaystyle{ Y_{ij} = \alpha_{i} + \beta_{ij} X_{ij} + e_{ij} }[/math]
 [math]\displaystyle{ Y_{ij} }[/math] refers to the score on the dependent variable for an individual observation at Level 1 (subscript i refers to individual case, subscript j refers to the group).
 [math]\displaystyle{ X_{ij} }[/math] refers to the Level 1 predictor.
 [math]\displaystyle{ \alpha_{i} }[/math] refers to the intercept of the dependent variable for individual case i.
 [math]\displaystyle{ \beta_{ij} }[/math] refers to the slope for individual case i for the relationship in group j (Level 2) between the Level 1 predictor and the dependent variable.
 [math]\displaystyle{ e_{ij} }[/math] refers to the random errors of prediction for the Level 1 equation (it is also sometimes referred to as [math]\displaystyle{ r_{ij} }[/math]).
[math]\displaystyle{ e_{ij} \sim \mathcal{N}(0,\sigma_3^2) }[/math]
At Level 1, both the intercepts and slopes in the groups can be either fixed (meaning that all groups have the same values, although in the real world this would be a rare occurrence), nonrandomly varying (meaning that the intercepts and/or slopes are predictable from an independent variable at Level 2), or randomly varying (meaning that the intercepts and/or slopes are different in the different groups, and that each have their own overall mean and variance).^{[2]}^{[4]}
When there are multiple level 1 independent variables, the model can be expanded by substituting vectors and matrices in the equation.
When the relationship between the response [math]\displaystyle{ Y_{ij} }[/math] and predictor [math]\displaystyle{ X_{ij} }[/math] can not be described by the linear relationship, then one can find some non linear functional relationship between the response and predictor, and extend the model to nonlinear mixedeffects model. For example, when the response [math]\displaystyle{ Y_{ij} }[/math] is the cumulative infection trajectory of the [math]\displaystyle{ i }[/math]th country, and [math]\displaystyle{ X_{ij} }[/math] represents the [math]\displaystyle{ j }[/math]th time points, then the ordered pair [math]\displaystyle{ (X_{ij},Y_{ij}) }[/math] for each country may show a shape similar to logistic function.^{[6]}^{[7]}
Level 2 regression equation
The dependent variables are the intercepts and the slopes for the independent variables at Level 1 in the groups of Level 2.
[math]\displaystyle{ \nu_{i} \sim \mathcal{N}(0,\sigma_1^2) }[/math]
[math]\displaystyle{ \tau_{ij} \sim \mathcal{N}(0,\sigma_2^2) }[/math]
[math]\displaystyle{ \alpha_{i} = \gamma + \nu_{i} }[/math]
[math]\displaystyle{ \beta_{ij} = \delta + \tau_{ij} }[/math]
 [math]\displaystyle{ \gamma }[/math] refers to the overall intercept. This is the grand mean of the scores on the dependent variable across all the groups when all the predictors are equal to 0.
 [math]\displaystyle{ \tau_{ij} }[/math] refers to the overall regression coefficient, or the slope, between the dependent variable and the Level 2 predictor.
 [math]\displaystyle{ \nu_{i} }[/math] refers to the deviation of case i from the overall intercept.
 [math]\displaystyle{ \delta }[/math] refers to the overall regression coefficient, or the slope, between the dependent variable and the Level 1 predictor.
Types of models
Before conducting a multilevel model analysis, a researcher must decide on several aspects, including which predictors are to be included in the analysis, if any. Second, the researcher must decide whether parameter values (i.e., the elements that will be estimated) will be fixed or random.^{[2]}^{[5]}^{[4]} Fixed parameters are composed of a constant over all the groups, whereas a random parameter has a different value for each of the groups.^{[4]} Additionally, the researcher must decide whether to employ a maximum likelihood estimation or a restricted maximum likelihood estimation type.^{[2]}
Random intercepts model
A random intercepts model is a model in which intercepts are allowed to vary, and therefore, the scores on the dependent variable for each individual observation are predicted by the intercept that varies across groups.^{[5]}^{[8]}^{[4]} This model assumes that slopes are fixed (the same across different contexts). In addition, this model provides information about intraclass correlations, which are helpful in determining whether multilevel models are required in the first place.^{[2]}
Random slopes model
A random slopes model is a model in which slopes are allowed to vary according to a correlation matrix, and therefore, the slopes are different across grouping variable such as time or individuals. This model assumes that intercepts are fixed (the same across different contexts).^{[5]}
Random intercepts and slopes model
A model that includes both random intercepts and random slopes is likely the most realistic type of model, although it is also the most complex. In this model, both intercepts and slopes are allowed to vary across groups, meaning that they are different in different contexts.^{[5]}
Developing a multilevel model
In order to conduct a multilevel model analysis, one would start with fixed coefficients (slopes and intercepts). One aspect would be allowed to vary at a time (that is, would be changed), and compared with the previous model in order to assess better model fit.^{[1]} There are three different questions that a researcher would ask in assessing a model. First, is it a good model? Second, is a more complex model better? Third, what contribution do individual predictors make to the model?
In order to assess models, different model fit statistics would be examined.^{[2]} One such statistic is the chisquare likelihoodratio test, which assesses the difference between models. The likelihoodratio test can be employed for model building in general, for examining what happens when effects in a model are allowed to vary, and when testing a dummycoded categorical variable as a single effect.^{[2]} However, the test can only be used when models are nested (meaning that a more complex model includes all of the effects of a simpler model). When testing nonnested models, comparisons between models can be made using the Akaike information criterion (AIC) or the Bayesian information criterion (BIC), among others.^{[1]}^{[2]}^{[5]} See further Model selection.
Assumptions
Multilevel models have the same assumptions as other major general linear models (e.g., ANOVA, regression), but some of the assumptions are modified for the hierarchical nature of the design (i.e., nested data).
 Linearity
The assumption of linearity states that there is a rectilinear (straightline, as opposed to nonlinear or Ushaped) relationship between variables.^{[9]} However, the model can be extended to nonlinear relationships.^{[10]} Particularly, when the mean part of the level 1 regression equation is replaced with a nonlinear parametric function, then such a model framework is widely called the nonlinear mixedeffects model.^{[7]}
 Normality
The assumption of normality states that the error terms at every level of the model are normally distributed.^{[9]}^{[disputed – discuss]} However, most statistical software allows one to specify different distributions for the variance terms, such as a Poisson, binomial, logistic. The multilevel modelling approach can be used for all forms of Generalized Linear models.
 Homoscedasticity
The assumption of homoscedasticity, also known as homogeneity of variance, assumes equality of population variances.^{[9]} However, different variancecorrelation matrix can be specified to account for this, and the heterogeneity of variance can itself be modeled.
 Independence of observations (No Autocorrelation of Model's Residuals)
Independence is an assumption of general linear models, which states that cases are random samples from the population and that scores on the dependent variable are independent of each other.^{[9]} One of the main purposes of multilevel models is to deal with cases where the assumption of independence is violated; multilevel models do, however, assume that 1) the level 1 and level 2 residuals are uncorrelated and 2) The errors (as measured by the residuals) at the highest level are uncorrelated.^{[11]}
 Orthogonality of regressors to random effects
The regressors must not correlate with the random effects, [math]\displaystyle{ u_{0j} }[/math]. This assumption is testable but often ignored, rendering the estimator inconsistent.^{[12]} If this assumption is violated, the randomeffect must be modeled explicitly in the fixed part of the model, either by using dummy variables or including cluster means of all [math]\displaystyle{ X_{ij} }[/math] regressors.^{[12]}^{[13]}^{[14]}^{[15]} This assumption is probably the most important assumption the estimator makes, but one that is misunderstood by most applied researchers using these types of models.^{[12]}
Statistical tests
The type of statistical tests that are employed in multilevel models depend on whether one is examining fixed effects or variance components. When examining fixed effects, the tests are compared with the standard error of the fixed effect, which results in a Ztest.^{[5]} A ttest can also be computed. When computing a ttest, it is important to keep in mind the degrees of freedom, which will depend on the level of the predictor (e.g., level 1 predictor or level 2 predictor).^{[5]} For a level 1 predictor, the degrees of freedom are based on the number of level 1 predictors, the number of groups and the number of individual observations. For a level 2 predictor, the degrees of freedom are based on the number of level 2 predictors and the number of groups.^{[5]}
Statistical power
Statistical power for multilevel models differs depending on whether it is level 1 or level 2 effects that are being examined. Power for level 1 effects is dependent upon the number of individual observations, whereas the power for level 2 effects is dependent upon the number of groups.^{[16]} To conduct research with sufficient power, large sample sizes are required in multilevel models. However, the number of individual observations in groups is not as important as the number of groups in a study. In order to detect crosslevel interactions, given that the group sizes are not too small, recommendations have been made that at least 20 groups are needed,^{[16]} although many fewer can be used if one is only interested in inference on the fixed effects and the random effects are control, or "nuisance", variables.^{[4]} The issue of statistical power in multilevel models is complicated by the fact that power varies as a function of effect size and intraclass correlations, it differs for fixed effects versus random effects, and it changes depending on the number of groups and the number of individual observations per group.^{[16]}
Applications
Level
The concept of level is the keystone of this approach. In an educational research example, the levels for a 2level model might be:
 pupil
 class
However, if one were studying multiple schools and multiple school districts, a 4level model could be:
 pupil
 class
 school
 district
The researcher must establish for each variable the level at which it was measured. In this example "test score" might be measured at pupil level, "teacher experience" at class level, "school funding" at school level, and "urban" at district level.
Example
As a simple example, consider a basic linear regression model that predicts income as a function of age, class, gender and race. It might then be observed that income levels also vary depending on the city and state of residence. A simple way to incorporate this into the regression model would be to add an additional independent categorical variable to account for the location (i.e. a set of additional binary predictors and associated regression coefficients, one per location). This would have the effect of shifting the mean income up or down—but it would still assume, for example, that the effect of race and gender on income is the same everywhere. In reality, this is unlikely to be the case—different local laws, different retirement policies, differences in level of racial prejudice, etc. are likely to cause all of the predictors to have different sorts of effects in different locales.
In other words, a simple linear regression model might, for example, predict that a given randomly sampled person in Seattle would have an average yearly income $10,000 higher than a similar person in Mobile, Alabama. However, it would also predict, for example, that a white person might have an average income $7,000 above a black person, and a 65yearold might have an income $3,000 below a 45yearold, in both cases regardless of location. A multilevel model, however, would allow for different regression coefficients for each predictor in each location. Essentially, it would assume that people in a given location have correlated incomes generated by a single set of regression coefficients, whereas people in another location have incomes generated by a different set of coefficients. Meanwhile, the coefficients themselves are assumed to be correlated and generated from a single set of hyperparameters. Additional levels are possible: For example, people might be grouped by cities, and the citylevel regression coefficients grouped by state, and the statelevel coefficients generated from a single hyperhyperparameter.
Multilevel models are a subclass of hierarchical Bayesian models, which are general models with multiple levels of random variables and arbitrary relationships among the different variables. Multilevel analysis has been extended to include multilevel structural equation modeling, multilevel latent class modeling, and other more general models.
Uses
Multilevel models have been used in education research or geographical research, to estimate separately the variance between pupils within the same school, and the variance between schools. In psychological applications, the multiple levels are items in an instrument, individuals, and families. In sociological applications, multilevel models are used to examine individuals embedded within regions or countries. In organizational psychology research, data from individuals must often be nested within teams or other functional units. They are often used in ecological research as well under the more general term mixed models.^{[4]}
Different covariables may be relevant on different levels. They can be used for longitudinal studies, as with growth studies, to separate changes within one individual and differences between individuals.
Crosslevel interactions may also be of substantive interest; for example, when a slope is allowed to vary randomly, a level2 predictor may be included in the slope formula for the level1 covariate. For example, one may estimate the interaction of race and neighborhood to obtain an estimate of the interaction between an individual's characteristics and the social context.
Applications to longitudinal (repeated measures) data
Alternative ways of analyzing hierarchical data
There are several alternative ways of analyzing hierarchical data, although most of them have some problems. First, traditional statistical techniques can be used. One could disaggregate higherorder variables to the individual level, and thus conduct an analysis on this individual level (for example, assign class variables to the individual level). The problem with this approach is that it would violate the assumption of independence, and thus could bias our results. This is known as atomistic fallacy.^{[17]} Another way to analyze the data using traditional statistical approaches is to aggregate individual level variables to higherorder variables and then to conduct an analysis on this higher level. The problem with this approach is that it discards all withingroup information (because it takes the average of the individual level variables). As much as 80–90% of the variance could be wasted, and the relationship between aggregated variables is inflated, and thus distorted.^{[18]} This is known as ecological fallacy, and statistically, this type of analysis results in decreased power in addition to the loss of information.^{[2]}
Another way to analyze hierarchical data would be through a randomcoefficients model. This model assumes that each group has a different regression model—with its own intercept and slope.^{[5]} Because groups are sampled, the model assumes that the intercepts and slopes are also randomly sampled from a population of group intercepts and slopes. This allows for an analysis in which one can assume that slopes are fixed but intercepts are allowed to vary.^{[5]} However this presents a problem, as individual components are independent but group components are independent between groups, but dependent within groups. This also allows for an analysis in which the slopes are random; however, the correlations of the error terms (disturbances) are dependent on the values of the individuallevel variables.^{[5]} Thus, the problem with using a randomcoefficients model in order to analyze hierarchical data is that it is still not possible to incorporate higher order variables.
Error terms
Multilevel models have two error terms, which are also known as disturbances. The individual components are all independent, but there are also group components, which are independent between groups but correlated within groups. However, variance components can differ, as some groups are more homogeneous than others.^{[18]}
Bayesian nonlinear mixedeffects model
Multilevel modeling is frequently used in diverse applications and it can be formulated by the Bayesian framework. Particularly, Bayesian nonlinear mixedeffects models have recently received significant attention. A basic version of the Bayesian nonlinear mixedeffects models is represented as the following threestage:
Stage 1: IndividualLevel Model
[math]\displaystyle{ {y}_{ij} = f(t_{ij};\theta_{1i},\theta_{2i},\ldots,\theta_{li},\ldots,\theta_{Ki} ) + \epsilon_{ij},\quad \epsilon_{ij} \sim N(0, \sigma^2), \quad i =1,\ldots, N, \, j = 1,\ldots, M_i. }[/math]
Stage 2: Population Model
[math]\displaystyle{ \theta_{li}= \alpha_l + \sum_{b=1}^{P}\beta_{lb}x_{ib} + \eta_{li}, \quad \eta_{li} \sim N(0, \omega_l^2), \quad i =1,\ldots, N, \, l=1,\ldots, K. }[/math]
Stage 3: Prior
[math]\displaystyle{ \sigma^2 \sim \pi(\sigma^2),\quad \alpha_l \sim \pi(\alpha_l), \quad (\beta_{l1},\ldots,\beta_{lb},\ldots,\beta_{lP}) \sim \pi(\beta_{l1},\ldots,\beta_{lb},\ldots,\beta_{lP}), \quad \omega_l^2 \sim \pi(\omega_l^2), \quad l=1,\ldots, K. }[/math]
Here, [math]\displaystyle{ y_{ij} }[/math] denotes the continuous response of the [math]\displaystyle{ i }[/math]th subject at the time point [math]\displaystyle{ t_{ij} }[/math], and [math]\displaystyle{ x_{ib} }[/math] is the [math]\displaystyle{ b }[/math]th covariate of the [math]\displaystyle{ i }[/math]th subject. Parameters involved in the model are written in Greek letters. [math]\displaystyle{ f(t ; \theta_{1},\ldots,\theta_{K}) }[/math] is a known function parameterized by the [math]\displaystyle{ K }[/math]dimensional vector [math]\displaystyle{ (\theta_{1},\ldots,\theta_{K}) }[/math]. Typically, [math]\displaystyle{ f }[/math] is a `nonlinear' function and describes the temporal trajectory of individuals. In the model, [math]\displaystyle{ \epsilon_{ij} }[/math] and [math]\displaystyle{ \eta_{li} }[/math] describe withinindividual variability and betweenindividual variability, respectively. If Stage 3: Prior is not considered, then the model reduces to a frequentist nonlinear mixedeffect model.
A central task in the application of the Bayesian nonlinear mixedeffect models is to evaluate the posterior density:
[math]\displaystyle{ \pi(\{\theta_{li}\}_{i=1,l=1}^{N,K},\sigma^2, \{\alpha_l\}_{l=1}^K, \{\beta_{lb}\}_{l=1,b=1}^{K,P},\{\omega_l\}_{l=1}^K  \{y_{ij}\}_{i=1,j=1}^{N,M_i}) }[/math]
[math]\displaystyle{ \propto \pi(\{y_{ij}\}_{i=1,j=1}^{N,M_i}, \{\theta_{li}\}_{i=1,l=1}^{N,K},\sigma^2, \{\alpha_l\}_{l=1}^K, \{\beta_{lb}\}_{l=1,b=1}^{K,P},\{\omega_l\}_{l=1}^K) }[/math]
[math]\displaystyle{ = \underbrace{\pi(\{y_{ij}\}_{i=1,j=1}^{N,M_i} \{\theta_{li}\}_{i=1,l=1}^{N,K},\sigma^2)}_{Stage 1: IndividualLevel Model} \times \underbrace{\pi(\{\theta_{li}\}_{i=1,l=1}^{N,K}\{\alpha_l\}_{l=1}^K, \{\beta_{lb}\}_{l=1,b=1}^{K,P},\{\omega_l\}_{l=1}^K)}_{Stage 2: Population Model} \times \underbrace{p(\sigma^2, \{\alpha_l\}_{l=1}^K, \{\beta_{lb}\}_{l=1,b=1}^{K,P},\{\omega_l\}_{l=1}^K)}_{Stage 3: Prior} }[/math]
The panel on the right displays Bayesian research cycle using Bayesian nonlinear mixedeffects model.^{[20]} A research cycle using the Bayesian nonlinear mixedeffects model comprises two steps: (a) standard research cycle and (b) Bayesianspecific workflow. Standard research cycle involves literature review, defining a problem and specifying the research question and hypothesis. Bayesianspecific workflow comprises three substeps: (b)–(i) formalizing prior distributions based on background knowledge and prior elicitation; (b)–(ii) determining the likelihood function based on a nonlinear function [math]\displaystyle{ f }[/math]; and (b)–(iii) making a posterior inference. The resulting posterior inference can be used to start a new research cycle.
See also
 Hyperparameter
 Mixeddesign analysis of variance
 Multiscale modeling
 Random effects model
 Nonlinear mixedeffects model
 Restricted randomization
References
 ↑ ^{1.0} ^{1.1} ^{1.2} ^{1.3} ^{1.4} Bryk, Stephen W. Raudenbush, Anthony S. (2002). Hierarchical linear models : applications and data analysis methods (2. ed., [3. Dr.] ed.). Thousand Oaks, CA [u.a.]: Sage Publications. ISBN 9780761919049.
 ↑ ^{2.00} ^{2.01} ^{2.02} ^{2.03} ^{2.04} ^{2.05} ^{2.06} ^{2.07} ^{2.08} ^{2.09} ^{2.10} ^{2.11} Fidell, Barbara G. Tabachnick, Linda S. (2007). Using multivariate statistics (5th ed.). Boston ; Montreal: Pearson/A & B. ISBN 9780205459384.
 ↑ Luke, Douglas A. (2004). Multilevel modeling (3. repr. ed.). Thousand Oaks, CA: Sage. ISBN 9780761928799.
 ↑ ^{4.0} ^{4.1} ^{4.2} ^{4.3} ^{4.4} ^{4.5} ^{4.6} Gomes, Dylan G.E. (20 January 2022). "Should I use fixed effects or random effects when I have fewer than five levels of a grouping factor in a mixedeffects model?". PeerJ 10: e12794. doi:10.7717/peerj.12794. PMID 35116198.
 ↑ ^{5.00} ^{5.01} ^{5.02} ^{5.03} ^{5.04} ^{5.05} ^{5.06} ^{5.07} ^{5.08} ^{5.09} ^{5.10} ^{5.11} Cohen, Jacob (3 October 2003). Applied multiple regression/correlation analysis for the behavioral sciences (3. ed.). Mahwah, NJ [u.a.]: Erlbaum. ISBN 9780805822236.
 ↑ Lee, Se Yoon; Lei, Bowen; Mallick, Bani (2020). "Estimation of COVID19 spread curves integrating global data and borrowing information". PLOS ONE 15 (7): e0236860. doi:10.1371/journal.pone.0236860. PMID 32726361. Bibcode: 2020PLoSO..1536860L.
 ↑ ^{7.0} ^{7.1} Lee, Se Yoon; Mallick, Bani (2021). "Bayesian Hierarchical Modeling: Application Towards Production Results in the Eagle Ford Shale of South Texas". Sankhya B 84: 1–43. doi:10.1007/s13571020002458.
 ↑ editor, G. David Garson (10 April 2012). Hierarchical linear modeling : guide and applications. Thousand Oaks, Calif.: Sage Publications. ISBN 9781412998857.
 ↑ ^{9.0} ^{9.1} ^{9.2} ^{9.3} Salkind, Samuel B. Green, Neil J. (2004). Using SPSS for Windows and Macintosh : analyzing and understanding data (4th ed.). Upper Saddle River, NJ: Pearson Education. ISBN 9780131465978. https://archive.org/details/usingspssforwind00samu.
 ↑ Goldstein, Harvey (1991). "Nonlinear Multilevel Models, with an Application to Discrete Response Data". Biometrika 78 (1): 45–51. doi:10.1093/biomet/78.1.45.
 ↑ ATS Statistical Consulting Group. "Introduction to Multilevel Modeling Using HLM 6". http://www.ats.ucla.edu/stat/hlm/seminars/hlm6/outline_hlm_seminar.pdf.
 ↑ ^{12.0} ^{12.1} ^{12.2} Antonakis, John; Bastardoz, Nicolas; Rönkkö, Mikko (2021). "On Ignoring the Random Effects Assumption in Multilevel Models: Review, Critique, and Recommendations" (in en). Organizational Research Methods 24 (2): 443–483. doi:10.1177/1094428119877457. ISSN 10944281. https://jyx.jyu.fi/bitstream/123456789/66704/2/Antonakisym.pdf.
 ↑ McNeish, Daniel; Kelley, Ken (2019). "Fixed effects models versus mixed effects models for clustered data: Reviewing the approaches, disentangling the differences, and making recommendations." (in en). Psychological Methods 24 (1): 20–35. doi:10.1037/met0000182. ISSN 19391463. PMID 29863377. http://doi.apa.org/getdoi.cfm?doi=10.1037/met0000182.
 ↑ Bliese, Paul D.; Schepker, Donald J.; Essman, Spenser M.; Ployhart, Robert E. (2020). "Bridging Methodological Divides Between Macro and Microresearch: Endogeneity and Methods for Panel Data" (in en). Journal of Management 46 (1): 70–99. doi:10.1177/0149206319868016. ISSN 01492063. http://journals.sagepub.com/doi/10.1177/0149206319868016.
 ↑ Wooldridge, Jeffrey M. (20101001) (in en). Econometric Analysis of Cross Section and Panel Data, second edition. MIT Press. ISBN 9780262296793. https://books.google.com/books?id=hSs3AgAAQBAJ&dq=info:T5fz2cmyyF8J:scholar.google.com&pg=PP1.
 ↑ ^{16.0} ^{16.1} ^{16.2} Leeuw, Ita Kreft, Jan de (1998). Introducing multilevel modeling. (Repr. ed.). London: Sage Publications Ltd. ISBN 9780761951414.
 ↑ Hox, Joop (2002). Multilevel analysis : techniques and applications (Reprint. ed.). Mahwah, NJ [u.a.]: Erlbaum. ISBN 9780805832198.
 ↑ ^{18.0} ^{18.1} Bryk, Anthony S.; Raudenbush, Stephen W. (1 January 1988). "Heterogeneity of variance in experimental studies: A challenge to conventional interpretations.". Psychological Bulletin 104 (3): 396–404. doi:10.1037/00332909.104.3.396.
 ↑ Lee, Se Yoon (2022). "Bayesian Nonlinear Models for Repeated Measurement Data: An Overview, Implementation, and Applications". Mathematics 10 (6): 898. doi:10.3390/math10060898.
 ↑ Lee, Se Yoon (2022). "Bayesian Nonlinear Models for Repeated Measurement Data: An Overview, Implementation, and Applications". Mathematics 10 (6): 898. doi:10.3390/math10060898.
Further reading
 Gelman, A.; Hill, J. (2007). Data Analysis Using Regression and Multilevel/Hierarchical Models. New York: Cambridge University Press. pp. 235–299. ISBN 9780521686891. https://books.google.com/books?id=lV3DIdV0F9AC&pg=PA235.
 Goldstein, H. (2011). Multilevel Statistical Models (4th ed.). London: Wiley. ISBN 9780470748657. https://books.google.com/books?id=mdwt7ibSGUYC.
 Hedeker, D.; Gibbons, R. D. (2012). Longitudinal Data Analysis (2nd ed.). New York: Wiley. ISBN 9780470889183.
 Hox, J. J. (2010). Multilevel Analysis: Techniques and Applications (2nd ed.). New York: Routledge. ISBN 9781848728455. https://books.google.com/books?id=jLPHBQAAQBAJ.
 Raudenbush, S. W.; Bryk, A. S. (2002). Hierarchical Linear Models: Applications and Data Analysis Methods (2nd ed.). Thousand Oaks, CA: Sage. This concentrates on education.
 Snijders, T. A. B.; Bosker, R. J. (2011). Multilevel Analysis: an Introduction to Basic and Advanced Multilevel Modeling (2nd ed.). London: Sage. ISBN 9781446254332. https://books.google.com/books?id=N1BQvcomDdQC.
 Swamy, P. A. V. B.; Tavlas, George S. (2001). "Random Coefficient Models". in Baltagi, Badi H.. A Companion to Theoretical Econometrics. Oxford: Blackwell. pp. 410–429. ISBN 9780631212546.
 Verbeke, G.; Molenberghs, G. (2013). Linear Mixed Models for Longitudinal Data. Springer. Includes SAS code
 Gomes, Dylan G.E. (20 January 2022). "Should I use fixed effects or random effects when I have fewer than five levels of a grouping factor in a mixedeffects model?". PeerJ 10: e12794. doi:10.7717/peerj.12794. PMID 35116198.
External links
Original source: https://en.wikipedia.org/wiki/Multilevel model.
Read more 