1

我确实使用库(三明治)生成了强大的标准错误,例如

library(sandwich)
cov.xxx <- vcovHC(xxx, type = "HC")
rob.std.err.xxx <- sqrt(diag(cov.xxx))

然后我想将它集成到一个 Latex 导出中,同时呈现几个模型。这“很好”。我用

library(texreg)
      texreg(list(x, xx, xxx, xxxx), 
              custom.model.names = c("1", "2", "3", "4"),
              custom.coef.names = c("Constant", 
              "A", "B", "C",  "D",  "E", "F", "G", "H","I", "J",  "K",  "L", 
              "Wrong Sequenced Element",
              "N", "O", "P", "Q", "R", "S", "T", "U",  "V", "W", "X"),
         groups = list("Var1" = 1:5, 
                       "Var2" = 6:14, 
                       "Var3" = 15:19, "Var4" = 20:25),
         sideways = T,
         single.row = T, 
        override.se = list(c(0.1,0.1, 0.1, 0.1, 0.1),
                           c(0.1, 0.1, 0.1, 0.1, 0.1,
                             0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1),
                           c(0.1, 0.1, 0.1, 0.1, 0.1,
                             0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 
                             0.1, 0.1, 0.1, 0.1, 0.1),
                           **rob.std.err.xxx   #values from robust regression**
                  ), 
         caption= "zzz", 
         return.string = TRUE, use.packages = FALSE , booktabs = TRUE, dcolumn = TRUE, digits = 3,
            file =  paste0(tablesdir, "xxxx.tex"))

但是 rob.std.err.xxx 中的属性顺序与我的 LM 模型不对应。而且我无法弄清楚如何正确更改 SE 覆盖。如果有人可以就 p 值的相同程序提供建议,将不胜感激!

我所需要的只是更改 rob.std.err.xxx 中的序列,然后将其导出到上面所示的目标结构中。

这个问题主要是指一个texreg“问题”。但正如评论中提到的,存在一个解决方案。

提前感谢您的任何支持。

4

0 回答 0