logo标签列表

将多列添加到我的texreg输出

rtexreg
9

我希望使用texreg制作一个带有分组列的表格。但是我只能看到分组行（groups）的选项。

这里有一个例子：

set.seed(01349)
DF <- data.frame(y = rnorm(100), x1A = rnorm(100), x2A = rnorm(100),
                 x1B = rnorm(100), x2B = rnorm(100))
regs <- lapply(paste0("x", 1:2, c("A", "A", "B", "B")), function(x)
          lm(paste0("y ~ ", x), data = DF))

这是我用普通的texreg尽可能接近的结果:

texreg(regs, custom.coef.names = c("Intercept", rep("x", 4)),
       custom.model.names = c("1", "2", "1", "2"))

使用LaTeX输出：

\begin{table}
\begin{center}
\begin{tabular}{l c c c c }
\hline
           & 1 & 2 & 1 & 2 \\
\hline
Intercept  & $-0.13$  & $-0.13$  & $-0.11$  & $-0.11$  \\
           & $(0.12)$ & $(0.12)$ & $(0.12)$ & $(0.12)$ \\
x          & $0.02$   & $0.07$   & $0.13$   & $-0.11$  \\
           & $(0.13)$ & $(0.12)$ & $(0.12)$ & $(0.13)$ \\
\hline
R$^2$      & 0.00     & 0.00     & 0.01     & 0.01     \\
Adj. R$^2$ & -0.01    & -0.01    & 0.00     & -0.00    \\
Num. obs.  & 100      & 100      & 100      & 100      \\
RMSE       & 1.18     & 1.17     & 1.17     & 1.17     \\
\hline
\multicolumn{5}{l}{\scriptsize{$^{***}p<0.001$, $^{**}p<0.01$, $^*p<0.05$}}
\end{tabular}
\caption{Statistical models}
\label{table:coefficients}
\end{center}
\end{table}

我希望增加一行（用%注释突出显示）：

\begin{table}
\begin{center}
\begin{tabular}{l c c c c }
\hline
%*************A HEADER LINE HERE*********************
 & \multicolumn{2}{c}{A} & \multicolumn{2}{c}{B} \\ %
%****************************************************
           & 1 & 2 & 1 & 2 \\
\hline
Intercept  & $-0.13$  & $-0.13$  & $-0.11$  & $-0.11$  \\
           & $(0.12)$ & $(0.12)$ & $(0.12)$ & $(0.12)$ \\
x          & $0.02$   & $0.07$   & $0.13$   & $-0.11$  \\
           & $(0.13)$ & $(0.12)$ & $(0.12)$ & $(0.13)$ \\
\hline
R$^2$      & 0.00     & 0.00     & 0.01     & 0.01     \\
Adj. R$^2$ & -0.01    & -0.01    & 0.00     & -0.00    \\
Num. obs.  & 100      & 100      & 100      & 100      \\
RMSE       & 1.18     & 1.17     & 1.17     & 1.17     \\
\hline
\multicolumn{5}{l}{\scriptsize{$^{***}p<0.001$, $^{**}p<0.01$, $^*p<0.05$}}
\end{tabular}
\caption{Statistical models}
\label{table:coefficients}
\end{center}
\end{table}

我是否有所遗漏，还是没有内置的方法可以做到这一点？

我的解决办法是：

x <- capture.output(texreg(
  regs, custom.coef.names = c("Intercept", rep("x", 4)),
  custom.model.names = c("1", "2", "1", "2")))

x[6] <- paste0("& \\multicolumn{2}{c}{A} & \\multicolumn{2}{c}{B} \\\\ \n", x[6])

cat(x, sep = "\n")

但是这显然有点像用胶带修补。
2个回答
6
这可能有些晚，但仍有用。 texreg的新版本（1.36.28）刚在GitHub上发布（还未在CRAN上发布，截至2020年5月22日）。它添加了选项custom.header，用于创建回归组。
应用于您的示例，它将创建：
library(texreg)

set.seed(01349)
DF <- data.frame(y = rnorm(100), x1A = rnorm(100), x2A = rnorm(100),
                 x1B = rnorm(100), x2B = rnorm(100))
regs <- lapply(paste0("x", 1:2, c("A", "A", "B", "B")), function(x)
  lm(paste0("y ~ ", x), data = DF))

screenreg(
  regs, 
  custom.header = list("A" = 1:2, "B" = 3:4),
  custom.coef.names = c("Intercept", rep("x", 4)),
  custom.model.names = c("1", "2", "1", "2"),
)

=============================================
                   A                 B       
           ----------------  ----------------
           1        2        1        2      
---------------------------------------------
Intercept   -0.13    -0.13    -0.11    -0.11 
            (0.12)   (0.12)   (0.12)   (0.12)
x            0.02     0.07     0.13    -0.11 
            (0.13)   (0.12)   (0.12)   (0.13)
---------------------------------------------
R^2          0.00     0.00     0.01     0.01 
Adj. R^2    -0.01    -0.01     0.00    -0.00 
Num. obs.  100      100      100      100    
=============================================
*** p < 0.001; ** p < 0.01; * p < 0.05

我使用了screenreg()来更轻松地显示输出结果，但它也可以与texreg()一起使用。

0

您可以通过在控制台中输入代码或从包texreg的github网站上阅读function的代码。

texreg
function (l, file = NULL, single.row = FALSE, stars = c(0.001, 
    0.01, 0.05), custom.model.names = NULL,...

这是4个模型表的texreg输出：

\begin{table}
\begin{center}
\begin{tabular}{l c c c c }
\hline
 & Model 1 & Model 2 & Model 3 & Model 4 \\

我查看了代码行 469，这是表格的开头：
 string <- paste0(string, "\\begin{tabular}{", coldef, 
            "}", linesep)

然后我加入了一些自己的修改：

string <- paste0(string, "\\begin{tabular}{", coldef, 
                 "}", linesep)
## Additions
string <- paste0(string, "\\\\[-1.8ex]\\hline", linesep)
string <- paste0(string, "& \\multicolumn{", length(l), 
                 "}{c}{\\textit{Dependent variable:}} \\\\", linesep)
string <- paste0(string, "\\cline{2-",length(modnames), "}", linesep)
string <- paste0(string, "\\\\[-1.8ex] & \\multicolumn{", length(l), 
                 "}{c}{", dep.var, "} \\\\", linesep)

然后，将该函数另存为不同的名称，例如：

texreg2 <- function (l, file = NULL, single.row = FALSE, ...)

现在，该函数需要从软件包中调用内部函数，因此您需要将自定义函数附加到软件包的命名空间中。非常简单：
environment(texreg2) <- asNamespace('texreg')

现在你可以调用你的新函数了。我的添加包括三行代码和一个依赖变量的名称，类似于stargazer。
texreg2(out, dep.var = "Normalize Citation Score")

\begin{table}
\begin{center}
\begin{tabular}{l c c c c }
\\[-1.8ex]\hline
& \multicolumn{4}{c}{\textit{Dependent variable:}} \\
\cline{2-5}
\\[-1.8ex] & \multicolumn{4}{c}{Normalize Citation Score} \\
\hline
 & Model 1 & Model 2 & Model 3 & Model 4 \\

最后，如果您不喜欢这种方法，您可以使用正则表达式来操作输出，查看此问题。
1
实际上这就是我一直以来添加“多列”所采用的方式。我发现这种方法有点太个人化了，没有必要将其编写成一个函数，而是逐个处理每种情况。感谢您的回复。 - MichaelChirico
网页内容由stack overflow 提供, 点击上面的
可以查看英文原文，
原文链接