我确实用库(三明治)生成强大的标准错误 e.g。
library(sandwich)
cov.xxx <- vcovHC(xxx, type = "HC")
rob.std.err.xxx <- sqrt(diag(cov.xxx))
然后我想将它整合到一个Latex导出中并行显示几个模型。这很有效&#34;嗯&#34;。我用
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&#34;问题&#34;。但正如评论中所提到的,存在一种解决方案。
提前感谢您的支持。