Linear regression 具有NeweyWest误差项的OLS回归

我需要用OLS方法估计一个线性回归。因为我假设误差项是相关的，所以我想解释误差项中的异方差和自相关。我写了一段代码，但由于我对R一窍不通，我怀疑它是否正确

到目前为止，我的代码是：

library(sandwich)
library(dynlm)
library(lmtest)
myregression <- dynlm(formula = x ~ y + z)
coeftest(myregression, vcov = NeweyWest(myregression))

库（三明治）
图书馆（dynlm）
图书馆（lmtest）
我的回归