r - 多子集延迟

Question

我相信 ddply 是我完成任务所需的工具，而且我在获得正确结果时遇到了一些困难。我已经阅读了几个小时关于 ddply 的文章，并尝试了不同的代码，但我自己没有进一步了解。
这是一个示例数据框

station <- c(rep("muc",13), rep("nbw", 17))
year <- c(rep(1994,4),rep(1995,4),rep(1996,5),rep(1994,5), rep(1995,4), rep(1996,4), rep(1997, 4))
depth <- c(rep(c("HUM","31-60","61-90","91-220"),2), rep(c("HUM","0-30", "31-60","61-90","91-220"),2),rep(c("HUM","0-30", "31-60","91-220"),1),rep(c("HUM","0-30", "31-60","61-90"),2))
doc <- c(80, 10, 3, 2,70, 15, 5, 5,70, 20, 5, 5, 2, 40, 10, 3, 2, 1,50, 15, 5, 2, 45, 20, 2, 1,35, 8, 2, 1)

df <-data.frame(station,year,depth,doc)
df

深度是指土壤深度(HUM=腐殖质层)，doc 是测量的土壤深度溶解有机碳 (doc)。请注意，并非每年都有针对文档的测量，并且缺少一些深度类。这很烦人，但经常出现在我的数据集中。
使用 ddply 我想向该数据框中添加一列，以便对于每个深度，返回上述土壤层的文档，并且应该给出 HUM NA，因为腐殖质层顶部没有任何东西。
举个例子:

depth   doc  doc_m1
HUM     80   NA
31-60   10   80
61-90   3    10
91-220  2    3

在数据框中，这当然应该针对每年和每个深度进行计算。我想避免 which 和 for 循环，似乎 ddply 适合于此，但是我没有运气获得滞后命令来使用 ddply。
这是我得到的代码(显然不是很远):

doc <- ddply(df, .(year), transform,
      doc_m1 = ????)

有人有建议吗？
提前致谢!

Answer 1

如果您的数据集中的深度已经按正确顺序排列(如您的示例中所示)，您可以这样做:

doc2 <- ddply(df, .(station, year), transform,
      doc_m1 = c(NA, doc[-length(doc)]))

注意我也在车站分组。这给出:

> head(doc2, 10)
   station year  depth doc doc_m1
1      muc 1994    HUM  80     NA
2      muc 1994  31-60  10     80
3      muc 1994  61-90   3     10
4      muc 1994 91-220   2      3
5      muc 1995    HUM  70     NA
6      muc 1995  31-60  15     70
7      muc 1995  61-90   5     15
8      muc 1995 91-220   5      5
9      muc 1996    HUM  70     NA
10     muc 1996   0-30  20     70

如果它们尚未按深度排序，则将深度作为具有正确顺序的级别的因子，然后据此进行排序。那么这种方法应该有效。