# 毎正時のタイムスタンプを作成
px <- seq(as.POSIXct('2023-01-01 01:00:00'),
          as.POSIXct('2023-01-01 12:00:00'), by = 60 * 60)

n <- length(px)

# ホワイトノイズ
set.seed(5963)
x <- rnorm(n, mean = 0, sd = 1)

# 時系列を扱うts型オブジェクトを作成
y <- ts(x) # 原系列

y.lag1 <- lag(y, k = -1) # 1次ラグ系列
y.lead1 <- lag(y, k = 1) # 1次リード系列

ts4 <- cbind(y, y.lag1, y.lead1)[1:n, ]

d <- as.data.frame(ts4)
d

##              y      y.lag1     y.lead1
## 1           NA          NA  1.21767531
## 2   1.21767531          NA -0.54813309
## 3  -0.54813309  1.21767531  0.16230527
## 4   0.16230527 -0.54813309  0.02054537
## 5   0.02054537  0.16230527  0.25239422
## 6   0.25239422  0.02054537  0.40647836
## 7   0.40647836  0.25239422 -0.22756814
## 8  -0.22756814  0.40647836 -0.92278941
## 9  -0.92278941 -0.22756814  0.40292174
## 10  0.40292174 -0.92278941 -1.55653902
## 11 -1.55653902  0.40292174 -0.35289239
## 12 -0.35289239 -1.55653902 -0.71134333

matplot(x = px, y = d, type = 'o', lty = 1:3, col = 1:3, pch = 1:3, 
        main = '原系列とラグ系列・リード系列', xlab = '時刻', ylab = '値')

grid()

legend('topright', lty = 1:3, col = 1:3, pch = 1:3,
       legend = c('原系列', '1次ラグ系列', '1次リード系列'))