分层广义线性模型

✍ dations ◷ 2025-09-06 14:57:08 #回归分析,统计模型

在统计学中,分层广义线性模型(hierarchical generalized linear models (HGLM))可视为广义线性模型的推广。在广义线性模型中,误差分量是统计独立的, 然而这一假设并非总是成立的。即在有些情况下,误差项之间有函数关系。分层广义线性模型允许有不同的误差分量,误差分量可以统计相关的,并不必要满足正态分布。当有不同的聚类存在时,同一聚类中的观测值是相关的,并且是正相关的。在这种情况下,广义线性模型是不适用的,忽略这些关联会引起造成一些问题 。

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

线性预测器形式为:

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

为了进行参数推断,有必要保证满足可识别性(英语:Identifiability) 。在以上模型中,v是不可识别的,因为

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

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

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

常见的模型总结如下:

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

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

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


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

相关

  • 风险比风险比率,正式的英文名称是Hazard Ratio。风险比率是两个风险率(Hazard Rate(英语:Hazard Rate))的比值。风险率是单位时间内发生的事件数占被试总体的百分比。瞬时风险率就是当时
  • 板甲板甲(英语:Plate armour)是以金属制成的一整块金属板盔甲、用于保护全身。从古典时期就开始使用,盛行于中古欧洲封建时代。古希腊罗马时期使用青铜材料,14世纪出现了用熟铁浇铸的
  • 世界粮食安全首脑会议世界粮食安全首脑会议于2009年11月16日至18日在意大利罗马举行。举办这次首脑会议的决定是根据粮农组织总干事雅克·迪乌夫博士的提议,由联合国粮食及农业组织理事会在2009年
  • 卡尔大公卡尔·路德维希大公(Archduke Karl Ludwig of Austria,1833年7月30日-1896年5月19日),奥匈帝国皇帝弗朗茨·约瑟夫一世之弟。侄子皇储鲁道夫死后,其长子弗朗茨·斐迪南大公成为皇
  • 塞加拉萨卡拉(阿拉伯语:سقارة‎;英语:Saqqara)是埃及境内一个古代大型墓地,位于开罗以南约30公里,面积约7公里乘1.5公里。该地现时仍屹立着最古老金字塔——阶梯金字塔,该金字塔在前2
  • 国家人权博物馆国家人权博物馆为中华民国设立的国家级人权博物馆,辖新北市新店区白色恐怖景美纪念园区及台东县绿岛乡白色恐怖绿岛纪念园区。其于2018年3月15日成立,同年5月18日揭牌。业务包
  • Psyllipsocidae叶啮虫科(学名:Psyllipsocidae),又名裸啮虫科、跳啮科,是昆虫纲啮虫目之下的一个科,含4属。根据Catalogue of Life、Dyntaxa 及TaiBIF,本科包括下列各物种: .mw-parser-output table
  • 内布拉斯加-林肯大学内布拉斯加大学林肯分校(University of Nebraska–Lincoln,简称UNL、NU,又译内布拉斯加-林肯大学)于1869年根据《土地拨赠法案》创建。是内布拉斯加大学的最主要成员及最早的分
  • 雷姆 (单位)雷姆(全称为人体伦琴当量,英语:roentgen equivalent man,符号为rem)为辐射剂量当量的单位,相当于一伦琴的X光射线或伽码射线。1956年国际放射防护委员会把放射性工作人员的防护标
  • 阿瑟·亚当斯阿瑟·亚当斯(英语:Arthur Adams,1820年-1878年)是英格兰医生和博物学家。亚当斯是1843年至1846年萨马朗号(英语:HMS Samarang (1822))调查亚洲东部群岛的随行外科医师。作为软体动