Question

我跑了一圈：

N = 10
alpha = 0.05
sigma = 0.01
for (i in 0:N){
Vt10[i] = ((sigma^2)/(alpha^2))*((N-i)+(2/alpha))*exp(-alpha*(N-i))-(1/(2*alpha))*exp(-2*alpha*(N-i))-(3/(2*alpha))
}

但是，Vt10最终只输出10个结果（不包括第一次迭代，i = 0时）。如何创建包含所有11次迭代的向量？

Answer 1

<强>答案：

N = 10
alpha = 0.05
sigma = 0.01
Vt10 = numeric(11)  ## this!
for (i in 0:N){
Vt10[i+1] = ((sigma^2)/(alpha^2))*((N-i)+(2/alpha))*exp(-alpha*(N-i))-(1/(2*alpha))*exp(-2*alpha*(N-i))-(3/(2*alpha))
}

此外，您需要Vt10[i+1]，因为您从0循环到10，而数组索引应该是1到11。

<强>注释：

您的代码看起来像C代码。不仅从0开始索引，而且还为此任务编写循环。

尝试：

N = 10
alpha = 0.05
sigma = 0.01
i = 0:10
Vt10 = ((sigma^2)/(alpha^2))*((N-i)+(2/alpha))*exp(-alpha*(N-i))-(1/(2*alpha))*exp(-2*alpha*(N-i))-(3/(2*alpha))

在这种情况下，不需要预定义矢量。 R知道你正在执行逐元素计算，并会自动为Vt10分配长度为11的向量。

保存向量R中的循环结果

1 个答案: