我使用动物园聚合函数来获取每日数据的月平均值:
monthlyMeanTemp <- aggregate(merged.precip.flow.and.T[,'E'], as.yearmon, mean, na.rm = TRUE) # ‘E’ is the column of temperature
以下是结果的头尾:
Jan 1979 Feb 1979 Mar 1979 Apr 1979 May 1979 Jun 1979
-14.05354839 -11.83078929 -7.32150645 -0.03214333 6.16986774 14.00944000
...
Apr 1997 May 1997 Jun 1997 Jul 1997 Aug 1997 Sep 1997
1.438547 7.421910 12.764450 15.086206 17.376026 10.125013`
是否有可能按月获得平均值(即所有1月值的平均值,所有2月值的平均值等),而不使用NA填充缺失月份,形成焦虑12矩阵(其中n是年数),然后使用colMeans函数?
答案 0 :(得分:1)
...刚刚找到答案:来自hydroTSM包:monthlyfunction(merged.precip.flow.and.T[,'E'], FUN=mean, na.rm=TRUE)