我对R中的二元逻辑回归感到困惑,其中(显然)是二分结果变量(编码0和1)和二分预测变量(编码0和1)。列联表表明结果是一个非常好的预测因子,但它在我的逻辑回归中并不显着。我发现虚拟问题也有同样的效果,所以我想知道当我使用“完美”预测器时,是否有人可以帮我发现问题?
outcome <- c(0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 1)
predictor <- outcome
model <- glm(outcome ~ predictor, family = binomial)
summary(model)
Call:
glm(formula = outcome ~ predictor, family = binomial)
Deviance Residuals:
Min 1Q Median 3Q
Max
-0.000006547293 -0.000006547293 -0.000006547293 0.000006547293 0.000006547293
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -24.56607 53484.89343 -0.00046 0.99963
predictor 49.13214 79330.94390 0.00062 0.99951
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 15.15820324648649020 on 10 degrees of freedom
Residual deviance: 0.00000000047153748 on 9 degrees of freedom
AIC: 4
Number of Fisher Scoring iterations: 23
我的问题是为什么“预测器”出现了p = .999而不是非常小的东西,因为它应该完全预测结果。提前谢谢。
编辑:如果我将主命令更改为结果~as.factor(预测器),则输出相同