18

本当に素朴な疑問!の線形回帰を多数実行y~xしており、コマンドで指定された標準エラー出力から手動で計算せずに、各回帰の分散を取得したいと考えていsummary.lmます。少し時間を節約するためです:-)。これを行うコマンドのアイデアはありますか? それとも、自分でそれを行うために関数を書く必要がありますか?

m<-lm(Alopecurus.geniculatus~Year)
> summary(m)

Call:
lm(formula = Alopecurus.geniculatus ~ Year)

Residuals:
    Min      1Q  Median      3Q     Max 
-19.374  -8.667  -2.094   9.601  21.832 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)  
(Intercept) 700.3921   302.2936   2.317   0.0275 *
Year         -0.2757     0.1530  -1.802   0.0817 .
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

Residual standard error: 11.45 on 30 degrees of freedom
  (15 observations deleted due to missingness)
Multiple R-squared: 0.09762,    Adjusted R-squared: 0.06754 
F-statistic: 3.246 on 1 and 30 DF,  p-value: 0.08168 

だから私は標準エラー出力を取得し、手動で計算せずに分散出力を取得したいと思っていました...

4

3 に答える 3

25

あなたが何の分散を望んでいるのかわかりません。

残差分散が必要な場合は、次のとおり(summary(m)$sigma)**2です。

勾配の分散が必要な場合は(summary(m)$coefficients[2,2])**2、 、またはvcov(m)[2,2]です。

于 2013-02-19T15:18:24.050 に答える
13

vcov(m)

係数の共分散行列–対角線上の分散を与えます。

于 2013-02-19T18:42:01.213 に答える
11

係数推定値の標準誤差を参照している場合、答えは

summary(m)$coef[,2] 

推定残差分散を参照している場合、それは

summary(m)$sigma

アクセスできるその他の情報についてはnames( summary(m) )、 とを入力します。names(m)

于 2013-02-19T15:17:52.750 に答える