我正在尝试利用 FitAR 包和自回归模型/AR(1)——参见下面的#A——将噪声(例如白噪声/随机噪声,参见下面的#B)与 lynx 计数进行比较。设置它一直令人困惑。我正在提取我遇到的随机噪声示例和来自 FitAR 的 lynx 数据。白噪声模型将有助于确定 lynx 数据中可能具有重要意义的内容。
#猞猁
install.packages(FitAR)
library(FitAR)
library(lattice)
library(leaps)
library(ltsa)
library(bestglm)
help("FitAR-package")
par(mfrow=c(1,2))
lynx <- (log(lynx))
ans <- FitAR((lynx),1)
z4<-Boot.FitAR(ans)
par(mfrow=c(2,1))
TimeSeriesPlot((lynx))
title(main="lynx")
TimeSeriesPlot(z4)
title(main="Simulated AR lynx")
#B 白噪声
install.packages("compositions")
library("compositions")
rnorm(n, mean = 0, sd = 1)
set.seed(100)
x <- NULL
x[1] <- 0
for (i in 2:100) {
x[i] <- x[i-1] + rnorm(1,0,1)
}
ts.plot(x, main = 'Random walk', xlab = '', ylab = '', col='blue', lwd = 2)