在報告迴歸結果時,有時包含統計顯着性恆星在常數項中包含統計顯着性星。是否有可能配置stargazer
讓逆變器保持恆星,但不能保持恆定期?Stargazer:省略僅用於常量的恆星
fit <- lm(rating ~ complaints, data=attitude)
stargazer(fit)
在報告迴歸結果時,有時包含統計顯着性恆星在常數項中包含統計顯着性星。是否有可能配置stargazer
讓逆變器保持恆星,但不能保持恆定期?Stargazer:省略僅用於常量的恆星
fit <- lm(rating ~ complaints, data=attitude)
stargazer(fit)
基本上,答案竟然是使用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)))
就像一個魅力。
你可以交替使用broom
包到擬合結果轉換爲數據幀,然後添加星星你的心臟的內容:
library("broom")
mod <- lm(mpg ~ wt + qsec, data = mtcars)
DF <- tidy(mod)
DF$stars <- c("", "***", "***") # inspect and add manually, or automate
而且xtable
包可以用來格式化對於LaTeX或其他。
stargazer的默認設置是乳膠。 OP應該說明他爲什麼不能編輯HTML或Latex輸出。 –
手動編輯LaTeX輸出是最後一個解決方案(我可能最終使用),但它並不理想 - 主要是因爲表格不能直接從代碼中重現。 –