我正在使用 step() 函数在 R 中进行反向消除。现在,我正在尝试查看每个自变量如何与它们的 AIC、F 和 P 值一起排名。
step(Mod1,direction="backward",test="F")
Df Deviance AIC F value Pr(>F)
<none> 6127.4 6215.4
- as.factor(var2) 3 6133.6 6215.6 2.6103 0.0497127 *
- as.factor(var28) 2 6131.7 6215.7 2.7292 0.0653326 .
- as.factor(var32) 2 6131.8 6215.8 2.7794 0.0621388 .
- as.factor(var30) 1 6130.3 6216.3 3.6075 0.0575550 .
- as.factor(var20) 1 6131.9 6217.9 5.7262 0.0167368 *
- as.factor(var9) 1 6133.5 6219.5 7.6627 0.0056507 **
- as.factor(var15) 1 6133.7 6219.7 7.8952 0.0049691 **
- as.factor(var10) 1 6133.8 6219.8 8.1314 0.0043621 **
- as.factor(var14) 1 6134.7 6220.7 9.2528 0.0023592 **
- as.factor(var33) 2 6137.1 6221.1 6.0993 0.0022552 **
- as.factor(var16) 1 6135.9 6221.9 10.6794 0.0010881 **
- as.factor(var19) 4 6142.5 6222.5 4.7684 0.0007674 ***
- as.factor(var23) 2 6138.9 6222.9 7.2488 0.0007158 ***
- as.factor(var24) 2 6139.0 6223.0 7.3060 0.0006761 ***
- as.factor(var13) 1 6139.3 6225.3 14.9746 0.0001099 ***
- as.factor(var11) 1 6141.0 6227.0 17.1558 3.480e-05 ***
- as.factor(var6) 2 6149.3 6233.3 13.8110 1.030e-06 ***
- as.factor(var22) 2 6150.6 6234.6 14.6341 4.534e-07 ***
- as.factor(var8) 4 6155.4 6235.4 8.8624 3.893e-07 ***
- as.factor(var3) 4 6172.7 6252.7 14.3214 1.189e-11 ***
- as.factor(var1) 1 6230.8 6316.8 130.7555 < 2.2e-16 ***
- as.factor(var5) 4 6245.6 6325.6 37.3782 < 2.2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
接下来,我想
i) 根据 p 值从最重要(顶部)到不太重要(底部)对变量进行排序
ii) 获取每个自变量的 R 平方并显示在最后一列