标准误差的公式

时间:2015-11-08 18:32:03

标签: r lm

所有

我试图找出计算下面回归因子的std.error的公式以及如何使用mean和sd函数计算它。 (标准误= 2.015)。请帮忙。

感谢。 瑞克

<script src="//ajax.googleapis.com/ajax/libs/jquery/2.1.1/jquery.min.js"></script>

<script type="text/javascript" src="//cdnjs.cloudflare.com/ajax/libs/jqueryui/1.11.4/jquery-ui.min.js"></script>
<link rel="stylesheet" type="text/css" href="//cdnjs.cloudflare.com/ajax/libs/jqueryui/1.11.4/jquery-ui.min.css">

<div id="divd" class="draggable ui-widget-content common-drag">Test</div>
<button id="button">Button</button>

1 个答案:

答案 0 :(得分:3)

尝试其中任何一个来获取std错误:

sqrt(sum(resid(lm1)^2)/(length(factor1) - nlevels(factor1)))

sqrt(deviance(lm1)/(length(factor1) - nlevels(factor1)))

summary.lm(lm1)$sigma

library(broom); glance(lm1)$sigma

如果您想要系数的标准误差,那么如果se是上述任何一个,那么:

sqrt(diag(vcov(lm1)))

se * sqrt(diag(solve(crossprod(model.matrix(lm1)))))

se * sqrt(diag(summary.lm(lm1)$cov))

coef(summary(lm1))[, 2]

library(broom); tidy(lm1)$std.error

注意(1)因为问题没有使用set.seed将随机数生成器设置为已知状态,所以数据不可重现;(2)如评论中所提到的summary.lm源代码将提供它如何做的细节,这可能与我们在这里展示的不完全相同,但除了数值误差外,它们是等效的。