Question

我目前被一些EddyData的后期处理所困扰。在一个示例（https://github.com/bgctw/REddyProc/blob/master/vignettes/aggUncertainty.md之后，我想到了group_by的意外结果，该结果是可重复的，但我不明白为什么。

Group_by仅返回一行，而聚合则给出预期的结果。

这是一个最小的示例：

library(tidyverse)

#create example data frame
date.time <- seq(from=as.POSIXct("2015-01-01 00:30:00"), to=as.POSIXct("2015-01-03 00:30:00"),by="30 mins")
nee <- runif(length(date.time),-200,200)
df <- data.frame(date.time, nee)

#calculate day of the year
df <- df %>% mutate(
  date.time = df$date.time
  , DoY = as.POSIXlt(date.time - 15*60)$yday # midnight belongs to the previous
)

#trying to summarise nee for each day
aggDay <- df %>% group_by(DoY) %>% summarise(nee=sum(nee))
aggDay
       nee
1 322.1195

aggDay仅返回一行，而在这种情况下aggregate将起作用

aggregate(df$nee, by=list(df$DoY), sum)
  Group.1          x
1       0  -25.15698
2       1  448.13960
3       2 -100.86310

不幸的是，原始代码涉及到一些进一步的计算，这就是我想呆在group_by上的原因。

#original code, not reproducible here
aggDay <- df %>% group_by(DoY) %>% 
  summarise(
    DateTime = first(DateTime)
    , nRec = sum( NEE_uStar_fqc == 0, na.rm = TRUE)
    , nEff = computeEffectiveNumObs(
      resid, effAcf = !!autoCorr, na.rm = TRUE)
    , NEE = mean(NEE_uStar_f, na.rm = TRUE)
    , sdNEE = if (nEff <= 1) NA_real_ else sqrt(
      mean(NEE_uStar_fsd^2, na.rm = TRUE) / (nEff - 1)) 
    , sdNEEuncorr = if (nRec == 0) NA_real_ else sqrt(
      mean(NEE_uStar_fsd^2, na.rm = TRUE) / (nRec - 1))
  )

Answer 1

我重新启动了RStudio，现在可以使用了。不要问我另一个加载的软件包一定存在问题。

Group_by仅返回一行，而合计则返回预期结果

1 个答案: