水文中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