你用一下summary

fit1<-dynlm(log(chnimp)~log(chempi)+log(gas)+log(rtwex)+befile6+affile6+afdec6, data=ts_barium)
#用 OLS 模型
library(sandwich) #载入尼威-韦斯特估计所需要的包
library(lmtest) #载入 coeftest()所需的包
coeftest(fit1, vcov = NeweyWest(fit1, lag =4)) # 求 4 阶滞后的尼威-韦斯特估计量