library(sandwich)
library(stargazer)
library(lmtset)
#回归
reg<-lm(y~x,data=df)
#求稳健标准误
rse1<-vcovHAC(reg)
rse2<-sqrt(diag(rse1))
#工作台显示
coeftest(reg,.vcov=rse1)
#输出到word文档显示(第一列同方差结果,第二列稳健结果)
stargazer(reg,reg,se=list(NULL,rse2),out="reg.doc",type="html")