平均县（fips）降水量并使用 dplyr 删除 stationID

Question

示例数据集：https://www.dropbox.com/s/f3l2uub1cttwmf2/test.csv?dl=0

我需要做的是逐年获取每个 fips 代码的平均 APPT，然后将其减少到每个 fips 代码仅 1 次观察。这应该给我每个 fips 代码的平均 APPT（每年平均降水量）。 gridNumber 变量是位于县内的站点 ID，因此通过对这些进行平均，我应该得到 fips 代码的 1 个平均降水量。

我正在尝试使用 dplyr 执行此操作，但任何事情都可以。这是一个包含美国所有县的相当大数据集的一个小样本。

这是我尝试过的：

df3 <- test %>%
     group_by(fips, Year) %>%
     summarize(ppt = mean(APPT))

感谢您的帮助

Answer 1

我假设您也必须加载plyr。请查看dplyr issue。一种方法是在新的 R 控制台上执行此操作或分离 plyr 甚至 dplyr::summarize 会有所帮助。 plyr 也有 summarize 函数，因此该函数可能正在使用 plyr::summarize。

df3 <- test %>%
            group_by(fips, Year) %>%
            dplyr::summarize(ppt = mean(APPT))

话虽如此，我建议不要将summarise 用于多个组，因为

 test %>%
      group_by(fips, Year) %>% 
      summarize(ppt=mean(APPT))
 #Source: local data frame [34 x 3]
 #Groups: fips #omitted the last 'Year`

虽然结果是一样的，但最近的一篇文章（没有链接）表明此分组函数与其他分组函数之间的结果存在差异。如果速度是一个问题，你可以使用data.table

library(data.table)
setDT(test)[, list(ppt=mean(APPT)), by=list(fips, Year)]

或使用base R

aggregate(APPT~fips+Year, test, mean)