是否以及如何在多元回归中使用单尾测试

机器算法验证 r
2022-03-27 09:06:00

我有片面的假设,我想用多重回归来检验。有谁知道如何进行单边假设检验?我听说过关于这个问题的不同立场(“不可能用回归测试单尾”和“对于单边假设,只需将 p 值减半”)。或者也许我的统计程序需要另一个命令(我使用 R)!?我不认为这是一个困难的问题,但找不到一致且令人满意的解决方案。预先感谢您的任何帮助意见。

2个回答

它只需要最少的手动计算来执行有关回归系数的单边假设检验。βi

两种可能的片面假设是:

H0:βi0(1)H1:βi<0

或者

H0:βi0(2)H1:βi>0

提供值用于两侧假设,计算为其中是检验统计量(即回归系数除以其标准误差),是剩余自由度。pR2P(Td|t|)Td

对于第一和第二单边假设,对应的单边值分别是pP(Tdt)P(Tdt)(1)(2)

以下是如何计算 中的单边pR

mod <- lm(Infant.Mortality~., data = swiss)
res <- summary(mod)

# For the two-sided hypotheses

2*pt(-abs(coef(res)[, 3]), mod$df)

(Intercept)   Fertility Agriculture Examination   Education    Catholic 
0.118504431 0.007335715 0.678267621 0.702498865 0.476305225 0.996338704 

# For H1: beta < 0

pt(coef(res)[, 3], mod$df, lower = TRUE)

(Intercept)   Fertility Agriculture Examination   Education    Catholic 
  0.9407478   0.9963321   0.3391338   0.6487506   0.7618474   0.5018306

# For H1: beta > 0

pt(coef(res)[, 3], mod$df, lower = FALSE)

(Intercept)   Fertility Agriculture Examination   Education    Catholic 
0.059252216 0.003667858 0.660866190 0.351249433 0.238152613 0.498169352

mod$df提取剩余自由度并coef(res)[, 3]提取检验统计量。

我不是统计学家,所以对我的回答持保留态度,但这里是。线性回归中报告的 p 值将使用F-test获得,通常用于评估最小二乘拟合问题,如线性回归)。F 检验没有二尾版本,因为 F 值的分布一开始是单尾的(因为 F=t^2)。在此处阅读更详细的答案