所有
我试图找出计算下面回归因子的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>
答案 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
源代码将提供它如何做的细节,这可能与我们在这里展示的不完全相同,但除了数值误差外,它们是等效的。