在R中获得mapply函数的均值

Question

我在这里有此代码：

mat <- cbind(shape,scale_ima)
mapply(rweibull, 10, mat[,"shape"], mat[,"scale_ima"])

mapply函数遍历mat矩阵上的每个条目，对于每个循环，从rweibull函数获得一个包含10个条目的向量。 我需要找到向量的均值。

之后，我需要将每个均值存储在一个新向量中，该怎么办？ 因此，对于矢量的每个均值，我需要从mapply中获取数据，我需要将数据存储在矢量中。

这不起作用：

mapply(mean(rweibull, 10, mat[,"shape"], mat[,"scale_ima"]))

Answer 1

您可以执行以下操作找到每个向量的均值并存储这些均值。

 result<-mapply(rweibull, 10, mat[,"shape"], mat[,"scale_ima"])

 #get and store the means
 meanStorage<-apply(result,2, mean)

为了解释为什么对mapply的调用不起作用，请查看mapply函数的参数。 第一个参数是一个函数，第二个至n（由...表示）参数是要传递给该函数的参数。 通过致电

mapply(mean(....))

您传递了均值函数以进行映射，但没有任何要均值的参数。

Answer 2

不确定您的预期结果：

set.seed(42)
colMeans(mapply(rweibull, n = 10L, shape = mat[,"shape"], scale = mat[,"scale_ima"]))
#[1] 0.08762851 0.11618909 0.10013489 0.24263038 0.12706183 0.18069583
#[7] 0.10953056 0.39062035 0.24496276 0.23653724

要么

set.seed(42)
mapply(function(x,y) mean(rweibull(10, x,y)), mat[,"shape"], mat[,"scale_ima"])
#[1] 0.08762851 0.11618909 0.10013489 0.24263038 0.12706183 0.18069583
# [7] 0.10953056 0.39062035 0.24496276 0.23653724

数据

mat <- structure(c(0.953866743, 0.939544872, 0.88055226, 0.937567804, 
0.902443856, 0.969984293, 0.953468872, 0.929905045, 0.889375987, 
0.910115923, 0.152704576, 0.168592082, 0.13059434, 0.153850643, 
0.172734767, 0.162162429, 0.172533372, 0.160826152, 0.190843263, 
0.156289128), .Dim = c(10L, 2L), .Dimnames = list(NULL, c("shape", 
"scale_ima")))