R randu 来自同余生成器 RANDU 的随机数

说明

400 个连续随机数的三元组取自在 VMS 1.5 下运行的 VAX FORTRAN 函数 RANDU。

用法

randu

格式

一个 DataFrame ，对名为 x 、 y 和 z 的 3 个变量进行 400 个观察，给出三元组中的第一个、第二个和第三个随机数。

细节

在三维显示中，很明显三元组落在 3 空间中的 15 个平行平面上。理论上这对于 RANDU 生成器的所有三元组来说都是正确的。

这些特定的 400 个三元组在序列中以 5 开始，即它们是 ((U[5i+1], U[5i+2], U[5i+3]), i= 0, ..., 399)，并且四舍五入到小数点后 6 位。

在 VMS 2.0 及更高版本中，此问题已得到修复。

例子

## We could re-generate the dataset by the following R code
seed <- as.double(1)
RANDU <- function() {
    seed <<- ((2^16 + 3) * seed) %% (2^31)
    seed/(2^31)
}
for(i in 1:400) {
    U <- c(RANDU(), RANDU(), RANDU(), RANDU(), RANDU())
    print(round(U[1:3], 6))
}

来源

大卫·多诺霍