对数线性模型

关于对数线性模型第1页，课件共43页，创作于2023年2月一、对数线性模型简介1、对数线性模型基本思想

对数线性模型分析是把列联表资料的网格频数的对数表示为各变量及其交互效应的线性模型，然后运用类似方差分析的基本思想，以及逻辑变换来检验各变量及其交互效应的作用大小

第2页，课件共43页，创作于2023年2月区别

方法列联表逻辑回归对数线性模型作用分析定类变量和定类变量之间有无关系分析尺度变量（也可引入类别变量）与二分类别变量之间的因果关系综合运用方差分析和逻辑回归中的建模方法，应用于纯粹定类变量之间，系统评价各变量间关系和交互作用大小的多元统计方法优缺点

不需要确定因变量和自变量。但是，卡方检验对三维和三维以上列联表资料的分析有一定困难，即对混杂变量的控制较难

解决了对混杂变量的控制的问题，而且，它能将因变量与自变量的关系用模型表示出来，清晰易理解。

但是，当模型中自变量较多，特别是名义变量较多，或名义变量的类别较多时，分析自变量之间的交互效应就很繁杂，可能需要建立很多哑变量

可以直接分析各种类型的分类变量，对于名义变量，也不需要事先建立哑变量，可以直接分析变量的主效应和交互效应。对数线性模型不仅可以解决卡方分析中常遇到的高维列联表的“压缩”问题，又可以解决logistic回归分析中多个自变量的交互效应问题

第3页，课件共43页，创作于2023年2月2、列联表的四种类型双向无序列联表；单向有序列联表；双向有序且属性不同的列联表；双向有序且属性相同的列联表第4页，课件共43页，创作于2023年2月3、列联表的优势约束条件少清晰可以快速准确进行判断第5页，课件共43页，创作于2023年2月4、列联表的劣势：对于多关系变量（两个以上）研究：不能被清晰解读失去了对多变量之间的交互联系的分析进行两变量间关联分析时缺乏统计控制不能准确定量描述一个变量对另一个变量的作用幅度第6页，课件共43页，创作于2023年2月5、对数线性模型：多维度列联表解决之道，以及模型自身特点

通过数学方法（方差分析+逻辑变换）来描述多元频数分布。

综合性：同时囊括多个变量于一个模型之中。

控制性：可以在控制其他变量的条件下研究两个分类变量之间的关联。

饱和性：将多元频数分布分解成具体的各项主效应和各项交互效应，以及高阶效应，不会漏项。（饱和模型与不饱和模型）定量性：以发生比的形式来表示自变量的类型不同反映在因变量频数分布上的差异。

可检验性：不仅可以对所有参数估计进行检验，使抽样数据可以推论总体，且能够通过不同模型的统计检验结果，对备选模型进行筛选和评价，进而确定具有最大解释能力且最简单的模型。消除抽样波动所带来的明显的不规则性第7页，课件共43页，创作于2023年2月6、对数线性模型的缺点

对数线性模型更强调的是变量之间的交互效应，它不能直接将因变量用自变量的函数表示出来。对数线性模型抽象复杂，特别是高维模型，不如线性回归模型易理解第8页，课件共43页，创作于2023年2月二、对数线性模型的基本原理1、与方差分析相关的

在多元方差分析中，以二元方差为例：每一个观测值yij=µ+Ai的效果+Bj的效果+(AB)ij交互作用+Ɛij

第9页，课件共43页，创作于2023年2月2、比数比

比数比是对数线性模型的基础，而比数比又是由比数计算而来。那么什么叫做比数呢？比数是一个事件发生的概率与其不发生概率之比,测量了一个事件发生的可能性。这个数值越高说明结果2相对于结果1发生的可能性就越高。

第10页，课件共43页，创作于2023年2月Fij代表某模型fij的期望值，令πij代表与单元格（i,j）有关的期望概率上表可转化为第11页，课件共43页，创作于2023年2月第12页，课件共43页，创作于2023年2月ɯ1=π12/π11ɯ2=π22/π21同理我们可以测量两个两个类别间的比值，称作比数比。Ɵ=ɯ1/ɯ2=π22π21/π12π21=F11F22/F12F21

一个大于1的比数比意味着行变量和列变量的第二个（或者第一个）存在正相关；等于1无关；小于1负相关。

