观星者：仅为常数省略星星

Question

在报告回归结果时，有时候在常数项中包含统计显着性星是很俗气的。是否可以配置stargazer以保留回归量的星号，但不能用于常数项？

fit <- lm(rating ~ complaints, data=attitude)
stargazer(fit)

Answer 1

您也可以使用broom包将拟合结果转换为数据框，然后将星星添加到您心中的内容中：

library("broom")
mod <- lm(mpg ~ wt + qsec, data = mtcars)
DF <- tidy(mod)
DF$stars <- c("", "***", "***") # inspect and add manually, or automate

xtable包可用于为LaTeX或其他任何方式格式化它。

Answer 2

基本上，答案是使用stargazer的{​​{1}}参数。从那里开始，我只需编写一系列（一系列）函数，这些函数采用回归拟合列表并返回p值向量列表。然后我手动将截距的p值更改为1，并且在截距上没有任何俗气的星星。此外，它没有手动LaTeX编辑的可重复性！

p

然后在观星者的电话中：

commarobust <- function(fit){
  require(sandwich)
  require(lmtest)
  coeftest(fit,vcovHC(fit, type="HC2"))
}

getrobustps <- function(fit){
  robustfit <- commarobust(fit)
  ps <- robustfit[,4]
  ps["(Intercept)"] <- 1
  return(ps)
}

makerobustpslist <- function(fitlist){
  return(lapply(fitlist, FUN=getrobustps) )
}

像魅力一样。