机器算法验证 - 如何使用逻辑回归进行事后测试？ - 吾爱随笔录

如何使用逻辑回归进行事后测试？

机器算法验证 spss 物流事后对比

2022-03-02 15:40:30

我查看了参与者在收到给定消息后是否回收。分为三组：对照组、A组和B组。

由于我的依赖度量是二分法的（回收与使用垃圾箱），我在 SPSS 中进行了逻辑回归。但是，SPSS 只允许我将各个组与对照组进行比较。我发现 A 组的回收率高于对照组，但 B 组没有。

我现在很好奇的是 A 组的回收率是否明显高于 B 组。

是否可以进行这样的事后测试？如果是这样，怎么做？

3个回答

不幸的是，我不知道 SPSS。也就是说，如果您想在 null 为 $H_0: \beta_{groupA} - \beta_{groupB} = 0$ 您可以向 SPSS 询问参数估计的方差/协方差矩阵并手动构建 Wald 检验。

在下面 $H_0$ 你的测试统计 $\chi^2_{obs}$ 分布为 $\chi^2$ 具有 1 个自由度的 rv

χ_{o b s}^{2} = \frac{({\hat{β}}_{g r o u p A} - {\hat{β}}_{g r o u p B})^{2}}{v a r [{\hat{β}}_{g r o u p A}] + v a r [{\hat{β}}_{g r o u p B}] - 2 * c o v [{\hat{β}}_{g r o u p A}, {\hat{β}}_{g r o u p B}]}

$\chi^2_{obs} = \frac{(\hat{\beta}_{groupA} - \hat{\beta}_{groupB})^2}{{\rm var}[\hat{\beta}_{groupA}]+{\rm var}[\hat{\beta}_{groupB}]-2*{\rm cov}[\hat{\beta}_{groupA},\hat{\beta}_{groupB}]}$

现在您可以计算您的 p 值。但我确信 SPSS 有一个命令可以对估计的参数执行特定的测试。

SPSS 仅让我将各个组与对照组进行比较。

实际上 SPSS Logistic Regression 有大约 6 种内置的对比类型。其中一个（指标）将每个组与一个对照组进行比较，您可以使用组的编号来指定。即，在编号为 1 到 4 并标记为 North、South、East 和 West 的组中，“indicator(3)”将 East 设置为对照组。另一种类型（偏差）显示每个组的 logit 如何偏离（未加权）平均组的 logit。进入通用帮助文件或命令语法参考（也可以在帮助中找到）以查找每种类型的定义很有用。就个人而言，我发现偏差和指标似乎是我所需要的。也许这让我成为了一名小联盟球员:-)

在 R 中，您可以为逻辑回归进行一般多重性调整对比。例如，请参见使用 R包的rms包的功能。contrast.rmsmultcomp

其它你可能感兴趣的问题

上一篇事后测试的奇怪结果下一篇比率的标准误