r - runif 以 .Machine$double.xmax 作为边界

Question

我想生成一个随机的实数（我猜是有理数）。

为此，我想使用runif(1, min = m, max = M)并且我的想法是设置m; M（绝对）尽可能大，以使间隔尽可能大。这让我想到了我的问题：

M <- .Machine$double.xmax
m <- -M
runif(1, m, M)
## which returns
[1] Inf

为什么它不返回一个数字？选择的间隔是否太大？

附言

> .Machine$double.xmax
[1] 1.797693e+308

Answer 1

正如 mt1022 所暗示的，原因在于runif C 源代码：

double runif(double a, double b)
{
    if (!R_FINITE(a) || !R_FINITE(b) || b < a)  ML_ERR_return_NAN;

    if (a == b)
    return a;
    else {
    double u;
    /* This is true of all builtin generators, but protect against
       user-supplied ones */
    do {u = unif_rand();} while (u <= 0 || u >= 1);
    return a + (b - a) * u;
    }
}

在return参数中，您可以看到a + (b - a) * u在用户提供的区间 [a, b] 中均匀变换 [0, 1] 生成随机值的公式。在您的情况下，它将是-M + (M + M) * u. 所以M + M万一它1.79E308 + 1.79E308生成Inf. 即finite + Inf * finite = Inf：

M + (M - m) * runif(1, 0, 1)
# Inf