使用均值向量调用 rnorm

Question

当我调用 rnorm 并将单个值作为平均值传递时，很明显会发生什么：从 Normal(10,1) 生成一个值。

y <- rnorm(20, mean=10, sd=1)

但是，我看到整个向量被传递给 rnorm （或 rcauchy 等）的示例；在这种情况下，我不确定 R 机器到底做了什么。例如：

a = c(10,22,33,44,5,10,30,22,100,45,97)
y <- rnorm(a, mean=a, sd=1)

有什么想法吗？

Answer 1

rnorm 生成的随机数数量等于 a 的长度。来自?rnorm：

n：观察次数。如果
'长度(n)> 1'，取长度
为所需数量。

要了解将 a 传递给 Mean 参数时发生的情况，如果我们更改示例会更容易：

a = c(0, 10, 100)
y = rnorm(a, mean=a, sd=1)
[1] -0.4853138  9.3630421 99.7536461

因此，我们生成具有平均值 的 length(a) 随机数a[i]。

Answer 2

一个更好的例子：

a <- c(0,10,100)
b <- c(2,4,6)
y <- rnorm(6,a,b)
y

结果

[1]  -1.2261425  10.1596462 103.3857481  -0.7260817   7.0812499  97.8964131

如您所见，对于y的第一个和第四个元素，rnorm采用<的第一个元素em>a 作为 mean 和 b 的第一个元素sd。

对于y的第二个和第五个元素，rnorm采用a的第二个元素> 作为 mean，b 的第二个元素作为 sd。

对于y的第三个和第六个元素，rnorm采用a的第三个元素> 作为mean，b的第三个元素作为sd。

您可以在 rnorm 的第一个参数中尝试不同的数字，看看发生了什么

例如，如果您在调用 rnorm 时使用 5 而不是 6 作为第一个参数，会发生什么？