r - MASS::lm.ridge 系数

Question

我正在尝试使用 lm.ridge 方法执行岭回归。我的问题是如何获得拟合模型的系数？我从调用中得到不同的结果:

• 模型$coef
• 系数(模型)

哪个是正确的？另外，为什么我从调用中得到不同的结果:

• coef(model) 并查看第一个系数，对比
• 系数(模型)[1]?

Answer 1

作为?lm.ridge说(在描述返回对象的 $coef 元素时)[强调]

coef: matrix of coefficients, one row for each value of ‘lambda’. Note that these are not on the original scale and are for use by the ‘coef’ method.

具体来说，这意味着 $coef元素不适用于最终用户(“如果你必须问......”)。 (如果您想查看 $coef 是如何翻译的，请检查 MASS:::coef.ridgelm 。)通常，最好使用访问器方法，例如 coef()。 ，当它存在时，而不是使用 $ 从返回对象的内部提取组件(或 @ 对于 S4 对象)——正是出于这个原因。包作者提供 coef()出于某种原因的方法......

我无法复制你的第二个问题。使用来自 ?lm.ridge 的模型, 除了打印的精度外，答案似乎完全相同......

> m1 <- lm.ridge(y ~ ., longley)
> coef(m1)
                            GNP    Unemployed  Armed.Forces    Population 
    <strong>2946.85636017</strong>    0.26352725    0.03648291    0.01116105   -1.73702984 
             Year      Employed 
      -1.41879853    0.23128785 
> coef(m1)[1]       
<strong>2946.856</strong>