第13页，课件共43页，创作于2023年2月比数比的不变性，不随1）总样本量2）行边缘分布3）列边缘分布的变化而变化。所以，只要关心比数比的估值，那么适用于简单随机样本的最大似然估计就可以被直接应用到分层样本中了。第14页，课件共43页，创作于2023年2月3、与逻辑变换有关的：对数线性模型的出现令R表示行，C表示列，fij表示第i行第j列的观测频次。那么期望频次Fij被设定为一个乘积的函数Fij=ƮƮRiƮCjƮRCijƮ代表概率里面的总概率值1，ƮR和ƮC分别代表R和C的边缘效应，ƮRC代表R与C的二维交互效应，而交互效应实质上测量的就是R与C之间的比数比，当ƮRCij=1的时候就是我们熟悉的独立模型。相乘形式的不好计算，我们将其取对数第15页，课件共43页，创作于2023年2月上两式的数学变换使各种效应项相乘的关系被转换成相加的关系，使各项效应独立化了。

常数效应；

A因素效应；B因素效应；（主效应）A、B两因素的交互效应；第16页，课件共43页，创作于2023年2月主效应和多元交互列表涉及因素数量相等；交互效应的总数则为所有因素各阶组合数之和。对数线性模型有一个限制条件：模型中每一项效应的各类参数之和等于0；如果每项效应中只有一类的参数未知，那么可以由已知参数推算出来。第17页，课件共43页，创作于2023年2月通过上组式子，我们可以计算出线性模型等式右侧的所有参数值。A因素效应是行平均值与总平均值之差B因素效应是列平均值与总平均值之差交互效应计算结果表示在除去所有其他分布效应之后两个因素之间的净关联。

第18页，课件共43页，创作于2023年2月常数项只受样本规模和交互单元数的影响；主效应项反映的是各因素内部类别频数分布的特征，是在总平均频数基础上的“补差”；如果模型中所有交互效应都等于0，我们将会看到虽然每行（列）频数不同，但行（列）频数分布比例却是相同的，都等于原来分类变量的类别分布比例。第19页，课件共43页，创作于2023年2月泊松分布多项分布乘积-多项分布所以我们不能直接应用最小二乘法对模型、总体、参数进行估计，但幸运的是，三个抽样模型下的极大似然估计是等同的。但是可以通过迭代再加权最小二乘法，可是运算起来比较繁琐。

4、分布第20页，课件共43页，创作于2023年2月5、估计参数估计通俗的来讲：根据抽样结果来合理地、科学的猜测一下总体的参数大概是什么？或者是在什么范围？点估计就是用样本计算出来的一个参数来估计未知参数；区间估计就是通过样本计算出来一个范围来对位置参数进行估计。第21页，课件共43页，创作于2023年2月极大似然法与最小二乘法的区别于联系最小二乘法所要解决的问题是：为了选出似的模型输出与系统输出尽可能接近的参数估计，用误差平方和即离差平方和的大小来表示接近程度。使离差平方和最小的参数值即为估计值。简单来说，已知点，自己拟合模型也即分布函数（概率密度函数的积分），进行预测。

极大似然估计所要解决的问题是：选择参数Ɵ，使已知数据在某种意义下最可能出现。某种意义指的是似然函数最大，此处似然函数就是概率密度函数。也就是经常提到的“模型已知，参数未定”。

第22页，课件共43页，创作于2023年2月二者的区别就是，后者需要知道概率密度函数。最小二乘法要的是求出最优的那个参数，而极大似然要求出概率最大（最可能出现的）参数。举个例子，生活中我们一个着眼最合理是哪一个，一个着眼于最可能的是哪一个（极大似然法）当总体服从正态分布时，二者是一样的。对于最小二乘法，当从模型总体随机抽取n组样本观测值后，最合理的参数估计量应该使得模型能最好地拟合样本数据；而对于最大似然法，当从模型总体随机抽取n组样本观测值后，最合理的参数估计量应该是使得从模型中抽取该n组样本观测值的概率最大。第23页，课件共43页，创作于2023年2月密度函数和似然函数（带着参数的密度函数）是相同的，但前者视参数是固定的且数据时变化的，后者视参数变化的且数据时固定的。（1）

写出似然函数；

（2）

对似然函数取对数，并整理；

（3）

求导数

；

（4）

解似然方程第24页，课件共43页，创作于2023年2月三、对数线性模型的假设检验1、假设检验的作用统计推论中包括参数估计与假设检验两部分，上面我们已经介绍了参数估计，那估计的可信度有多少，还要经过假设检验。不经过统计检验，研究者便不能肯定得到的参数估计是不是仅仅源于抽样误差，因而不能肯定在总体中是否存在相同情况。所有结论只能限于这个样本之内，不能肯定再抽一个样本能否得到类似结果。第25页，课件共43页，创作于2023年2月2、统计量似然卡方比，根据相关计算，看原假设是否成立。贝叶斯信息标准，不同模型而言越小的BIC越好。

