stargazer, xtable, etc .: Ausgabe von benutzerdefinierten Standardfehlern mit Sternen, mit Ausnahme der Steuervariablen

Question

Das Bild unten zeigt, was ich aus R in LaTeX ausgeben möchte. Ich bin bereit, jeden Ansatz/jedes Paket zu verwenden, das funktioniert.

Wie in der Abbildung zu sehen ist, hat die gewünschte Ausgangs mehrere Modelle (die Spalten) und mehrere Regressionen (Zeilen).

Ich habe drei Hürden:

Erstens, wie kann ich Ausgang Standardfehler von zwei verschiedenen Modellen unter einer Punktschätzung? Keine der SE, die ich vorstellen möchte, sind die konventionellen SEs; Beide sind modifizierte Cluster-robuste SEs berechnet nach Ich führe die Regression mit benutzerdefinierten coeftest() Funktionen, und die Ausgabe eines coeftest Objekts.

Zweitens, wie kann ich die Sterne präsentieren? Ich habe einen Workaround in R entwickelt, um die Punktschätzung mit den beiden darunter liegenden SE-Berechnungen auszugeben, aber nicht mit den Sternen, die automatisch übertragen werden, wie dies etwa bei xtable oder stargazer der Fall ist.

Drittens möchte ich nur die Punktschätzung und Standardfehler für die Behandlungsvariable ausgeben. Wie Sie am Ende der Tabelle sehen können, gibt es in den Modellen (2) & (4) Steuervariablen, aber ich möchte keine weiteren Informationen über sie anzeigen.

Auch ist es erwähnenswert, dass mein Ausgang keine lm Objekt ist, sondern ein coeftest Objekt, das ist stargazer -kompatible aber nicht xtable -kompatible.

Answer 1

Werfen Sie einen Blick auf texreg, müssen Sie den Latext-Code ändern, den gewünschten Code einfügen.

library("texreg");library(lmtest);library("sandwich") 
library(nlme) 

m1<-lm(distance ~ age, data = Orthodont) 
coeftest(m1, vcov=sandwich) 

m2<-lm(distance ~ age + Sex, data = Orthodont) 
coeftest(m2, vcov=sandwich) 

test1=texreg(list(m1,m2),caption="Models", 
       label="Compmod",stars=c(0.05, 0.01)) 
#cat(test1) 

model.1 <- lme(distance ~ age, data = Orthodont, random = ~ 1) 
model.2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1) 
texreg(list(m1,m2,model.1, model.2),caption="Models", 
     label="Compmod", booktabs = TRUE, dcolumn = TRUE)