計算ACF手動

Question

我對ACF計算下面的公式，我我想知道我使用R中的功能ACF，但手動計算它下面的數據，然後繪製結果使用GGPLOT2

data<-rweibull(168,2,3) 
time<-seq(from=as.POSIXct("2014-01-01 00:00"),to=as.POSIXct("2014-01-07 23:00"),by="hour") 
df<-data.frame(time,data) 

其中k是時滯，μ是時間序列的均值（，Yt = 1,2，... .N）。考慮到時間滯後100小時，ACF在此計算時間序列。在我的最終輸出中，我想要有一個具有100個值的向量。在這個等式中，我不知道如何設置R中方程的Yt + k部分100個時滯。任何幫助，將不勝感激。

現在我有代碼但我不知道如何設置代碼以獲得100滯後結果在一個向量中。

x <- data - mean(data) 
sum(x[1:99] * x[2:100])/sum(x^2)# TIME LAG OF 100 

Answer 1

，如果你想你可以試試這個...

data<-rweibull(168,2,3) 
time<-seq(from=as.POSIXct("2014-01-01 00:00"),to=as.POSIXct("2014-01-07 23:00"),by="hour") 
df<-data.frame(time,data) 
x <- data - mean(data) 

maxLag<-100 
acfResults<- data.frame(k=c(1:maxLag), acf=NA, stringsAsFactors = FALSE) 


for (k in 1: maxLag) { 
    N <-nrow(df)-k 
    acf<- sum(x[1:N]*x[(1+k):(N+k)])/sum(x^2) 
    acfResults[k,1] <-k 
    acfResults[k,2] <-acf 
}