使用聚合函数的嵌套循环中的错误

Question

我有一个 data.frame，我想在其中计算嵌套在每年和每个站点中的每个图的年度数据的最小值、平均值和最大值以及周期（具有正值的月数）。我正在尝试使用嵌套循环和聚合函数来执行此操作，但一直遇到错误。

这是我迄今为止尝试过的

DT5 <- NULL
for(i in levels(DT4$Site)) {
  for(j in levels(DT4$Year)) {
    tmp <- subset(subset(DT4, Site == i), Year== j)
    min <- aggregate(tmp[,5], list(tmp$Plot), min)
    mean <- aggregate(tmp[,5], list(tmp$Plot), mean)
    max <- aggregate(tmp[,5], list(tmp$Plot), max)
    per <- sum(tmp[,5] > 0)
    tmp <- cbind(rep(i, nrow(mean)), rep(j, nrow(mean)), tmp$Year, mean)
  }
  if(is.null(DT5)){DT5<-tmp} else {DT5<-rbind(DT5,tmp)} 
}

其中 DT4 是包含 Site、Year、Month、Plot 和 WatLev 列的 data.frame

最终，我想实现一个 data.frame，其中包含 Site、Year、Plot、min、mean、max 和 per 为每个地块计算的列，在每个站点中每年。

这是一个示例数据集：

Site <- rep(rep(c("SiteA", "SiteB", "SiteC"), each=144),times=3)
Year <- rep(rep(2001:2003, each=48),times=3)
Month <- rep(rep(rep(c("Jan", "Feb", "Mar", "Apr", "May", "Jun", "Jul", "Aug", "Sep", "Oct", "Nov", "Dec"), each=4), times=3),3)
Plot <- rep(rep(c("A", "B", "C", "D"), times=36),3)
WatLev <- runif(1296, -50, 5)
DT4 <- cbind(Site, Year, Month, Plot, WatLev)

Answer 1

好的，试试这个：

library(dplyr)
DT4 %>% as_tibble() %>%
  mutate(WatLev = as.numeric(WatLev)) %>% 
  group_by(Site, Year) %>% 
  summarise(min = min(WatLev, na.rm = TRUE),
            mean = mean(WatLev, na.rm = TRUE),
            max = max(WatLev, na.rm = TRUE),
            per = sum(WatLev > 0))

输出是：

# A tibble: 9 x 6
# Groups:   Site [3]
  Site  Year    min  mean   max   per
  <chr> <chr> <dbl> <dbl> <dbl> <int>
1 SiteA 2001  -49.8 -24.3  4.76    11
2 SiteA 2002  -50.0 -23.6  4.98     9
3 SiteA 2003  -49.2 -20.7  4.76    16
4 SiteB 2001  -49.3 -21.4  4.85    11
5 SiteB 2002  -49.7 -19.0  4.72    19
6 SiteB 2003  -50.0 -22.0  4.89    17
7 SiteC 2001  -49.2 -20.4  4.97    14
8 SiteC 2002  -49.6 -22.2  4.28    13
9 SiteC 2003  -49.0 -23.3  4.85    14