如何将汇总函数应用于时段子集?
例如:
r['T16:00/T17:00']$Value
如何在每天的采样时间内为价值function (x) quantile(x, c(.90))
应用{/ p}。
答案 0 :(得分:3)
在完成时间子集后,您可以使用apply.daily
将函数应用于每天的数据。
rt <- r['T16:00/T17:00','Value']
rd <- apply.daily(rt, function(x) xts(t(quantile(x,0.9)), end(x)))
您可以看到我需要执行一些后空翻,以确保apply.daily
可以处理从您的函数返回的对象。主要是,它必须是一行多列xts对象。