第7章对数线性模型.pptVIP

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