Question

我想知道如何使用以下 R 代码中的xtable函数获取LATEX表。

Block <- gl(8, 4)
A <- factor(c(0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,
          0,1,0,1,0,1,0,1,0,1,0,1))
B <- factor(c(0,0,1,1,0,0,1,1,0,1,0,1,1,0,1,0,0,0,1,1,
          0,0,1,1,0,0,1,1,0,0,1,1))
C <- factor(c(0,1,1,0,1,0,0,1,0,0,1,1,0,0,1,1,0,1,0,1,
          1,0,1,0,0,0,1,1,1,1,0,0))
Yield <- c(101, 373, 398, 291, 312, 106, 265, 450, 106, 306, 324, 449,
       272, 89, 407, 338, 87, 324, 279, 471, 323, 128, 423, 334,
       131, 103, 445, 437, 324, 361, 302, 272)
aovdat <- data.frame(Block, A, B, C, Yield)

summary(aov(Yield~Block+A*B+Error(A*Block), data=aovdat))


Error: A
  Df Sum Sq Mean Sq
A  1 3465.3  3465.3

Error: Block
      Df Sum Sq Mean Sq
Block  7   4499  642.71

Error: A:Block
          Df Sum Sq Mean Sq F value Pr(>F)
B          1  41616   41616  3.5354 0.1091
Residuals  6  70628   11771               

Error: Within
          Df Sum Sq Mean Sq F value  Pr(>F)  
B          1 119568  119568  7.3751 0.01673 *
A:B        1     28      28  0.0017 0.96734  
Residuals 14 226975   16213                  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

使用以下代码

xtable(summary(aov(Yield~Block+A*B+Error(A*Block), data=aovdat)))

产生以下错误

Error in rbind(deparse.level, ...) : 
  numbers of columns of arguments do not match

是否可以获得这样的LaTeX表？

SOV       Df Sum Sq Mean Sq F value  Pr(>F)  
A          1 3465.3  3465.3
Block      7   4499  642.71
B          1  41616   41616  3.5354 0.1091
Residuals  6  70628   11771               
B          1 119568  119568  7.3751 0.01673 *
A:B        1     28      28  0.0017 0.96734  
Residuals 14 226975   16213                  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Answer 1

summary在列表中生成了四个表，所以使用lapply就像这样：

mod <- aov(Yield ~ Block + A * B + Error(A * Block), data=aovdat)
lapply(summary(mod), xtable)

将生成一个LaTeX表列表。

带有for的{{1}}循环将仅生成LaTeX代码，表格名称为LaTeX标题。

print

如果您希望将所有表都作为一个表，那么您需要做一些工作，因为for (i in 1:4) print(xtable(summary(mod)[[i]], names(summary(mod)[i])))不喜欢具有不同列数的表。 plyr包中的rbind将完成这项工作

rbind.fill

但我们必须手动将rownames作为数据添加回来，因为a）library(plyr) aovtab <- do.call(rbind.fill, lapply(1:length(summary(mod)), function(x) unclass(summary(mod)[[x]])[[1]]))删除它们，b）它们仍然是重复的，rbind.fill不喜欢。

xtable

Answer 2

在玩了一会儿之后，以下似乎以LaTeX形式产生结果。

xtable(as.data.frame(unlist(coef(anova.mod))))

这给出了块的LaTeX表及其效果，如下所示。

\begin{table}[ht]
\begin{center}
\begin{tabular}{rr}
  \hline
 & unlist(coef(anova.mod)) \\ 
  \hline
(Intercept).(Intercept) & 291.59 \\ 
  A.A1 & 20.81 \\ 
  Block.Block2 & -7.50 \\ 
  Block.Block3 & 5.50 \\ 
  Block.Block4 & -14.25 \\ 
  Block.Block5 & -0.50 \\ 
  Block.Block6 & 11.25 \\ 
  Block.Block7 & -11.75 \\ 
  Block.Block8 & 24.00 \\ 
  A:Block.B1 & 144.25 \\ 
  Within.B1 & 139.00 \\ 
  Within.A1:B1 & 4.33 \\ 
   \hline
\end{tabular}
\end{center}
\end{table}

这是否足以满足您的目的？

通常，xtable更喜欢矩阵或数据框输入，而R经常返回列表，所以我通常做的是提取我感兴趣的部分然后转换为矩阵或数据帧。

为丑陋的TeX道歉，这里没有对它的支持。

LaTeX表通过xtable

2 个答案: