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))
}

來源

大衛·多諾霍