stargazer, xtable и т. д .: вывод пользовательских стандартных ошибок со звездочками, исключая контрольные переменные

На изображении ниже показано, что я пытаюсь вывести из R в LaTeX. Я готов использовать любой подход / пакет, который работает.

Желаемый вывод TeX

Как вы можете видеть на изображении, желаемый результат имеет несколько моделей (столбцы) и несколько регрессий (строки).

У меня три препятствия:

Во-первых, как я могу вывести стандартные ошибки двух разных моделей ниже одной точечной оценки? Ни один из SE, который я пытаюсь представить, не является обычным SE; оба являются модифицированными кластерно-устойчивыми SE, рассчитанными после, я запускаю регрессию, используя пользовательские функции coeftest() и выводя объект coeftest.

Во-вторых, как я могу представить звезды? Я разработал обходной путь в R для вывода точечной оценки с двумя вычислениями SE под ней, но не с автоматическим переносом звезд, как что-то вроде xtable или stargazer.

В-третьих, я хотел бы вывести только точечную оценку и стандартные ошибки для переменной обработки. Как вы можете видеть внизу таблицы, в моделях (2) и (4) есть управляющие переменные, но я не хочу отображать какую-либо дополнительную информацию о них.

Также стоит отметить, что мои выходные данные не являются объектом lm, а вместо этого являются объектом coeftest, который совместим с stargazer, но не совместим с xtable.


person Dr. Beeblebrox    schedule 09.05.2016    source источник


Ответы (1)


Взгляните на texreg, вам нужно будет изменить латекстный код, вставив нужный код.

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)

введите здесь описание изображения

person Robert    schedule 09.05.2016