【发布时间】:2013-09-13 09:42:01
【问题描述】:
我有一个包含半小时值的大型气象数据集。我正在使用 ddply 函数来创建每日值(平均值、总和等),但是我在获取最小值和最大值时遇到了问题!对于某些值,例如气温我想提取每日平均值、最小值和最大值,但是返回的数据框给了我与平均值相同的最小值和最大值。有人可以看看我写的命令,看看他们是否发现任何错误?
HFTda <- ddply (HFT, .(Year, DOY),summarise,
Temp = mean (Temp, na.rm = F),
Tmax = max (Temp, na.rm = F),
Tmin = min (Temp, na.rm = F),
RHmax = max (RH.x, na.rm = F),
RHmin = min (RH.x, na.rm = F),
VPD = mean (VPD.x, na.rm = F),
Grad = mean (GradInc2, na.rm = F),
Gex = mean (Gex, na.rm = F),
Prec = sum (Prec, na.rm = F),
H = mean (HFill, na.rm = F),
LE = mean (LEFill, na.rm = F),
NEE = mean (co2_fluxFill, na.rm = F),
ET = mean (h2o_fluxFill, na.rm = F),
SWout = mean (YB5214h, na.rm = F))
【问题讨论】:
-
注:每日拆分基于年和日数(DOY)
-
请给出一些可重现的例子。我的建议是将您的调用汇总到一个匿名函数中,并在其中放置一个浏览器调用,以便您可以检查所有值
-
我第二个@KarlForner 将
dput( head( HFT ) )的输出添加到代码块中的问题中。