In statistics, hierarchical generalized linear models extend generalized linear models by relaxing the assumption that error components are independent. This allows models to be built in situations where more than one error term is necessary and also allows for dependencies between error terms. The error components can be correlated and not necessarily follow a normal distribution. When there are different clusters, that is, groups of observations, the observations in the same cluster are correlated. In fact, they are positively correlated because observations in the same cluster share some common features. In this situation, using generalized linear models and ignoring the correlations may cause problems.
Overview and model
Model In a hierarchical model, observations are grouped into clusters, and the distribution of an observation is determined not only by common structure among all clusters but also by the specific structure of the cluster where this observation belongs. So a random effect component, different for different clusters, is introduced into the model. Let y {\displaystyle y} be the response, u {\displaystyle u} be the random effect, g {\displaystyle g} be the link function, η = X β {\displaystyle \eta =X\beta } , and v = v ( u ) {\displaystyle v=v(u)} is some strictly monotone function of u {\displaystyle u} . In a hierarchical generalized linear model, the assumption on y | u {\displaystyle y|u} and u {\displaystyle u} need to be made: y ∣ u ∼ f ( θ , ϕ ) {\displaystyle y\mid u\sim \ f(\theta ,\,\phi )} and u ∼ f u ( α ) . {\displaystyle u\sim \ f_{u}(\alpha ).}
The linear predictor is in the form:
g ( E ( y ) ) = g ( μ ) = η = X β + v {\displaystyle g(E(y))=g(\mu )=\eta =X\beta +v\,}
where g {\displaystyle g} is the link function, μ = E ( y ) {\displaystyle \mu =E(y)} , η = X β + v {\displaystyle \eta =X\beta +v} , and v = v ( u ) {\displaystyle v=v(u)} is a monotone function of u {\displaystyle u} . In this hierarchical generalized linear model, the fixed effect is described by β {\displaystyle \beta } , which is the same for all observations. The random component u {\displaystyle u} is unobserved and varies among clusters randomly. So v {\displaystyle v} takes the same value for observations in the same cluster and different values for observations in different clusters.
Identifiability Identifiability is a concept in statistics. In order to perform parameter inference, it is necessary to make sure that the identifiability property holds. In the model stated above, the location of v is not identifiable, since
X β + v = ( X β + a ) + ( v − a ) {\displaystyle X\beta +v=(X\beta +a)+(v-a)\,}
for constant a {\displaystyle a} . In order to make the model identifiable, we need to impose constraints on parameters. The constraint is usually imposed on random effects, such as E ( v ) = 0 {\displaystyle E(v)=0} .
… excerpt ends here. Continue reading the full article.
