我真的很接近,但我错过了最后一步,你能告诉我我需要改变下面吗?
我的数据集(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天。我想要每月的平均值,而不是总数)。
我相信这是一个简单的解决办法,但我不知道如何做到这一点?任何建议都是赞赏欢呼桑德拉
1条答案
按热度按时间tp5buhyn1#
假设您有所有年份(*2006年至2016年 *)的数据,这将给予所需的输出:
这是一个快速的解决办法。只要除以年数,你就得到了平均值。