在R中模拟时间序列随机变量？

Question

我必须承认我对R来说是全新的。因此，我的问题可能非常简单。

对于作业，我需要模拟随机游走系列。初始位置是实线上的固定点a。然后以长度X1进行第一步，并将随机游走过程的当前位置改变为S（1）= a + X1。这个过程一直持续到n = 1000。

在生成随机变量X~N（1,10 ^ 2）1000次之后，我需要报告S（n）的值并显示该随机游走系列的图。我还需要报告系列的均值和方差。

这是我到目前为止所做的：

set.seed(1234)
x<-rnorm(1000,mean=1,sd=10)
a<--2

如何生成一个S，以便每n到1000可以得到一个S（n）？ （如果这是一个非常基本的问题，我很抱歉。）

Answer 1

x的每个元素都是一个步骤，因此，行进的总距离，即S（N）的值是元素的总和。此外，每个S(p)行进的距离等于随机行走向量x到索引p的总和。

使用plot和cumsum，您可以生成显示S（n）随时间变化的函数的图表：

不要太粗鲁，但我认为这些基本问题通常是你的主管或老师最能帮助你的。

Answer 2

M=1000; x = rnorm(M); plot(cumsum(x), type='l')