根据另外两个矩阵中的元素计算矩阵中的每个元素

Question

我正在尝试根据存储在另外两个矩阵中的平均值和标准差来填充来自随机正态分布的i x j个条目的矩阵。有没有办法使用rnorm从两个“数据”矩阵（带有均值和标准差的两个矩阵）中拉出每个条目而不使用循环？

Answer 1

当然，就这样做：

means <- matrix(1:4, 2, 2)
sds <- matrix((1:4)/1000, 2, 2)
result <- matrix(rnorm(4, mean = means, sd = sds), 2, 2)

或（根据以下Frank的评论）

result <- array(rnorm(length(means), mean = means, sd = sds),
                dim = dim(means))