r语言 - 二进制运算符的非数字参数，AR(1) 模型

Question

我有一个练习要做，我必须运行以下 AR（1） 模型：

xi =c+φxi−1+ηi (i=1,...,T)

我知道 ni ~ N（

0,1） ; x0 ~ N（c/（1-φ），1/（1-φˆ2））; c= 2 ; φ = 0.6我正在尝试做一个 for 循环。我的代码如下：

n <- rnorm(T, 0, 1)
c <- 2
phi <- 0.6
x_0 <- rnorm(1,c/(1-phi), 1/(1-phi**2))
v <- vector("numeric", 0)
#for (i in 2:T){
  name <- paste("x", i, sep="_")
  v <- c(v,name)
  v[1] <- c + phi*x_0 + n[1]
  v[i] <- c + phi*v[i-1] + n[i]
}

但是，我不断收到此错误：

Error in phi * v[i - 1] : non-numeric argument to binary operator

我了解此错误是什么，但我找不到任何解决方案来解决它。有人可以开导我吗？如何将数值分配给名称向量？

谢谢！

Answer 1

您将v定义为数字向量，但随后v <- c(v, name) v转换为字符向量，因为name是字符。这就是导致错误的原因。

如果我没记错的话，您的目的是为数字向量中的值分配名称。没关系，你只需要一种不同的方法。

n <- rnorm(t)
c <- 2
phi <- 0.6
x_0 <- rnorm(1, c/(1-phi), 1/(1-phi^2))
v <- c + phi*x_0 + n[1]
for (i in 2:t) {
    v[i] <- c + phi*v[i-1] + n[i]
}
names(v) <- paste("x", 1:t, sep="_")

R 中的向量没有静态大小;它们会根据需要动态调整大小。因此，即使我们使用标量值初始化v，它也会增长以适应循环中的每个新值。

最后一步是给v一个名字列表。这可以使用names(v) <- .现在看看v——它有名字！

顺便说一句

，由于 T 是 R 中TRUE的同义词，因此最好不要使用 T 作为变量名。因此，我在这里使用了t。

Answer 2

我想你似乎需要以下内容。它将生成 11 个元素，包括初始 x 值。您可以稍后将其排除。

set.seed(1237)
t <- 10
n <- rnorm(t, 0, 1)
c <- 2
phi <- 0.6
x0 <- rnorm(1, c/(1-phi), 1/(1-phi**2))
v <- c(x0, rep(0, t))
for(i in 2:length(v)) {
  v[i] <- c + phi * v[i-1] + n[i-1]
}
v
[1] 4.967833 4.535847 2.748292 2.792992 5.389548 6.173001 4.526824 3.790483 4.307981 5.442913 4.958193