r - stargazer, xtable, etc.: outputting custom standard errors with stars, excluding control variables -

the image below demonstrates trying output r latex. i'm willing use approach/package works.

as can see in image, desired output has multiple models (the columns) , multiple regressions (the rows).

i have 3 hurdles:

first, how can output standard errors 2 different models below 1 point estimate? neither of ses seek present conventional ses; both modified cluster-robust ses calculated after run regression using custom coeftest() functions, , outputting coeftest object.

second, how can present stars? have devised workaround in r output point estimate 2 se calculations below it, not stars automatically transferred xtable or stargazer does.

third, output point estimate , standard errors treatment variable. can see @ bottom of table, in models (2) & (4) there control variables, not want display further information them.

also, worth noting output not lm object, instead coeftest object, stargazer-compatible not xtable-compatible.

take @ texreg, need modify latext code, inserting code want.

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)