让我们说a_i是iid标准普通r.v.的i = 1,...,100
我试图想象时间序列z_t =从i = 1到t的总和(从j = 1到(a_i)的i之和)
所以,这相当于z_t = z_(t-1)+从i = 1到a的总和__
这是我到目前为止所做的:
set.seed(42)
a<-c(rnorm(100))
我的想法是创建一个函数z:{1,... 100} - &gt; R将执行这个双和。然后我可以在t = 1到t = 100
之间对函数进行矢量化和绘图更新,这是我最终做的事情
a<-c(rnorm(100))
z<-function(t){
sum(cumsum(a)[1:t])
}
y<-1:100
for(i in 1:100){
y[i]<-z(i)
}