[英]save results of loop in a vector R
我运行一个循环:
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个迭代的向量?
回答:
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的向量。
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.