【发布时间】:2022-02-19 23:47:22
【问题描述】:
如何在 R 中模拟概率分布为 P(Xn=1)=1/2, P(Xn=-1)=1/2 的随机变量 Xn 序列? (用于模拟随机游走)。
【问题讨论】:
-
你可以像
ifelse(runif(1000) < 0.5,-1,1)一样使用ifelse
标签: r simulation probability random-walk
如何在 R 中模拟概率分布为 P(Xn=1)=1/2, P(Xn=-1)=1/2 的随机变量 Xn 序列? (用于模拟随机游走)。
【问题讨论】:
ifelse(runif(1000) < 0.5,-1,1)一样使用ifelse
标签: r simulation probability random-walk
您可以简单地采样所需数量的元素:
n_decisions <- 100
set.seed(1)
random_walk <- sample(c(1, -1), size = n_decisions, replace = TRUE)
table(random_walk)
random_walk
-1 1
51 49
cumsum(random_walk)
[1] 1 0 1 2 1 2 3 4 3 2 3 4 5 6 7 6 5 4 3 4 5 6 7 8 9
[26] 10 9 10 11 10 9 8 9 8 9 10 9 10 9 8 7 6 7 6 5 4 3 2 3 4
[51] 3 4 3 2 3 4 3 2 1 2 3 2 1 0 -1 -2 -3 -2 -3 -4 -5 -6 -5 -4 -3
[76] -4 -5 -4 -3 -4 -5 -6 -5 -4 -3 -4 -3 -4 -3 -4 -3 -2 -3 -4 -3 -2 -1 -2 -3 -2
绘制 N = 150.000
【讨论】: