7

我正在使用 dplyr，我想知道是否可以在一行中计算组之间的差异。如下面的小示例所示，任务是计算 A 组和 B 组标准化“分”变量之间的差异。

library(dplyr)
# creating a small data.frame
GROUP <- rep(c("A","B"),each=10)
NUMBE <- rnorm(20,50,10)
datf <- data.frame(GROUP,NUMBE)

datf2 <- datf %.% group_by(GROUP) %.% mutate(cent = (NUMBE - mean(NUMBE))/sd(NUMBE))

gA <- datf2 %.% ungroup() %.% filter(GROUP == "A") %.% select(cent)
gB <- datf2 %.% ungroup() %.% filter(GROUP == "B") %.% select(cent)

gA - gB

通过创建不同的对象当然没有问题 - 但是是否有更“内置”的方式来执行此任务？更像这样的东西在下面不起作用的幻想代码？

datf2 %.% summarize(filter(GROUP == "A",select(cent)) - filter(GROUP == "B",select(cent)))

谢谢！

4

2 回答 2

6

鉴于我们每组有 10 个，添加一个索引 1:10, 1:10 并对其进行总结：

> datf2$entry=c(1:10,1:10)
> datf2 %.% ungroup() %.% group_by(entry) %.% summarize(d=cent[1]-cent[2])
Source: local data frame [10 x 2]

   entry          d
1      1 -0.8272879
2      2 -0.9159827
3      3 -0.5064762
4      4  0.4211639
5      5  1.3681720
6      6  3.3430289
7      7  1.0086822
8      8 -0.6163907
9      9 -0.7325220
10    10 -2.5423875

相比：

> gA - gB
         cent
1  -0.8272879
2  -0.9159827
3  -0.5064762
4   0.4211639
5   1.3681720
6   3.3430289
7   1.0086822
8  -0.6163907
9  -0.7325220
10 -2.5423875

有没有办法将entry字段注入数据或dplyr调用？我不确定，它似乎依赖于对数据了解太多的函数......

于 2014-03-23T12:42:33.773 回答
4

感谢您的启发。我进一步开发了这个解决方案：

mutate(datf2,diffence = filter(datf2, GROUP == "A")$cent - filter(datf2, GROUP == "B")$cent)

这会将结果添加为 data.frame 中的列。

于 2014-03-24T18:45:37.917 回答