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 回答