对数线性模型.ppt

第1页，共16页，星期日，2025年，2月5日对数线性模型Logit模型描述的是概率与协变量之间的关系；对数线性模型用来描述期望频数与协变量之间的关系；考虑期望频数m的取值范围在0到无穷之间，故需要进行对数变换为，使它的取值在之间；对数线性模型具有以下形式：不过，与logit不同的是，对数模型中没有解释变量，是用行列因子的效应参数来表示。第2页，共16页，星期日，2025年，2月5日二维列联表的对数线性模型设它的对数线性模型就是对进行分解，分解的方法与方差分析中效应分解的方法完全相同。于是有，其中，是总的平均，和分别是属性A在Ai时和属性B在Bj时的效应，而是属性A和B的交互作用(关联项或关联参数)。以上模型是二维列联表的饱和模型，其期望频数的估计就是实际频数。第3页，共16页，星期日，2025年，2月5日二维列联表的对数线性模型二维列联表的非饱和模型为：其中没有交互效应，经过转换可以得到可见，列联表独立性成立，故称为独立性对数线性模型。其mij的估计为一般认为，在对数线性模型中，当低阶效应为0时，其高阶效应也为0.因此，非饱和模型除以上形式外，还有另外两种情况：第4页，共16页，星期日，2025年，2月5日二维列联表的对数线性模型分别为：其中，属性A与B独立，或仅有属性A的效应，或仅有属性B的效应，且期望频数的估计分别为：独立对数线性模型是否成立的似然比检验和皮尔逊卡方检验统计量为：第5页，共16页，星期日，2025年，2月5日二维列联表的对数线性模型独立性对数线性模型也可以理解为：根据独立性条件，在两边取对数就得到一个相加关系，即依赖于一个与样本量n有关的项，一个与落入第i行的概率有关的项，和一个与落入第j列的概率有关的项，因此独立性具有形式：其中，三个参数分别对应样本量的总效应、行因素的效应和列因素的效应。因此，独立性的原假设等价于该模型的原假设，独立性卡方和似然比检验，就是该模型的拟合优度检验。第6页，共16页，星期日，2025年，2月5日二维列联表的对数线性模型从独立性模型可知，列联表的对数线性模型并不区分响应变量和解释变量，对单元频数的建模，把行和列都看作响应变量。对数线性模型的优点在于：能定量表示属性A在Ai时和属性B在Bj时的效应，以及它们之间的交互效应。【例】对给出的二维列联表(表7.1)，构建对数线性模型。首先根据原列联表可以计算出期望频数估计值，进行独立性检验，,df=(3-1)(3-1)=4,p0.3,故独立性成立。第7页，共16页，星期日，2025年，2月5日二维列联表的对数线性模型在模型检验通过后，来进行参数估计，方法与方差分析完全相同总的效应平均为：属性A在A1时的效应(行效应)为：类似地，可得到属性A在A2，A3时的效应分别为：属性B在B1,B2,B3时的效应(列效应)分别为：第8页，共16页，星期日，2025年，2月5日二维列联表的对数线性模型根据以上独立模型的结果可知，在属性A的效应中，A1的行效应最大为0.48，A3的行效应最小为-0.49；在属性B的效应中，B1和B3的列效应相同为-0.31，B2的列效应最大为0.62；因此，对数线性模型可以定量给出各种效应的值。从以上对数线性模型的应用可以看出，对数模型有假设前提：上例是假定每个因子的效应参数和等于0.在一些软件中通常是假定每个属性的最后一个类别的参数等于0.第9页，共16页，星期日，2025年，2月5日高维列联表的对数线性模型第10页，共16页，星期日，2025年，2月5日高维列联表的对数线性模型与二维列联表的情况类似，高维列联表的对数线性模型也分饱和模型、非饱和模型；对于三维列联表的饱和模型为：其中，除四个单因子项外，三个双因子项，一个三因子项；而它的非饱和对数模型可能有8种情况：对于模型，等价于三维列联表中A,B,C相互独立，记为(A,B,C)；对于模型,等价于A与BC,B与AC,C与AB相互独立；第11页，共16页，星期日，2025年