SPSS学习系列28二元Logistic回归.docx

28. 二元Logistic回归二元或多元线性回归的因变量都是连续型变量，若因变量是分类变量（例如：患病与不患病；不重要、重要、非常重要），就需要用Logistic回归。Logistic回归分析可以从统计意义上估计出在其它自变量固定不变的情况下，每个自变量对因变量取某个值的概率的数值影响大小。Logistic回归模型有“条件”与“非条件”之分，前者适用于配对病例对照资料的分析，后者适用于队列研究或非配对的病例-对照研究成组资料的分析。对于二分类因变量，y=1表示事件发生；y=0表示事件不发生。事件发生的条件概率P{ y=1 | xi } 与 xi 之间是非线性关系，通常是单调的，即随着xi 的增加/减少，P{ y=1 | xi } 也增加/减少。Logistic函数F(x)=，图形如下图所示：该函数值域在(0,1)之间，x趋于-∞时，F(x)趋于0；x趋于+∞时，F(x)趋于1. 正好适合描述概率P{ y=1 | xi }. 例如，某因素x导致患病与否：x在某一水平段内变化时，对患病概率的影响较大；而在x较低或较高时对患病概率影响都不大。记事件发生的条件概率P{ y=1 | xi } = pi，则pi = =记事件不发生的条件概率为1- pi = 则在条件xi下，事件发生概率与事件不发生概率之比为 = 称为事件的发生比，简记为odds. 对odds取自然对数得到上式左边（对数发生比）记为Logit(y), 称为y的Logit变换。可见变换之后的Logit(y)就可以用线性回归，计算出回归系数α和β值。若分类因变量y与多个自变量xi有关，则变换后Logit(y)可由多元线性回归：或 一、简单的二元Logistic回归出现某种结果的概率与不出现的概率之比，称为优势比OR.问题1：研究“低体重出生儿”与“孕妇是否吸烟”之间的关系有数据文件：因变量low：是否“低体重出生儿”（0=正常，1=低体重）；自变量smoke：是否吸烟（0=不吸烟，1=吸烟）【分析】——【回归】——【二元Logistic】，打开“Logistic回归”窗口，将变量“low”选入【因变量】框，变量“smoke”选入【协变量】框；点【确定】，得到因变量编码初始值内部值正常0低出生体重1块 0: 起始块分类表a,b已观测已预测低出生体重儿百分比校正正常低出生体重步骤 0低出生体重儿正常1300100.0低出生体重590.0总计百分比68.8a. 模型中包括常量。b. 切割值为 .500 若模型只含常数项，预测正确率为68.8%（=130/189）；方程中的变量BS.E,WalsdfSig.Exp (B)步骤 0常量-.790.15725.3271.000.454B=-0.79为模型常数项估计值，S.E为B的标准误；Wals为Wald卡方检验，原假设H0：回归系数=0；Exp(B)=0.454（表示患病率与未患病率之比：(1-68.8%)/68.8%）；不在方程中的变量得分dfSig.步骤 0变量smoke4.9241.026总统计量4.9241.026 引入变量后的得分，以及该变量的回归系数是否为0的检验，原假设H0：回归系数=0；（主要针对逐步引进多个变量时的变量筛选）块 1: 方法 = 输入模型系数的综合检验卡方dfSig.步骤 1步骤4.8671.027块4.8671.027模型4.8671.027 似然比卡方值，上一模型（常数项模型）与当前模型似然比值之差，检验两个模型有无差异，原假设H0：无差异。模型汇总步骤-2 对数似然值Cox Snell R 方Nagelkerke R 方1229.805a.025.036a. 因为参数估计的更改范围小于 .001，所以估计在迭代次数 4 处终止。 当前模型的对数似然比值=229.805，两个R方从不同角度反映了当前模型自变量解释因变量的变异占因变量总变异的比例。分类表a已观测已预测低出生体重儿百分比校正正常低出生体重步骤 1低出生体重儿正常1300100.0低出生体重590.0总计百分比68.8a. 切割值为 .500 分类表，概率值0.5判定为患病；与常数项模型结果相同，预测正确率仍为68.8%方程中的变量BS.E,WalsdfSig.Exp (B)步骤 1asmoke.704.3204.8521.0282.022常量-1.087.21525.6271.000.337a. 在步骤 1 中输入的变量: smoke.自变量的偏回归系数及标准误、Wald卡方、自由度、P值、OR值Exp(B). 从而得到Logistic回归模型：Logit(P)=-1.087+0.704*smoke模型结果解释：产妇在妊娠期间是否吸烟* 低出生体重儿 交叉制表低出生体重儿合计正常低出生体重产妇在妊娠期间是否吸烟不吸烟计数8629115产