logistic回归1讲义.ppt

有层次结构的数据 组别 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 甲组 0/7 6/9 3/8 1/7 1/8 0/6 2/7 4/5 3/7 0/8 0/7 1/7 7/7 0/7 1/6 乙组 5/7 0/8 2/7 1/7 2/8 0/8 0/8 0/6 0/7 0/9 2/6 不同处理的孕鼠所产仔鼠骨骼畸形发生情况 目的：比较甲、乙两组的仔鼠畸形发生率 （仔鼠中畸形发生数/该窝仔鼠总数） 杨岷，医学和公共卫生研究常用多水平统计模型 数据具有个体嵌套在更高水平单位里的分级结构特点 同一窝别的仔鼠发生畸形的概率趋于相同，而不同窝别的仔鼠发生畸形的概率则趋于不同——试验中各个仔鼠发生畸形的概率不是完全独立的 拟合单水平Logistic回归模型 β0 Β1（处理因素） Χ2值 P值 -1.749 0.773 4.110 0.043 拒绝H0，接受H1，有差别 拟合2水平Logistic回归模型（考虑了窝别效应） β0 Β1（处理因素） Χ2值 P值 -2.539 1.078 1.177 0.278 不拒绝H0 town1 village1 town2 town3 village2 family1 … … family2 family3 … … … … individual1 individual2 individual3 … … Multilevel models, from 1980s 传统方法与logistic回归 四格表资料与logistic的关系 分层四格表资料与logistic的关系 例 四格表资料 探讨妇女使用雌激素与患子宫内膜癌关系的病例-对照研究 雌激素 使用过 未使用过 病例 55 128 对照 19 164 四格表资料的logistic回归 X=1 表示使用过雌激素 X=0 表示未使用过雌激素 例 分层四格表资料与logistic回归 年龄组 X=1, 使用过 X=0, 未使用过 Y=1 Y=0 Y=1 Y=0 1: 25~29 4 62 2 224 2: 30~34 9 33 12 390 3: 35~39 4 26 33 330 4: 40~44 6 9 65 362 5: 45~49 6 5 93 301 按年龄分层的心肌梗死与近期口服避孕药的关系 Mantel-Haenszel法结果 logistic回归结果 传统方法与logistic回归的关系 单因素病例对照研究的OR，与logistic回归等价； 分层病例对照研究的ORMH，与logistic回归结果近似。 logistic回归是传统方法的扩展。 条件Logistic回归 ? 处理匹配资料 ? 使用COX回归分析方法 ? 只能做危险因素分析，不能做预测 和判别用 流程图：统计推断方法 开始 仅一个变量？ 是 单样本问题 是 否 两样本问题 是 ND？ 是 否 分布为二项分布 是 独立样本？ 否 是 1 McNemar 检验 否 ND？ 是 否 分类数据？ 是 否 R × C列联表 否 2 流程图：统计推断方法 是 是 2 × 2列联表？ 是 否 是 否 符合χ2检验条件？ 否 Fisher精确检验 1 无混杂用两样本χ2检验 若存在混杂用Mantel-Haenszel检验 2 × k列联表？ 研究k个二项比例的趋势？ 是 趋势χ2检验 否 R× C列联表？ R× C表χ2检验 R× C表χ2检验 流程图：统计推断方法 2 研究两变量间关系？ 是 两个连续变量？ 是 一个连续变量， 一个分类变量？ 否 是 有序变量？ 否 是 否 两个分类变量 研究关联性 R× C表χ2检验 否 两个以上变量 结果变量为连续变量 或二分类变量？ 连续变量 二分类变量 事件发生的时间 是否重要？ 是 Logistic回归 否 Generate hypotheses Determine data needed Design data collection form collection data manage data Description numbers, distribution, proportions, rate Simple crosstabulations measures of association, tests of significance confidence intervals Stratified Analysis evaluate for confounding and effect modification Multivariate Analysis as needed Interpret for causality 分别讨论偏回归系数的