分层广义线性模型

✍ dations ◷ 2024-12-23 02:01:03 #回归分析,统计模型

在统计学中,分层广义线性模型(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)} } 中的参数。有多种方法进行分层广义线性模型中的参数估计。如果只对固定效应估计函数感兴趣,可以采用总体平均模型。如果要推断个体,就需要估计随机效应。 拟合分层广义线性模型有多种技术。

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

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


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

相关

  • 教授教授,是一种高等教育体系中的职称。在中国汉、唐的大学中即设有此官职;在现代汉语、日语、及韩语的语境中,多作为英语“Professor”一词的同义语使用,指在现代高等教育机构(例如
  • 细胞色素c氧化酶细胞色素c氧化酶(英语:Cytochrome c oxidase)是一种氧化还原酶,通用名为“细胞色素-c氧化酶”,系统名称为“亚铁细胞色素-c:氧气氧化还原酶”(EC 1.9.3.1)。它是一种存在于细菌或线
  • 自发裂变自发裂变(英语:Spontaneous fission)是一种放射性衰变,只发生于原子量高的化学元素。由于元素的核结合能在原子量约为58个原子质量单位(u)时最高,因此更高质量的原子核会自发性分解
  • 基因库基因库或基因池(英语:gene pool)是指同一种群中全部个体所含有的全部基因的集合。因此,每个个体所拥有的基因都是该族群各个成员所共享的。比如黑人、白人、红人、黄种人皆共享
  • 台积电台湾积体电路制造股份有限公司(英语:Taiwan Semiconductor Manufacturing Company, Limited),通称台积电、台积或TSMC,与旗下公司合称时则称做台积电集团,是全球最大的晶圆代工半
  • ɲ硬颚鼻音(palatal nasal)是辅音的一种,在国际音标中以⟨ɲ⟩表示,在X-SAMPA中以⟨J⟩表示。吴语苏州话的声母gn就是此音。粤语钦廉方言相对广州话来说,部分“j”被“ɲ”替代,如“
  • 萨布素萨布素(?-1701年),富察氏,隶属满洲镶黄旗,清朝早期军事人物,曾参与对俄罗斯帝国的战事。萨布素来自军人世家,由四世祖充顺巴本开始,世代为岳克通鄂城主。充顺巴本有一名叫哈木都的后代
  • 海乌姆诺海乌姆诺(Chełmno,德语:Culm oder Kulm)是位于波兰北部维斯瓦河沿岸的城市。人口约有20,000人。自1999年起,海乌姆诺属库亚维-波美拉尼亚省。在1975至1998年期间,海乌姆诺属托伦
  • 美国网络司令部美国网络司令部(英语:United States Cyber Command USCYBERCOM)是美国军方机构,负责开展网络军事行动及保护军方电脑系统,属于美国国防部的一体化作战司令部。2009年6月23日,由美
  • 天鹅山天鹅山(Swan Hill)是位于澳大利亚维多利亚州墨累河南岸的一个城市。据2011年人口普查,天鹅山有人口9894人。天鹅山在1965年正式设市。