web-dev-qa-db-ja.com

線形モデルでRを使用してR-2乗値を抽出する

summaryを使用すると、これを手動で行うことができますが、R平方値のトンを計算する必要があります。したがって、コンピューターを抽出する必要があります。以下に簡単な例を示します。

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

どうやってやるの?

35
Zander

Rの2乗、調整されたRの2乗、およびサマリーに表示されるその他のすべての値には、サマリーオブジェクト内からアクセスできます。 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 2乗値を取得するには、summary(M.lm)$r.squaredまたはsummary(M.lm)$adj.r.squaredと入力します

65
Andrew

1つの予測子では、cor(salarygov$MaxSalary ,salarygov$Score)^2を使用できます。または、summary(M.lm)$r.squared

7
Roland

どちらに興味があるかによります:

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