在 R 中没有 lm 对象的情况下计算 Newey-West 标准误差

发布于 2024-10-07 12:30:06 字数 611 浏览 4 评论 0原文

更新 - 我关闭了这个问题并 发布在 crossvalidated.com

我发现了一些有关使用 sandwich 包和 NeweyWest() 函数的好信息找出异方差自相关一致 (HAC) 标准误差。

NeweyWest() 只接受 lm 对象。

> library(sandwich)
> NeweyWest(rnorm(100))
Error in UseMethod("estfun") : 
  no applicable method for 'estfun' applied to an object of class "c('double', 'numeric')"
> 

我经常得到与线性回归无关的回报向量,我想找到 HAC 标准误差。有什么想法吗?我应该自己写吗?谢谢!

Update -- I closed this question and posted on crossvalidated.com.

I have found some good information on using the sandwich package and the NeweyWest() function to find heteroskedastic autocorrelation consistent (HAC) standard errors.

But NeweyWest() only takes lm objects.

> library(sandwich)
> NeweyWest(rnorm(100))
Error in UseMethod("estfun") : 
  no applicable method for 'estfun' applied to an object of class "c('double', 'numeric')"
> 

I frequently get vectors of returns unassociated with a linear regression for which I would like to find HAC standard errors. Any ideas? Should I write my own? Thanks!

如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。

扫码二维码加入Web技术交流群

发布评论

需要 登录 才能够评论, 你可以免费 注册 一个本站的账号。

评论(1

二智少女猫性小仙女 2024-10-14 12:30:06

有一点小小的误会。我在考虑残差,但你问的是平均值的标准误差。这可以通过对截距建模向量来轻松获得,或者:

NeweyWest(lm(rnorm(100)~1))

对于标准差:

x <- rnorm(100)
NeweyWest(lm(x~1))*length(x)

抱歉造成误解,我的错。

There's been a slight misunderstanding. I was thinking in terms of residuals, but what you asked is the standard error of the mean. That's easily obtained by modelling your vector against the intercept, or :

NeweyWest(lm(rnorm(100)~1))

For the standard deviation :

x <- rnorm(100)
NeweyWest(lm(x~1))*length(x)

Sorry for the misunderstanding, my bad.

~没有更多了~
我们使用 Cookies 和其他技术来定制您的体验包括您的登录状态等。通过阅读我们的 隐私政策 了解更多相关信息。 单击 接受 或继续使用网站,即表示您同意使用 Cookies 和您的相关数据。
原文