觀星者：只為恆星省略星號

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的p參數。 從那里開始，我只需要編寫一個（系列）函數即可，該函數采用了一系列的回歸擬合並返回了一個p值向量的列表。 然后，我手動將截距的p值更改為1，並且預先確定，截距上沒有發粘的星星。 另外，它無需手動LaTeX編輯就可重現！

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) )
}

然后在觀星者電話中：

stargazer(fit_1, fit_2, fit_3, fit_4, fit_5, 
          p=makerobustpslist(list(fit_1, fit_2, fit_3, fit_4, fit_5)))

奇跡般有效。