多水平统计模型演示文稿.pptVIP

同一医院所诊疗的三名患者的协差阵： 当前第30页\共有74页\编于星期三\2点 医院1：3名患者 ， 医院2：2名患者 应变量向量 Y 总的协方差阵： ? 当前第31页\共有74页\编于星期三\2点 固定与随机参数估计 迭代广义最小二乘算法(Iterative Generalized Least Squares，IGLS) (Goldstein，1986) 限制性迭代广义最小二乘法(Restricted Iterative Generalized Least Squares，RIGLS) (Goldstein，1989) 当前第32页\共有74页\编于星期三\2点 两水平方差成份模型： 当前第33页\共有74页\编于星期三\2点 当前第34页\共有74页\编于星期三\2点 假定已知方差 的值，则可直接构造 残差的协方差矩阵 采用广义最小二乘法(Generalized Least Squares GLS)可获得固定系数的估计： 当前第35页\共有74页\编于星期三\2点 假定方差 未知 在初始阶段，假定 = 0 ，则给出固定系数通常的 OLS 估计 ，得到粗残差： 当前第36页\共有74页\编于星期三\2点 对应于 2 个医院，一个诊疗 3 名患者，另一个诊疗 2 名患者，则 和 均具有 9+4=13 个元素。因为 的期望为 可将这些向量间关系表达为以下线性模型： 当前第37页\共有74页\编于星期三\2点 = +R= + +R 当前第38页\共有74页\编于星期三\2点 为一个残差向量 模型右边包含两个已知的解释变量，其系数即待估计的随机参数 和 通过 GLS 方法获得 和 的估计 回到初始模型则获得固定系数新的估计 在随机与固定参数估计间反复迭代直至收敛 当前第39页\共有74页\编于星期三\2点 分析软件 ML3 (1994) / MLN (1996) / MLwiN (1999) HLM (Hierarchical Linear Model) SAS (Mixed) 当前第40页\共有74页\编于星期三\2点 * 高级卫生统计学 * 例25-1 检索有关吸烟与肺癌关系研究的文献共49篇，各项研究独立且研究假设相同，可以得到的数据包括各项研究的OR值或Ln(OR)值及其标准误、研究水平上的有关解释变量包括样本含量、是否进行性别调整以及设计类型等。 三、应用实例 当前第41页\共有74页\编于星期三\2点 * 高级卫生统计学 * 表25-1 49篇吸烟与肺癌文献的研究结果 三、应用实例 当前第42页\共有74页\编于星期三\2点 * 高级卫生统计学 * 水平 2：文献研究结果（效应尺度） 水平1：个体数据不可得 三、应用实例 为第j个研究的Ln(OR)值 当前第43页\共有74页\编于星期三\2点 * 高级卫生统计学 * 总方差： 三、应用实例 变换为： 定义： 当前第44页\共有74页\编于星期三\2点 * 高级卫生统计学 * 全部Ln(OR)的平均估计值：1.44 三、应用实例 当前第45页\共有74页\编于星期三\2点 * 高级卫生统计学 * 三、应用实例 引入解释变量： ：性别调整 ：设计类型 当前第46页\共有74页\编于星期三\2点 * 高级卫生统计学 * 性别调整系数不显著，但降低了随机效应 效应尺度Ln(OR)与性别调整有关，符号为负 性别调整后Ln(OR)值降低 三、应用实例 当前第47页\共有74页\编于星期三\2点 * 高级卫生统计学 * 三、应用实例 当前第48页\共有74页\编于星期三\2点 * 高级卫生统计学 * 第二节 随机系数模型(Random Coefficient Model) 当前第49页\共有74页\编于星期三\2点 与方差成份模型的区别在于 斜率是随机的： 协变量对反应变量的效应在不同的水平 2 单位间是不同的 当前第50页\共有74页\编于星期三\2点 方差成份模型： 斜率 固定 协变量 对反应变量的效应 在各个医院间不变 随机系数模型：斜率 随机 协变量 对反应变量的效应 在各个医院间是不同的 当前第51页\共有74页\编于星期三\2点 截距 与方差成份模型一致 斜率