1

我有如下所示的数据框：

set.seed(123)
df <- data.frame(factor1 = rep(c("A", "B"),50),
                 factor2 = rep(c("X","X", "Y", "Y"),25),
                 value = rnorm(100))

我想计算 factor1:factor2 对的一些汇总值。我使用以下方法计算了平均值和标准差：

summary <- as.matrix(cast(df, factor1~factor2, mean))
summary.sd <- as.matrix(cast(df, factor1~factor2, sd))
summary.table <- t(rbind(summary, summary.sd))
colnames(summary.table) <- c("A.mean", "B.mean", "A.sd", "B.sd")

但是我想将比较 A 与 B 的 t.test 中的 p 值添加到 summary.table 中。到目前为止，我已经这样做了，但这不仅没有写入 summary.table，而且我无法获得名称与它一起打印的 factor2 变量：

for (measurement in levels(df$factor2)) print(t.test(value~factor1, data=subset(df, factor2==measurement)))

我认为必须有一些简单的方法来做到这一点，或者也许我不知道的一个包会使这更简单。

4

1 回答 1

3

我会这样做：

一、getmean和sdsummaries using ddplyfrom plyrusingsummarise

require(plyr)
require(reshape2)
o1 <- ddply(df, .(factor1, factor2), summarise, mean = mean(value), sd=sd(value))

#   factor1 factor2       mean        sd
# 1       A       X 0.03746854 0.8730525
# 2       A       Y 0.18352432 0.7635439
# 3       B       X 0.10317706 1.0494930
# 4       B       Y 0.03745372 0.9876173

然后，p-values为t-testNULL和级别mean(A) = mean(B)获取：XYfactor2

o2 <- ddply(df, .(factor2), summarise, pval=t.test(value ~ factor1)$p.value)

#   factor2      pval
# 1       X 0.8108754
# 2       Y 0.5614256

然后，使用reshape2'smelt并dcast转换o1为所需的格式。

o1.mc <- dcast(melt(o1, c("factor1", "factor2")), factor2 ~ variable + factor1)

#   factor2     mean_A     mean_B      sd_A      sd_B
# 1       X 0.03746854 0.10317706 0.8730525 1.0494930
# 2       Y 0.18352432 0.03745372 0.7635439 0.9876173

现在，将其与o2：

merge(o1.mc, o2)

#   factor2     mean_A     mean_B      sd_A      sd_B      pval
# 1       X 0.03746854 0.10317706 0.8730525 1.0494930 0.8108754
# 2       Y 0.18352432 0.03745372 0.7635439 0.9876173 0.5614256
于 2013-02-10T15:21:53.510 回答