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

Question

6

我必须承认我对R完全不熟悉。因此，我的问题可能非常简单。

对于一项任务，我需要模拟一个随机游走序列。初始位置是实线上的一个固定点a。然后采取长度为X1的第一步，并将随机游走过程的当前位置更改为S(1)=a+X1。这个过程一直持续到n=1000。

生成了1000次的随机变量X~N(1,10^2)之后，我需要报告S(n)值并展示这个随机游走序列的图表。我还需要报告序列的均值和方差。

这就是我目前的进展：

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

我该如何生成一个S，以便我可以为每个小于1000的n拥有一个S(n)?（如果这是一个非常基础的问题，那么我很抱歉。）

- Jess

3

你可能需要使用函数 plot 和 cumsum。 - sebastian-c

2个回答

3

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

- Wes McClintick

网页内容由stack overflow 提供, 点击上面的

可以查看英文原文，
原文链接

- Paul Hiemstra · Accepted Answer

x的每个元素都代表一个步骤，因此总行程即S(N)的值是所有元素的和。此外，在每个S(p)的行程中，该点等于随机行走向量x累加到索引p的和。

使用plot和cumsum可以生成显示S(n)随时间变化的图形：

enter image description here

不想失礼，但我认为这些基本问题通常最适合由您的主管或老师来帮助解决。