8

我使用以下代码来总结我的数据，按 Compound、Replicate 和 Mass 分组。

summaryDataFrame <- ddply(reviewDataFrame, .(Compound, Replicate, Mass), 
  .fun = calculate_T60_Over_T0_Ratio)

一个不幸的副作用是生成的数据帧按这些字段排序。我想这样做并保持 Compound、Replicate 和 Mass 的顺序与原始数据框中的顺序相同。有任何想法吗？我尝试将顺序整数的“排序”列添加到原始数据中，但是我当然不能将其包含在 .variables 中，因为我不想对其进行“分组”，因此它不会在摘要数据帧。

谢谢您的帮助。

4

2 回答 2

11

这在plyr邮件列表中出现了一段时间（由@kohske 提出），这是 Peter Meil​​strup 为有限情况提供的解决方案：

#Peter's version used a function gensym to
# create the col name, but I couldn't track down
# what package it was in.
keeping.order <- function(data, fn, ...) { 
  col <- ".sortColumn"
  data[,col] <- 1:nrow(data) 
  out <- fn(data, ...) 
  if (!col %in% colnames(out)) stop("Ordering column not preserved by function") 
  out <- out[order(out[,col]),] 
  out[,col] <- NULL 
  out 
} 

#Some sample data 
d <- structure(list(g = c(2L, 2L, 1L, 1L, 2L, 2L), v = c(-1.90127112738315, 
-1.20862680183042, -1.13913266070505, 0.14899803094742, -0.69427656843677, 
0.872558638137971)), .Names = c("g", "v"), row.names = c(NA, 
-6L), class = "data.frame") 

#This one resorts
ddply(d, .(g), mutate, v=scale(v)) #does not preserve order of d 

#This one does not
keeping.order(d, ddply, .(g), mutate, v=scale(v)) #preserves order of d

请阅读有关 Hadley 注释的线程，了解为什么此功能可能不够通用，无法转入ddply，特别是因为它可能适用于您的情况，因为您可能每件返回的行数较少。

编辑以包含针对更一般情况的策略

如果ddply正在输出按您不喜欢的顺序排序的内容，则基本上有两个选择：使用有序因子预先指定拆分变量的所需排序，或者在事后手动对输出进行排序。

例如，考虑以下数据：

d <- data.frame(x1 = rep(letters[1:3],each = 5), 
                x2 = rep(letters[4:6],5),
                x3 = 1:15,stringsAsFactors = FALSE)

暂时使用字符串。ddply将对输出进行排序，在这种情况下将需要默认的词法排序：

> ddply(d,.(x1,x2),summarise, val = sum(x3))
  x1 x2 val
1  a  d   5
2  a  e   7
3  a  f   3
4  b  d  17
5  b  e   8
6  b  f  15
7  c  d  13
8  c  e  25
9  c  f  27


> ddply(d[sample(1:15,15),],.(x1,x2),summarise, val = sum(x3))
  x1 x2 val
1  a  d   5
2  a  e   7
3  a  f   3
4  b  d  17
5  b  e   8
6  b  f  15
7  c  d  13
8  c  e  25
9  c  f  27

如果生成的数据框没有以“正确”的顺序结束，那可能是因为您确实希望其中一些变量成为有序因子。假设我们真的想要x1并x2像这样订购：

d$x1 <- factor(d$x1, levels = c('b','a','c'),ordered = TRUE)
d$x2 <- factor(d$x2, levels = c('d','f','e'), ordered = TRUE)

现在，当我们使用 时ddply，结果排序将如我们所愿：

> ddply(d,.(x1,x2),summarise, val = sum(x3))
  x1 x2 val
1  b  d  17
2  b  f  15
3  b  e   8
4  a  d   5
5  a  f   3
6  a  e   7
7  c  d  13
8  c  f  27
9  c  e  25

这里故事的寓意是，如果ddply以您不想要的顺序输出某些东西，这是一个好兆头，表明您应该对要拆分的变量使用有序因子。

于 2011-08-29T20:30:23.057 回答
1

我最终在原始数据框中添加了一个“索引”列。它由两列pasted组成sep="_"。然后我制作了另一个数据框，仅unique由“索引”列的成员和一个计数器组成1:length(df)。我ddply()对返回排序数据框的数据进行了处理。然后为了让事情恢复到原来的顺序，我做merge()了结果数据框和索引数据框（确保列被命名为相同的东西，这样更容易）。最后，我做了order并删除了无关的列。

不是一个优雅的解决方案，而是一个有效的解决方案。

感谢您的协助。它让我朝着正确的方向思考。

于 2011-08-30T15:57:54.490 回答