R语言 每组平均计数

gpnt7bae  于 2023-04-03  发布在  其他
关注(0)|答案(1)|浏览(117)

我真的很接近,但我错过了最后一步,你能告诉我我需要改变下面吗?
我的数据集(mydata)看起来像这样(前几个条目丢失):

date_daily station_date air_min air_avg air_max rain_mm temp_davg_c temp_dmax_c temp_dmin_c       date months years
1 2006/01/01         <NA>    <NA>    <NA>    <NA>      NA          NA          NA          NA 2006-01-01     01  2006
2 2006/01/02         <NA>    <NA>    <NA>    <NA>      NA          NA          NA          NA 2006-01-02     01  2006
3 2006/01/03         <NA>    <NA>    <NA>    <NA>      NA          NA          NA          NA 2006-01-03     01  2006
4 2006/01/04         <NA>    <NA>    <NA>    <NA>      NA          NA          NA          NA 2006-01-04     01  2006
5 2006/01/05         <NA>    <NA>    <NA>    <NA>      NA          NA          NA          NA 2006-01-05     01  2006
6 2006/01/06         <NA>    <NA>    <NA>    <NA>      NA          NA          NA          NA 2006-01-06     01  2006

每天的数据范围从2006年到2016年,所以我感兴趣的是平均月降雨量以及每月降雨量〉2毫米的平均天数。
所以我尝试了以下方法:

library(dplyr)
stats <- mydata %>% group_by(months) %>%
  summarise(days_w_rain = length(rain_mm[!is.na(rain_mm)]),
            avg_rain = mean(rain_mm, na.rm=TRUE),
            days_2mm = sum(rain_mm > 2, na.rm=TRUE),
            p2mm = days_2mm / days_w_rain*100)

这给了我以下结果:

# A tibble: 12 x 5
   months days_w_rain   avg_rain days_2mm      p2mm
    <chr>       <int>      <dbl>    <int>     <dbl>
 1     01         304 11.8078520      148 48.684211
 2     02         310  7.4946903      125 40.322581
 3     03         334  9.2680090      156 46.706587
 4     04         327  3.7469113       92 28.134557
 5     05         310  1.6073032       39 12.580645
 6     06         300  1.1322100       39 13.000000
 7     07         310  0.6212226       25  8.064516
 8     08         333  1.6641862       47 14.114114
 9     09         308  4.2796169       71 23.051948
10     10         310 10.7842742      155 50.000000
11     11         300 18.4062900      189 63.000000
12     12         287 20.6697038      155 54.006969

这是每个月平均降雨量的正确数字,但它给了我所有一月降雨的总天数(即304天),而不是一月份的平均值,与降雨量〉2毫米的天数相同(即2006年至2016年1月所有月份的148天。我想要每月的平均值,而不是总数)。
我相信这是一个简单的解决办法,但我不知道如何做到这一点?任何建议都是赞赏欢呼桑德拉

tp5buhyn

tp5buhyn1#

假设您有所有年份(*2006年至2016年 *)的数据,这将给予所需的输出:

stats <- mydata %>% group_by(months) %>%
          summarise(days_w_rain = length(rain_mm[!is.na(rain_mm)]/length(unique(years))),
                    avg_rain = mean(rain_mm, na.rm=TRUE),
                    days_2mm = sum(rain_mm > 2, na.rm=TRUE)/length(unique(years)),
                    p2mm = days_2mm / days_w_rain*100)

这是一个快速的解决办法。只要除以年数,你就得到了平均值。

相关问题