如何在R中的10天時間視窗中重新排列時間序列資料？-有解無憂

我目前正在研究從 2002-10-04 到 2012-10-28 的時間序列流資料，我需要做的是使用 10 天時間視窗資料開發回歸模型。更具體地說，我需要使用 2002 年至 2012 年所有年份的 Oct-04 至 Oct 13 來構建回歸模型。然后我需要使用從 2002 年到 2012 年的所有年份的 Oct-05 到 10 月 14 日來構建另一個回歸模型，然后使用所有年份的 Oct-06 到 10 月 15 日來構建下一個模型，如此反復直到最后。

這就是我的資料的樣子。

> head(CFbasin)
        Data      Qkl     Qllt     Qhaw      Qdp     Qlit
1 2002-10-04 19.25546 8.353470 4.502379 2.217209 1.985011
2 2002-10-05 19.56694 8.126935 4.615646 1.622555 1.628219
3 2002-10-06 19.73684 7.560598 4.389111 1.251605 1.492298
4 2002-10-07 18.12278 7.079212 3.992675 1.158159 1.413011
5 2002-10-08 18.12278 6.824360 3.794457 1.070377 1.393189
6 2002-10-09 17.83961 6.739409 3.369705 1.073208 1.353545

> tail(CFbasin)
           Data      Qkl     Qllt      Qhaw      Qdp     Qlit
3673 2012-10-23 24.89051 16.67862  8.608321 2.477724 1.432832
3674 2012-10-24 25.00378 16.48040 11.638224 2.820358 1.393189
3675 2012-10-25 25.37189 16.99011  7.758816 3.001586 1.322397
3676 2012-10-26 26.07982 16.87684  6.484558 2.814695 1.279921
3677 2012-10-27 27.41071 17.07506  4.813864 3.086536 1.228951
3678 2012-10-28 28.88318 17.16001  5.635052 3.114853 1.220456

我只試過一次，這是我的代碼：

CFbasin %>% filter(month(Date) == 10 & day(Date) >= 4 & day(Date) <= 14)

這使我可以獲取 2002 年至 2012 年 10 月 4 日至 10 月 13 日的所有資料，然后進行線性回歸。但我不確定如何讓它在整個資料集上作業，然后進行線性回歸，我正在考慮forloop 和 function rollapply()，但真的不清楚如何安排我的資料集。任何建議和推薦將不勝感激，在此先感謝您！

uj5u.com熱心網友回復：

這是一個使用末尾注釋中的資料和寬度為 5 的示例。

library(zoo)
coefs <- function(x) coef(lm.fit(cbind(1, x[, -1]), x[, 1]))
rollapplyr(CFbasin[, -1], 5, coefs, by.column = FALSE, fill = NA)

添加

這將 NA 用于第一個輸出 4 行，然后是與所有年份的 1:5 行對應的行以形成下一個回歸，所有年份的 2:6 行形成之后的回歸，依此類推。閏年不使用 12 月 31 日。

yday <- as.POSIXlt(CFbasin$Data)$yday
coefs <- function(ix) {
  x <- CFbasin[yday %in% ix, -1]
  if (NROW(x) == 0) NA else coef(lm(Qkl ~., x))
}
rollapplyr(0:364, 5, coefs, fill = NA)

筆記

Lines <- "
        Data      Qkl     Qllt     Qhaw      Qdp     Qlit
1 2002-10-04 19.25546 8.353470 4.502379 2.217209 1.985011
2 2002-10-05 19.56694 8.126935 4.615646 1.622555 1.628219
3 2002-10-06 19.73684 7.560598 4.389111 1.251605 1.492298
4 2002-10-07 18.12278 7.079212 3.992675 1.158159 1.413011
5 2002-10-08 18.12278 6.824360 3.794457 1.070377 1.393189
6 2002-10-09 17.83961 6.739409 3.369705 1.073208 1.353545"
CFbasin <- read.table(text = Lines)

轉載請註明出處，本文鏈接：https://www.uj5u.com/qiye/353849.html

標籤：r 数据框时间序列

上一篇：組資料框和單元格值與中值計算之間的計算比率

下一篇：從資料框中僅獲取唯一的行