我正在尝试使用参数为 1 的指数分布生成 1000 个数字。
将种子值设置为 1 后,我尝试了rexp(1000, 1)和replicate(1000, rexp(1, 1)),但是得到的两个向量的中位数不同。
我预计这两个表达式生成的向量是相同的,因为它们都是从相同种子值下的相同指数分布中采样的。
rexp(1000, 1)和 和有什么区别replicate(1000, rexp(1, 1))?在实践中我应该使用哪个?
这是我尝试过的代码:
> options(digits = 2)
> set.seed(1)
>
> a <- rexp(1000, 1)
> b <- replicate(1000, rexp(1, 1))
>
> median(a)
[1] 0.73
> median(b)
[1] 0.68
Run Code Online (Sandbox Code Playgroud)