在stargazer表中显示自由度

Question

在使用 Sweave 和 R 构建文档时，我会利用 stargazer 库来制作表格。

在使用 stargazer 时，是否有机制可以显示使用 glm 构建的模型中与残差偏差相关的自由度？

最小代码：

library(stargazer)
set.seed(1234)
data <- data.frame(x=1:10)
data$y <- data$x + rnorm(10, 0, 0.2)
model <- glm(y~x, data=data, family=gaussian)
summary(model)
stargazer(model,title="A test", align=T,label="Tab:test",style="all2")

生成的stargazer表将包含观测值、对数似然、AIC、残差偏差和零偏差，但不包括自由度。我可以计算自由度，但认为这些应该直接显示出来。另请参见：https://sites.google.com/site/marekhlavac/stargazer 更新#1：

感谢Marek的回复。为了让其他人受益，以下是让您绕过此问题的过程：

1. 从http://cran.r-project.org/src/contrib/Archive/stargazer/获取版本4.0（而不是4.5-我会回到这个问题）。
2. 在R下的包目录结构中，根据下面答案中的说明编辑“stargazer-internal.R”。
3. 确保您的R会话中没有加载库
4. 确保您已删除任何现有的stargazer lib
5. 安装编辑后的stargazer软件包。
6. 像往常一样重新加载库并编译。

以下是命令：

detach("package:stargazer", unload=TRUE)
remove.packages("stargazer")

从命令行操作：

R CMD INSTALL -l <path to library directory> stargazer

最终（假设你手头有几个模型），

library(stargazer)
stargazer(model6,model7,model8, title="Logistic model summary",align=T,label="Tab:logmod1",  font.size="footnotesize", style="all2")

结果：

% Table created by stargazer v.4.0 by Marek Hlavac, Harvard University. E-mail: hlavac at fas.harvard.edu
% Date and time: Tue, Sep 24, 2013 - 17:17:17
% Requires LaTeX packages: dcolumn 
\begin{table}[!htbp] \centering 
  \caption{Logistic model summary} 
  \label{Tab:logmod1} 
\footnotesize 
\begin{tabular}{@{\extracolsep{5pt}}lD{.}{.}{-3} D{.}{.}{-3} D{.}{.}{-3} } 
\\[-1.8ex]\hline 
\hline \\[-1.8ex] 
 & \multicolumn{3}{c}{\textit{Dependent variable:}} \\ 
\cline{2-4} 
\\[-1.8ex] & \multicolumn{3}{c}{whalesighted} \\ 
\\[-1.8ex] & \multicolumn{1}{c}{\textit{logistic}} & \multicolumn{1}{c}{\textit{probit}} & \multicolumn{1}{c}{\textit{glm: binomial}} \\ 
 & \multicolumn{1}{c}{\textit{}} & \multicolumn{1}{c}{\textit{}} & \multicolumn{1}{c}{\textit{link = cloglog}} \\ 
\\[-1.8ex] & \multicolumn{1}{c}{(1)} & \multicolumn{1}{c}{(2)} & \multicolumn{1}{c}{(3)}\\ 
\hline \\[-1.8ex] 
 visibility & 0.392^{***} & 0.226^{***} & 0.216^{***} \\ 
  & (0.051) & (0.027) & (0.026) \\ 
  Constant & -1.251^{***} & -0.745^{***} & -1.149^{***} \\ 
  & (0.246) & (0.144) & (0.182) \\ 
 \hline \\[-1.8ex] 
Observations & \multicolumn{1}{c}{232} & \multicolumn{1}{c}{232} & \multicolumn{1}{c}{232} \\ 
Log Likelihood & \multicolumn{1}{c}{-110.485} & \multicolumn{1}{c}{-110.888} & \multicolumn{1}{c}{-112.694} \\ 
Akaike Inf. Crit. & \multicolumn{1}{c}{224.970} & \multicolumn{1}{c}{225.775} & \multicolumn{1}{c}{229.388} \\ 
Residual Deviance (df = 230) & \multicolumn{1}{c}{220.970} & \multicolumn{1}{c}{221.775} & \multicolumn{1}{c}{225.388} \\ 
Null Deviance (df = 231) & \multicolumn{1}{c}{310.759} & \multicolumn{1}{c}{310.759} & \multicolumn{1}{c}{310.759} \\ 
\hline 
\hline \\[-1.8ex] 
\textit{Note:}  & \multicolumn{3}{r}{$^{*}$p$<$0.1; $^{**}$p$<$0.05; $^{***}$p$<$0.01} \\ 
\normalsize 
\end{tabular} 
\end{table} 

回到我在实施基于4.5代码的解决方法时遇到的错误。实际上，当我从mac二进制文件（版本4.5.1）安装时，我会遇到相同的错误（http://cran.r-project.org/web/packages/stargazer/index.html），只需尝试使用stargazer，如下所示。

> install.packages("stargazer")
trying URL 'http://cran.ms.unimelb.edu.au/bin/macosx/contrib/3.0/stargazer_4.5.1.tgz'
Content type 'application/x-tar' length 332917 bytes (325 Kb)
opened URL
==================================================
downloaded 325 Kb

> stargazer(model6,model7,model8,      
+           title="Logistic model summary",
+           align=T,
+           label="Tab:logmod1", 
+           font.size="footnotesize",
+           style="all2")
Error in `rownames<-`(`*tmp*`, value = "visibility") : 
  length of 'dimnames' [1] not equal to array extent

Marek，供您参考，我将向您发送结果的traceback()。祝好。

- t-student

---

我不是stargazer的用户，但如果有最少量的代码提供，我通常会“试一试”解决问题。 - IRTFM

感谢您的详细跟进。我会检查您的问题，并在下一个版本中进行必要的修复。请注意，在不实施解决方法的情况下，4.5.1可以正常工作（只是不会默认输出自由度以获取空值和残差偏差）。 - user1953965

Answer 1

Stargazer的作者在这里。看起来该软件包默认不会输出残差和零偏差的自由度。我将考虑在未来版本中更改默认设置。

现在，您可以尝试使用源包（来自CRAN），并修改stargazer-internal.R中的函数.adjust.settings.style以包含以下内容作为快速修复：

   if (style == "all") {
      .format.table.parts <<- c("=!","dependent variable label","dependent variables","models","columns","numbers","-","coefficients","-","omit","-","additional","N","R-squared","adjusted R-squared","max R-squared","log likelihood","sigma2","theta(se)*(p)", "SER(df)","F statistic(df)*(p)","chi2(df)*(p)","Wald(df)*(p)","LR(df)*(p)","logrank(df)*(p)","AIC","BIC","UBRE","rho(se)*(p)","Mills(se)*(p)","residual deviance(df)*","null deviance(df)*","=!","notes")  
      .format.coefficient.table.parts <<- c("variable name","coefficient*","standard error","t-stat","p-value")  
    }

    else if (style == "all2") {
      .format.table.parts <<- c("=!","dependent variable label","dependent variables","models","columns","numbers","-","coefficients","-","omit","-","additional","N","R-squared","adjusted R-squared","max R-squared","log likelihood","sigma2","theta(se)*(p)", "SER(df)","F statistic(df)*(p)","chi2(df)*(p)","Wald(df)*(p)","LR(df)*(p)","logrank(df)*(p)","AIC","BIC","UBRE","rho(se)*(p)","Mills(se)*(p)","residual deviance(df)*","null deviance(df)*","=!","notes")  
      .format.coefficient.table.parts <<- c("variable name","coefficient*","standard error")  
    }

请注意，这里唯一的更改是我在“残差偏差”和“空模型偏差”中添加了“(df)*”。

- user1953965