第26页，课件共43页，创作于2023年2月3、对数线性模型的统计检验四种主要检验：1、对于假设模型的整体检验；2、分层效应的检验；3、单项效应的检验；4、单个参数估计的检验。第27页，课件共43页，创作于2023年2月对数线性模型的统计检验1、对于假设模型的整体检验采用似然比卡方检验（likelihood-ratiochi-squaretest，标为L2）在样本量较大时，L2与皮尔逊卡方统计量的值十分接近。L2优越性： 1、期望频数采用似然估计方法，因而更加稳健； 2、可以被分解成若干部分，即各项效应都有对应的似然卡方值，并且它们的似然卡方值之和等于整个模型的似然卡方比值。第28页，课件共43页，创作于2023年2月对数线性模型的统计检验公式：

其中为估计交互频数。原假设：检验模型的频数估计与观测频数无差异，也可以理解为检验模型和饱和模型无差异。（无关假设）第29页，课件共43页，创作于2023年2月对数线性模型的统计检验饱和对数线性模型可以完美无缺的再现观测频数，因此不需要对饱和模型进行整体性检验。DF等于0，意味着所检验的模型与饱和模型之间的效应项目没有差别。

第30页，课件共43页，创作于2023年2月对数线性模型的统计检验真正有意义的是检验非饱和模型（简略模型，reducedmodel）如果简略模型仍然可以比较准确的拟合观测数据（其拟合程度与饱和模型无显著差异），说明剔除的效应对于拟合意义不大。（科学的简约性原则）研究目的：不是为了再现观测频数，而是通过在模型中加入和减少交互效应项的试验，以寻求真正重要的因素。从饱和模型开始逐步剔除不重要的交互效应项，在保证拟合程度不受较大影响的前提下，直到形成效应项最少的模型。（找到最关键因素）第31页，课件共43页，创作于2023年2月对数线性模型的统计检验举例说明：由图可知，自由度变为1，L2由0增大到10.284，显著性水平α为0.01（P）（拒绝原假设），说明简略模型和饱和模型存在十分显著的差异，即拟合程度受到很大影响。显著=不能剔除该交互因素在因素很多的复杂饱和模型中，通过此方法删减多个不显著效应项来形成简略模型。第32页，课件共43页，创作于2023年2月对数线性模型的统计检验整体检验的不足之处：

整体检验显著只能说明撤销的效应项中起码有一项是有显著作用的，但不能确定是哪一项显著。所以，整体检验在实际对数线性模型分析中，主要服务于整个检验模型的检验情况，而确定各项效应时则是通过单项效应的检验。

且，对于一个多阶多项效应的复杂模型，采用整体检验方式就意味着逐项效应的剔除测试，这样分析过程效率太低。第33页，课件共43页，创作于2023年2月对数线性模型的统计检验2、分层效应检验当研究中涉及的因素较多时，不仅主效应项会增加，交互效应项增加得更快。例如，四个因素的模型，主效应4个，二阶交互效应6项，三阶交互效应4项，四阶交互效应1项。如此，逐项检验筛选重要目标就太繁琐了。且，在一般情况下，高阶交互效应不太容易显著。因此采用按阶次集体检验交互效应项的方法十分间接有效。第34页，课件共43页，创作于2023年2月对数线性模型的统计检验分层效应检验有两种：一、某一阶及更高阶所有交互效应项的集体检验，它的检验是否显著表明这一阶及以上各阶中是否至少有一项是重要的；二、某一阶所有交互效应的集体检验，它的检验是否显著表明这一阶所有交互效应中是否至少有一项是重要的。前者检验比后者综合性更强。第35页，课件共43页，创作于2023年2月对数线性模型的统计检验案例

二阶以上(简略模型)一阶以上一阶二阶第36页，课件共43页，创作于2023年2月对数线性模型的统计检验分层检验提供了模型L2的分解。第一种分层检验中，一阶及以上所有效应都从模型中删除，就会使简略模型的L2增加到13.142，而第二种分层检验告诉我们，这个L2的增量是一阶效应L22.858与二阶效应L210.284之和。第37页，课件共43页，创作于2023年2月对数线性模型的统计检验分层效应检验的不足：

整体检验或分层检验的结果只能说明所有效应中或某一组效应中至少有一项效应具有显著重要影响。但并不能明确知道究竟是哪一项显著。

为了了解到底是哪些具体