Question

我有这个自定义功能：

f <- function(N, zero_per, sed) {
  set.seed(sed)
  M <- N-(N*zero_per)
  y_g <- rgamma(M, shape=rnorm(M,1,2), rate=rnorm(M,1,2))
  df_g <- data.frame(y=y_g)
  df_l <- data.frame(y=rep(0, N-M))
  df <- rbind(df_g, df_l)
  return(df)
}

这个函数是循环的一部分，我在第二个参数中使用了一个序列：

> sq <- seq(0.05, 0.7, 0.05)
> sq
 [1] 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 0.70

第12位的值为0.6：

> sq[12]
[1] 0.6

然而，当我调用函数时，一次使用数字参数0.6，一次使用sq [12]，我会得到不同的结果：

> head(f(5000, sq[12], 1))
          y
1       NaN
2       NaN
3       NaN
4 0.7547084
5 1.0643216
6       NaN
Warning message:
In rgamma(M, shape = rnorm(M, 1, 2), rate = rnorm(M, 1, 2)) : NAs produced

> head(f(5000, 0.6, 1))
          y
1       NaN
2       NaN
3       NaN
4 0.8990203
5 1.0719659
6       NaN
Warning message:
In rgamma(M, shape = rnorm(M, 1, 2), rate = rnorm(M, 1, 2)) : NAs produced

我使用相同的种子。

奇怪的是，当我使用0.5时，这不会发生：

> head(f(5000, sq[10], 1))
          y
1       NaN
2       NaN
3       NaN
4 0.4426849
5       NaN
6       NaN
Warning message:
In rgamma(M, shape = rnorm(M, 1, 2), rate = rnorm(M, 1, 2)) : NAs produced
> head(f(5000, 0.5, 1))
          y
1       NaN
2       NaN
3       NaN
4 0.4426849
5       NaN
6       NaN
Warning message:
In rgamma(M, shape = rnorm(M, 1, 2), rate = rnorm(M, 1, 2)) : NAs produced

警告可以忽略，它们只是因为我剥夺了我的真实功能而发生。

这里发生了什么？

Answer 1

问题在于：

sq[12]==0.6
#[1] FALSE

虽然

sq[10]==0.5
#[1] TRUE

要查看，请执行以下操作：

options(digits = 20)

sq[12]
#[1] 0.60000000000000009

sq[10]
#[1] 0.5

使用索引向量参数调用自定义函数

1 个答案: