2012-03-06 7 views
3

からエラー値を取得します。私が使って解く線形回帰問題を抱えている線形回帰関数LM

しかし、私は標準誤差(フィッティングエラー)を取得したいと思います。私は値を見ることができましたが、私はデータフレームの中にそれを格納するためにそれを得る方法を知らない。

私はsummary(m)を使用して値を取得し、結果はこのようなものです:

Call: 
lm(formula = value ~ mean, data = d) 

Residuals: 
    Min  1Q Median  3Q  Max 
-25.000 -15.909 -2.124 14.596 44.697 

Coefficients: 
       Estimate Std. Error t value Pr(>|t|)  
(Intercept) 2.500e+01 1.064e+00 23.49 <2e-16 *** 
mean  -1.759e-06 1.536e+00 0.00  1  
--- 
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

Residual standard error: 16.85 on 1298 degrees of freedom 
Multiple R-squared: 1.01e-15, Adjusted R-squared: -0.0007704 
F-statistic: 1.311e-12 on 1 and 1298 DF, p-value: 1 

ので、質問は:どのように私は、これらの値へのアクセスを得ることができます?

はあなたにresid(m)を使用して

答えて

8

関数summaryは、Rリストを返します。

##Generate some dummy data 
x = runif(10);y = runif(10) 
m = summary(lm(y ~ x)) 

通常のリスト構文を使用して、必要なものを抽出することができます。例えば、

m[[4]] 

は、モデルのデータフレームが

R> m[[4]] 
      Estimate Std. Error t value Pr(>|t|) 
(Intercept) 0.44265  0.2443 1.8123 0.1075 
x   0.07066  0.4460 0.1584 0.8781 

m[[6]]Residual standard error

R> m[[6]] 
[1] 0.2928 

coefficients(m)

+0

だから、それは私が欲しいエラー関数の結果ですか? – ifreak

+0

私はそうだと思いますが、おそらくあなたはそれを評価するより良い場所にいるでしょう。 – csgillespie

+0

だから、最初の部分の合計(0からaまで)と2番目の部分の和(aから長さまで)であるエラー関数を作ったと考えることができます。合計で我々は傍受と傾斜のいくつかの機能を持っています。この値はこの関数の結果でなければなりません。 – ifreak

2

アクセス残差に感謝します。

EDIT:コメントから、あなたはsum(resid(m)^2)が必要なようです。

+0

などの周りのいくつかの便利な機能がありますを返すに適合返します私は住んでみた(m)と私は多くの行のようになった、何これは意味ですか?私はフィッティングエラーが欲しい。 – ifreak

+0

フィッティングエラーが欲しいです。 だから、基本的にはエラーの値です。例えば、線形回帰は誤差関数を最小にすることに作用する。だから私はこのエラー関数の値が必要です。 – ifreak