如何制作一个带有字母的anova表，表示groupwise在R->latex上的显著差异？

Question

在R中有没有一个包可以生成这样的表：

基本上，给定一个因素和变量的数据集，产生一个格式良好的表格，其中包含成对的t检验和字母，表明差异的显著性？

虚拟数据：

var = c(rnorm(100,mean=1,sd=2),rnorm(30,mean=-1,sd=2),rnorm(50,mean=0,sd=4))
factor = as.factor(c(rep(1,100),rep(2,30),rep(3,50))

我如何获取虚拟数据并将其转换为示例表的(一列)？

Answer 1

免责声明

不是一个确切的答案，但这可能是一个好的开始！使用tables包。

library(Hmisc)
library(tables)
stderr <- function(x) sd(x)/sqrt(length(x))
latex(
  object = tabular((Species+1) ~ All(iris)* PlusMinus(mean, stderr, digits=1), data=iris)
 , title = "Test"
 , file=""
 , size = "small"
 , cdot = 3
 , here = TRUE
 , booktabs=TRUE
 , center="centering"
 )

输出

​

​

编辑的

您可以使用multcompView包获取信件。

代码

library(multcompView)

Sepal.Length.fm <- aov(Sepal.Length~Species, data=iris)
Sepal.Length.Letters <- data.frame("Letters"=multcompLetters(extract_p(TukeyHSD(Sepal.Length.fm)$"Species"))$"Letters")

Sepal.Width.fm <- aov(Sepal.Width~Species, data=iris)
Sepal.Width.Letters <- data.frame("Letters"=multcompLetters(extract_p(TukeyHSD(Sepal.Width.fm)$"Species"))$"Letters")

Petal.Length.fm <- aov(Petal.Length~Species, data=iris)
Petal.Length.Letters <- data.frame("Letters"=multcompLetters(extract_p(TukeyHSD(Petal.Length.fm)$"Species"))$"Letters")

Petal.Width.fm <- aov(Petal.Width~Species, data=iris)
Petal.Width.Letters <- data.frame("Letters"=multcompLetters(extract_p(TukeyHSD(Petal.Width.fm)$"Species"))$"Letters")

Letters <- cbind(Sepal.Length.Letters, Sepal.Width.Letters, Petal.Length.Letters, Petal.Width.Letters)

输出

           Letters Letters Letters Letters
versicolor       a       a       a       a
virginica        b       b       b       b
setosa           c       c       c       c

Answer 2

您可以使用aggregate函数来完成此操作。

df = as.data.frame(cbind(var,factor))
aggregate(df$var, by=list(df$factor), t.test)

  Group.1         x
1       1  4.939821
2       2 -2.128924
3       3 -1.431482

如果希望多个测试显示为新列，可以创建一个新函数作为最后一个参数。