分层广义线性模型

模型

在分层模型中，观测值可进行聚类，并且观测值的分布不仅由所有聚类的共同结构决定，也由聚类的具体结构决定。于是，模型要引入随机效应分量，并且不同的聚类有不同的随机效应分量。设${\displaystyle y}$ 为响应变量，${\displaystyle u}$ 为随机效应，${\displaystyle g}$  为连结函数。在分层广义线性模型中，需要假设${\displaystyle y|u}$  和 ${\displaystyle u}$ 满足：^[3] ${\displaystyle y\mid u\sim \ f(\theta ,\,\phi )}$  and ${\displaystyle u\sim \ f_{u}(\alpha ).}$ 

线性预测器形式为：

${\displaystyle g(E(y))=g(\mu )=\eta =X\beta +v\,}$ 

其中，${\displaystyle \mu =E(y)}$ ， ${\displaystyle \eta =X\beta +v}$ ， ${\displaystyle v=v(u)}$  为${\displaystyle u}$ 的严格单调函数。在分层广义线性模型中，固定效应为${\displaystyle \beta }$ ，对所有观测值都相同。随机分量${\displaystyle u}$  是不可观测的，不同聚类对应的随机分量取值是随机的。于是，同一聚类的观测值对应的${\displaystyle u}$ 的取值相同，不同聚类的观测值对应的${\displaystyle u}$ 的取值也不同。^[2]

可识别性

为了进行参数推断，有必要保证满足可识别性（英语：Identifiability）^[4] 。在以上模型中，v是不可识别的，因为

${\displaystyle X\beta +v=(X\beta +a)+(v-a)\,}$ 

其中${\displaystyle a}$ 为常数。^[3] 要使模型可识别，需要对参数另加约束。约束常加在随机效应上，比如${\displaystyle E(v)=0}$ 。^[3]

不同分布函数和链接函数的模型

假设不同的分布函数${\displaystyle y\mid u}$  和 ${\displaystyle u}$ ，采用不同的链接函数${\displaystyle g}$  和 '${\displaystyle v}$ ，可以得到不同的模型。另外，广义线性混合模型（英语：generalized linear mixed model）是分层广义线性模型的一个特例。在分层广义线性模型中， 随机效应的分布函数${\displaystyle u}$ 不必要满足正态分布。如果${\displaystyle u}$ 的分布为正态分布，${\displaystyle v}$ 的链接函数为恒等函数，此时的分层广义线性模型即为广义线性混合模型。^[3]

${\displaystyle y\mid u}$  和${\displaystyle u}$ 的分布可取为共轭分布，此时分层广义线性模型有一些特殊的性质，并且易于计算和理解^[3]。比如，如果${\displaystyle y\mid u}$  的分布为平均值一定的泊松分布，${\displaystyle u}$ 的分布为伽玛分布，并取标准对数连接函数，则此时分层广义线性模型为泊松共轭分层广义线性模型。如果${\displaystyle y\mid u}$ 的分布为平均值一定的二项分布，${\displaystyle u}$  的分布为共轭贝塔分布，采用标准罗吉特连接函数，则此时分层广义线性模型为贝塔共轭模型。另外，广义线性混合模型其实就是正态共轭分层广义线性模型。^[3]

常见的模型总结如下：^[5]

分层广义线性模型适用条件是观测值可归为不同的聚类。估计函数有两类：固定效应估计函数和随机效应估计函数，分别相应于 ${\displaystyle \eta =\mathbf {x} {\boldsymbol {\beta }}}$  和 ${\displaystyle \mathbf {v(u)} }$  中的参数。有多种方法进行分层广义线性模型中的参数估计。如果只对固定效应估计函数感兴趣，可以采用总体平均模型。如果要推断个体，就需要估计随机效应。^[2] 拟合分层广义线性模型有多种技术。

分层广义线性模型在实际生活中有诸多应用。

工程

这一模型可用于分析半导体制造中相互关联的过程形成的负载的层级过程^[6][7]。工程师可以应用此模型发现和分析重要的次过程，同时评估这些次过程对最终性能的影响 ^[6]。

商业

市场问题也可以用分层广义线性模型来分析。研究者应用此模型研究了一国范围内的消费者，以解决国际市场研究中的嵌套数据结构问题^[8]。