Я генерирую надежные стандартные ошибки с библиотекой (бутербродом), например.
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. Но, как упоминалось в комментарии, решение существует.
Заранее спасибо за любую поддержку.