我在 R 中运行了这个序数逻辑回归:
mtcars_ordinal <- polr(as.factor(carb) ~ mpg, mtcars)
我得到了这个模型的总结:
summary(mtcars_ordinal)
Re-fitting to get Hessian
Call:
polr(formula = as.factor(carb) ~ mpg, data = mtcars)
Coefficients:
Value Std. Error t value
mpg -0.2335 0.06855 -3.406
Intercepts:
Value Std. Error t value
1|2 -6.4706 1.6443 -3.9352
2|3 -4.4158 1.3634 -3.2388
3|4 -3.8508 1.3087 -2.9425
4|6 -1.2829 1.3254 -0.9679
6|8 -0.5544 1.5018 -0.3692
Residual Deviance: 81.36633
AIC: 93.36633
我可以像这样得到系数的对数几率mpg
:
exp(coef(mtcars_ordinal))
mpg
0.7917679
阈值的对数几率如下:
exp(mtcars_ordinal$zeta)
1|2 2|3 3|4 4|6 6|8
0.001548286 0.012084834 0.021262900 0.277242397 0.574406353
有人可以告诉我我对这个模型的解释是否正确:
随着
mpg
增加一个单位,从类别 1 移动carb
到其他 5 个类别中的任何一个的几率降低 -0.23。如果对数赔率超过 0.0015 的阈值,那么汽车的预测值将是 2 类carb
。如果对数几率超过 0.0121 的阈值,那么汽车的预测值将是 3 类carb
,依此类推。