如何获得每个预测器对 lm 中最终回归预测的贡献答案

【问题标题】：How can I get the contribution by each predictor to the final regression prediction in lm如何获得每个预测器对 lm 中最终回归预测的贡献
【发布时间】：2017-05-09 09:16:11
【问题描述】：

当我使用rlm 或lm 时使用R 我想得到模型每个预测变量的贡献。

当我认为它们不在lm 对象中时，会出现问题

下面是样本数据（我正在寻找一种可以推广到任意数量的预测变量的方法）

样本数据：

set.seed(1)                                              
y <- rnorm(10)                                           
m <- data.frame(v1=rnorm(10), v2=rnorm(10), v3=rnorm(10))
lmObj <- lm(formula=y~0+v1*v3+v2*v3, data=m)               
betaHat <- coefficients(lmObj)                           

betaHat
      v1       v3       v2    v1:v3    v3:v2
 0.03455 -0.50224 -0.57745  0.58905 -0.65592

# How do I get the data.frame or matrix with columns (v1,v3,v2,v1:v3,v3:v2) 
# worth [M$v1*v1, ... , (M$v3*M$v2)*v3:v2]

【问题讨论】：

我建议考虑relaimpo 包：cran.r-project.org/web/packages/relaimpo/relaimpo.pdf
@MarcoSandri 我想要矩阵，那里没有分析，纯粹是如何从 lm 对象中获取矩阵
@LyzandeR：不，这不是我想要的 $[M$v1*v1, ... , (M$v3*M$v2)*v3:v2]$

标签： r dataframe lm coefficients

【解决方案1】：

我认为通过“贡献”你想要解释每个术语的方差（ANOVA 表有帮助），而实际上你想要术语预测：

predict(lmObj, type = "terms")

见?predict.lm。

【讨论】：

是的，这正是我想要的，太糟糕了，它不适用于rlm，所以对于rlm，我的回答是好的（见predict.lmRob）

【解决方案2】：

其实我是从lm自己那里得到的，诀窍是要求x=TRUE

lmObj <- lm(formula=y~0+v1*v3+v2*v3, data=m, x=TRUE)               
lmObj$x %*% diag(lmObj$coefficients) 
         [,1]     [,2]     [,3]     [,4]     [,5]
1   0.0522305 -0.68238 -0.53066  1.20993 -0.81898
2   0.0134687  0.05162 -0.45164 -0.02360  0.05273
3  -0.0214632 -0.19470 -0.04306 -0.14187 -0.01896
4  -0.0765156  0.02702  1.14875  0.07019 -0.07021
5   0.0388652  0.69161 -0.35792 -0.91250  0.55985
6  -0.0015524  0.20843  0.03241  0.01098 -0.01528
7  -0.0005594  0.19803  0.08996  0.00376 -0.04029
8   0.0326086  0.02979  0.84928 -0.03298 -0.05722
9   0.0283723 -0.55248  0.27611  0.53213  0.34500
10  0.0205187 -0.38330 -0.24134  0.26699 -0.20921

【讨论】：