水文中r的每小时时间序列分析

Question

我正在尝试将事件与流数据分开。 我有每小时的数据。 我已经运行了代码

dailyMQ <- data.frame(Date=seq(from=as.Date("01.01.2000", format="%d.%m.%Y"), 
to=as.Date("01.01.2004", format="%d.%m.%Y"), by="days"), 
discharge=rbeta(1462,2,20)*100)

用于每日数据。 但我正在尝试每小时数据但遇到错误。 谁能建议我如何为每小时数据编写代码？

谢谢

Answer 1

Date格式不能直接拆分为小时。

您可以使用POSIXct日期时间格式：

HourlyMQ <- data.frame(Date=seq(from=as.POSIXct("01.01.2019", format="%d.%m.%Y"), to=as.POSIXct("11.12.2019", format="%d.%m.%Y"),by="hours"),discharge=rbeta(8257,2,20)) 
HourlyMQ
#>                     Date    discharge
#> 1    2019-01-01 00:00:00 0.2452214482
#> 2    2019-01-01 01:00:00 0.0620291334
#> 3    2019-01-01 02:00:00 0.0608788870
#> 4    2019-01-01 03:00:00 0.0697449808
#> 5    2019-01-01 04:00:00 0.0302780135