使用R语言从线性模型中提取R平方值

Question

我知道使用 summary 可以帮助我手动计算，但是我需要计算大量的R平方值。因此，我需要电脑为我提取它。以下是一个简单的例子：

library(alr3)
M.lm=lm(MaxSalary~Score,data=salarygov)
#Here you will see the R square value
summary(M.lm)

我该怎么做？

- Zander

Answer 1

您在摘要中看到的R-squared，调整后的R-squared和所有其他值都可以从摘要对象中访问。您可以使用str(summary(M.lm))查看所有内容：

> str(summary(M.lm))  # Truncated output...
List of 11
 $ call         : language lm(formula = MaxSalary ~ Score, data = salarygov)
 $ terms        :Classes 'terms', 'formula' length 3 MaxSalary ~ Score
 ...
 $ residuals    : Named num [1:495] -232.3 -132.6 37.9 114.3 232.3 ...
 $ coefficients : num [1:2, 1:4] 295.274 5.76 62.012 0.123 4.762 ...
 $ aliased      : Named logi [1:2] FALSE FALSE
 $ sigma        : num 507
 $ df           : int [1:3] 2 493 2
 $ r.squared    : num 0.817
 $ adj.r.squared: num 0.816
 $ fstatistic   : Named num [1:3] 2194 1 493
 $ cov.unscaled : num [1:2, 1:2] 1.50e-02 -2.76e-05 -2.76e-05 5.88e-08

要获得R平方值，请键入summary(M.lm)$r.squared或summary(M.lm)$adj.r.squared。

- Andrew

---

当只有一个左侧变量时，这似乎才有效。 - John

当有两个左侧变量时，摘要变量仍应全部显示。请参见str(summary(lm(price + carat ~ x + y + z, data=ggplot2::diamonds)))。 - Andrew

Answer 2

取决于您感兴趣的是哪一个:

 # adjusted R²
 summary(M.lm)$adj.r.squared
 # R²
 summary(M.lm)$r.squared

- droopy

Answer 3

如果只有一个预测变量，您可以简单地使用 cor(salarygov$MaxSalary ,salarygov$Score)^2。或者，summary(M.lm)$r.squared。

- Roland