3

考虑以下示例

inds <- c('var1','','var2','')
model1 <- c(10.2,0.00,0.02,0.3)
model2 <- c(11.2,0.01,0.02,0.023)

df = df=data.frame(inds,model1,model2)
df
 inds model1 model2
 var1  10.20 11.200
        0.00  0.010
 var2   0.02  0.020
        0.30  0.023

在这里，您可以获得带有系数和 P 值的自定义回归模型的输出（如果需要，我实际上可以显示任何其他统计数据，例如系数的标准误差）。

有两个变量，var1和var2。

例如，在模型 1 中，var1系数为10.2，P 值为 ，0.00而var2系数为0.02，P 值为0.30。

是否有一个包可以自动处理这些（自定义）表格，并且可以创建一个带有星号的整洁 Latex 表格以表示重要性？

谢谢！

4

1 回答 1

5

这是一个使用texreg.

请注意texreg>= 1.36.18 是必需的。

您在数据框中提供的信息（系数和 p 值）可以在数据框中以任意方式排列。因此，我们需要编写代码，从数据框中的适当位置选择这些数据并使用它们来创建texreg对象。当您请求一个通用的（并且可能是可重用的）解决方案时，我们应该将代码包装在一个可重用的函数中。我会调用这个函数extractFromDataFrame。所以这是一个函数，它从数据框中提取信息并texreg为不同的模型创建一个对象列表：

require("texreg")

extractFromDataFrame <- function (dataFrame) {
  coef.row.indices <- seq(1, nrow(dataFrame) - 1, 2)
  pval.row.indices <- seq(2, nrow(dataFrame), 2)
  texregObjects <- list()
  for (i in 2:ncol(dataFrame)) {
    coefs <- dataFrame[coef.row.indices, i]
    coefnames <- as.character(dataFrame[coef.row.indices, 1])
    pvalues <- dataFrame[pval.row.indices, i]
    tr <- createTexreg(coef = coefs, coef.names = coefnames, pvalues = pvalues)
    texregObjects[i - 1] <- list(tr)
  }
  return(texregObjects)
}

在这个函数中，我们首先定义系数存储在数据帧的哪些行中，以及 p 值存储在哪些行中。然后我们创建了一个存储texreg对象的空列表。我们遍历所有列，但第一列仅包含标签。在每个模型列中，我们保存系数、它们的名称和 p 值，然后将它们交给createTexreg构造函数，这是一个texreg根据数据为我们创建对象的函数。我们将texreg对象添加到列表中。最后，我们返回texreg对象列表。

我们现在可以将该函数应用于任何看起来像问题中提供的数据框，具有任意数量的列 (> 1)。在这种情况下，在将函数应用到df对象之后，如果我们想确保我们做的一切都正确，我们可能想要打印列表的内容：

tr <- extractFromDataFrame(df)
tr

事实上，结果包含相关数据：

[[1]]

No standard errors were defined for this texreg object.
No decimal places were defined for the GOF statistics.

     coef.   p
var1 10.20 0.0
var2  0.02 0.3

No GOF block defined.

[[2]]

No standard errors were defined for this texreg object.
No decimal places were defined for the GOF statistics.

     coef.     p
var1 11.20 0.010
var2  0.02 0.023

No GOF block defined.

现在我们可以简单地将texreg对象列表交给screenreg，例如 ，screenreg(tr)结果如下：

========================
      Model 1    Model 2
------------------------
var1  10.20 ***  11.20 *
var2   0.02       0.02 *
========================
*** p < 0.001, ** p < 0.01, * p < 0.05

或htmlreg用于创建 HTML 表格。或者，按照原始问题中的要求，texreg用于创建 LaTeX 表。的输出texreg(tr, single.row = TRUE)如下所示：

\begin{table}
\begin{center}
\begin{tabular}{l c c }
\hline
 & Model 1 & Model 2 \\
\hline
var1 & $10.20^{***}$ & $11.20^{*}$ \\
var2 & $0.02$        & $0.02^{*}$  \\
\hline
\multicolumn{3}{l}{\scriptsize{$^{***}p<0.001$, $^{**}p<0.01$, $^*p<0.05$}}
\end{tabular}
\caption{Statistical models}
\label{table:coefficients}
\end{center}
\end{table}

可以修改此解决方案以适应标准误差、置信区间或拟合优度统计量。

texreg可以使用各种参数来自定义输出，例如，包括使用booktabs包或十进制对齐方式dcolumn。

请注意，您不应调用数据框df，因为该对象名称已在stats包中定义。

于 2016-10-22T23:15:47.823 回答