二分类广义线性模型
—因变量是二分类变量
【例】用苯氰对小白鼠腹腔注射进行毒性实验,试分析苯氰的剂量和小白鼠死亡率之间的关系,实验结果如下表所示。(本例引用自陈希孺.广义线性模型(一)[J].数理统计与管理, 2002,21(5):54-61)
表1 小白鼠腹腔注射苯氰的毒性实验结果
操作步骤01输入原始数据
图1 原始数据输入结果
02调用命令图2 广义线性模型的命令调用
03选择相关项目首先因本例小白鼠死亡与否是二分类数据,服从二项分布,可用probit回归分析,因而在模型类型对话框中选择“二元概率”。
图3 模型类型选项
敲黑板讲重点
不同类型、不同分布的数据所使用的广义线性模型大不相同。如二分类logistic回归模型,在模型类型上选择“二元Logistic”;等级资料的logistic回归模型,则应选择“有序Logistic”。
在响应对话框中,因变量选择死亡数,试验变量为实验动物。
图4 响应的选项
预测变量的选项框中将“剂量”放入协变量中。
图5 预测变量的选项
模型选项框中,选择剂量作为主效应。
图6 模型效应的选项
图7 广义线性模型估算的选项
图8 广义线性模型统计的选项
结果解读首先输出的结果为模型中的分类变量和连续变量的描述信息和拟合优度情况。
系统根据数据建立的模型,似然比卡方值为97.371,P值小于0.001,模型成立。
模型效应的检验中,变量剂量的P值小于0.001,说明其对因变量有影响。
从参数估计的结果来看,变量剂量的β估计值为0.01,P值小于0.001,因而保留在模型中。Probit模型为:Probit(p)=-3.70+0.01x。
论文表述研究表明,苯氰的剂量和小白鼠的死亡存在关系,其广义线性模型表达为:死亡数=-3.70+0.01剂量(p<0.001)。
微信编辑:汤蝶
原创人员:何少平