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