2015-04-06 63 views
0

在報告迴歸結果時,有時包含統計顯着性恆星在常數項中包含統計顯着性星。是否有可能配置stargazer讓逆變器保持恆星,但不能保持恆定期?Stargazer:省略僅用於常量的恆星

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

回答

0

基本上,答案竟然是使用stargazerp的說法。從那裏開始,我只需要編寫一系列(函數)函數,它們列出了迴歸擬合列表,並返回了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))) 

就像一個魅力。

0

你可以交替使用broom包到擬合結果轉換爲數據幀,然後添加星星你的心臟的內容:

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

而且xtable包可以用來格式化對於LaTeX或其他。

+0

stargazer的默認設置是乳膠。 OP應該說明他爲什麼不能編輯HTML或Latex輸出。 –

+0

手動編輯LaTeX輸出是最後一個解決方案(我可能最終使用),但它並不理想 - 主要是因爲表格不能直接從代碼中重現。 –