在R

Question

我想了解如何在R中从头开始计算自相关函数。如果cor(x=y, y=lag(x=y, k=2))是ACF对象，如何使用y获取ts？

我已经尝试了use参数[use = "complete.obs" # c("everything", "all.obs", "complete.obs", "na.or.complete", "pairwise.complete.obs")]的所有选择。

set.seed(1)
y <-ts(data = rnorm(20), start = c(2010, 1), frequency = 4)
y

# Qtr1        Qtr2        Qtr3        Qtr4
# 2010  0.91897737  0.78213630  0.07456498 -1.98935170
# 2011  0.61982575 -0.05612874 -0.15579551 -1.47075238
# 2012 -0.47815006  0.41794156  1.35867955 -0.10278773
# 2013  0.38767161 -0.05380504 -1.37705956 -0.41499456
# 2014 -0.39428995 -0.05931340  1.10002537  0.76317575

acf(x=y, plot=FALSE)
# Autocorrelations of series ‘y’, by lag
# 
# 0.00   0.25   0.50   0.75   1.00   1.25   1.50   1.75   2.00   2.25   2.50   2.75   3.00   3.25 
# 1.000 -0.122 -0.185 -0.049  0.147 -0.283 -0.255  0.212  0.097 -0.120 -0.181  0.286 -0.063  0.094 

cor(
      x   = y
    , y   = lag(x=y, k=2)
    , use = "complete.obs" # c("everything", "all.obs", "complete.obs", "na.or.complete", "pairwise.complete.obs")
    )

# [1] 1

Answer 1

看起来你有两个问题。

第一个问题：cor(x, lag(x, k = k))将始终为1，因为这两个系列不是在时间上对齐的。 在使用ts.union

之前，您需要使用cbind.ts或cor

X <- ts.union(yt = y, yt2 = lag(x = y, k =  2))
head(X)
##            yt      yt2
## [1,]       NA -0.62645
## [2,]       NA  0.18364
## [3,] -0.62645 -0.83563
## [4,]  0.18364  1.59528
## [5,] -0.83563  0.32951
## [6,]  1.59528 -0.82047
tail(X)
##              yt      yt2
## [17,]  1.124931 -0.01619
## [18,] -0.044934  0.94384
## [19,] -0.016190  0.82122
## [20,]  0.943836  0.59390
## [21,]  0.821221       NA
## [22,]  0.593901       NA

问题在于没有时间索引，原始数据是相同的（具有时间偏移）。你可以自己测试一下

x1 <- as.vector(y)
x2 <- as.vector(lag(y, k = 2))
all.equal(x1, x2)
## [1] TRUE

因此，如果要计算时间序列与其滞后之间的相关系数，可以使用X（使用ts.union创建）

cor(X[, 1], X[, 2], use = "complete.obs")
## [1] -0.19018

结果仍与acf(y, plot = FALSE)$acf[3]

不同

acf(y, plot = FALSE)$acf[3]
## [1] -0.18521

这导致我们第二个原因为什么您无法使用cor来计算acf：

Acf的数学定义假设至少一个二阶平稳性（每个滞后的均值和方差相等）：

但是如果你使用cor的标准实现，对于每个系列，将计算不同的均值和方差值（滞后的值将与原始系列不同）。

c0 <- var(y)
m <- mean(y)
n <- length(y)
ct <- sum((X[, 1] - m) * (X[, 2] - m), na.rm = TRUE) / (n - 1)
(rt <- ct / c0)
## [1] -0.18521

all.equal(rt, acf(y, plot = FALSE)$acf[3])
## [1] TRUE

Answer 2

根据@khashaa提供的链接得到正确答案。感谢

acf(x=y, plot=FALSE)
# Autocorrelations of series ‘y’, by lag
# 
# 0.00   0.25   0.50   0.75   1.00   1.25   1.50   1.75   2.00   2.25   2.50   2.75   3.00   3.25 
# 1.000 -0.122 -0.185 -0.049  0.147 -0.283 -0.255  0.212  0.097 -0.120 -0.181  0.286 -0.063  0.094 

sum((y-mean(y))*(lag(x=y, k=0)-mean(y)))/sum((y-mean(y))^2)
[1] 1
sum((y-mean(y))*(lag(x=y, k=1)-mean(y)))/sum((y-mean(y))^2)
[1] -0.1222859
sum((y-mean(y))*(lag(x=y, k=2)-mean(y)))/sum((y-mean(y))^2)
[1] -0.1852114
sum((y-mean(y))*(lag(x=y, k=3)-mean(y)))/sum((y-mean(y))^2)
[1] -0.04940401