在函数中使用 mapply(...)
Using mapply inside a function(...)
我正在尝试生成多个随机数系列,每个系列都有一组不同的参数。然后我想复制所有这些。我想要一个函数来执行此操作,它似乎应该是一个 1-liner,带有 rnorm()
之类的东西来生成数字,mapply()
来为几个参数集生成这些数字,并且 replicate()
重复整个过程。
为了实现这一点,我写了一个函数,它的 ...
个参数我想对应于 mapply
中的 ...
,但似乎 ...
被忽略了通过 mapply
。我有一些示例可以阐明我想要的输出以及问题:
# Some values to be passed to the ... of mapply
mus <- c(-30, 0, 30)
sds <- c(0.1, 1, 10)
# =================================
# = What I want, if it would work =
# =================================
f <- function(n=3, ...){
replicate(2, mapply(rnorm, MoreArgs=list(n=n), ...))
}
f(mean=mus, sd=sds) # Doesn't work :(
>, , 1
>
> [,1]
>[1,] -0.4901243
>[2,] 0.8268027
>[3,] -0.4829781
>
>, , 2
>
> [,1]
>[1,] -0.02025903
>[2,] -1.57011537
>[3,] 0.49234503
# ==========================================================
# = What I can currently get to work, but not style I want =
# ==========================================================
# I can't understand why f.inside is needed!
f2 <- function(n=3, ...){
f.inside <- function() mapply(rand.gen, MoreArgs=list(n=n), ...)
replicate(2, f.inside())
}
f2(mean=mus, sd=sds) # Desired output
>, , 1
>
> [,1] [,2] [,3]
>[1,] -29.83762 -0.06138165 9.956601
>[2,] -30.04880 1.39123405 33.036675
>[3,] -29.94070 -1.15741513 19.337497
>[4,] -29.92360 0.74300731 38.741367
>[5,] -29.81723 0.84565813 22.261605
>
>, , 2
>
> [,1] [,2] [,3]
>[1,] -30.01407 -0.5198845 20.85942
>[2,] -29.77586 -0.1705062 22.06274
>[3,] -29.96901 -0.4412471 21.42849
>[4,] -30.04079 0.4230790 28.35480
>[5,] -30.04794 0.3000821 50.09012
我想我需要把 ...
包裹起来;我曾随意尝试过 alist()
、list()
、bquote()
、expression()
、as.call()
和其他几种失败的方法。
1) 为什么 f()
中的 mapply
似乎完全忽略了 ...
?
编辑:这是一个辅助问题,答案是"dots don't work with replicate"。好的,进入核心问题...
# =====================================
# = OK, Function for Refined Question =
# =====================================
f.edit <- function(n=3, ...){
l <- list(...)
replicate(2,mapply(rnorm, MoreArgs=list(n=n), l))
}
f.edit(mean=mus, sd=sds) # Doesn't work :(
2) 如何将对象的元素拆分为 length(object) 元素,并将它们作为单独的项目传递给 mapply
的 ...
?
我需要改进我的省略号。我真的被这个难住了。
以您正在尝试的方式 f.edit
...
您可以构建参数列表并执行 do.call
和 mapply
function(n=3, ...){
dots<-list(...);
replicate(2,do.call(mapply,c(list(FUN=rnorm),dots,list(MoreArgs=list(n=n)))))
}
或者,提取参数直接调用
function(n=3, ...){
dots<-list(...);
replicate(2,mapply(rnorm,dots$mean,dots$sd,MoreArgs=list(n=n)))
}
在任何一种情况下,您都在正确的环境中强制计算 ...
,就像在 f2
和 中所做的那样。
f.edit
的问题在于 rnorm
被称为 rnorm(n,list(means=mus,sd=sds))
而不是 rnorm(n,means=mu,sd=sds)
。
但是,在这种情况下,您可能会考虑使用 rnorm
已经在其 mean 和 sd 参数上矢量化的事实
function(n,reps,mean,sd) {
m<-length(mean);
aperm(array(rnorm(n*m*reps,mean,sd),c(m,n,reps)),c(2,1,3))
}
问题是您不能将 ...
与 replicate
一起使用,如前所述 here。
我认为您已经找到了最易读的替代方法,其中包含帮助函数 f.inside()
以获取 ...
术语。
另一种选择是避免 replicate
并使用 lapply
(您最终忽略的索引值)
f3 <- function(n=3, ...){
lapply(1:2, function(i, ...) mapply(rnorm, MoreArgs=list(n=n), ...), ...)
}
f3(mean=mus, sd=sds)
我正在尝试生成多个随机数系列,每个系列都有一组不同的参数。然后我想复制所有这些。我想要一个函数来执行此操作,它似乎应该是一个 1-liner,带有 rnorm()
之类的东西来生成数字,mapply()
来为几个参数集生成这些数字,并且 replicate()
重复整个过程。
为了实现这一点,我写了一个函数,它的 ...
个参数我想对应于 mapply
中的 ...
,但似乎 ...
被忽略了通过 mapply
。我有一些示例可以阐明我想要的输出以及问题:
# Some values to be passed to the ... of mapply
mus <- c(-30, 0, 30)
sds <- c(0.1, 1, 10)
# =================================
# = What I want, if it would work =
# =================================
f <- function(n=3, ...){
replicate(2, mapply(rnorm, MoreArgs=list(n=n), ...))
}
f(mean=mus, sd=sds) # Doesn't work :(
>, , 1
>
> [,1]
>[1,] -0.4901243
>[2,] 0.8268027
>[3,] -0.4829781
>
>, , 2
>
> [,1]
>[1,] -0.02025903
>[2,] -1.57011537
>[3,] 0.49234503
# ==========================================================
# = What I can currently get to work, but not style I want =
# ==========================================================
# I can't understand why f.inside is needed!
f2 <- function(n=3, ...){
f.inside <- function() mapply(rand.gen, MoreArgs=list(n=n), ...)
replicate(2, f.inside())
}
f2(mean=mus, sd=sds) # Desired output
>, , 1
>
> [,1] [,2] [,3]
>[1,] -29.83762 -0.06138165 9.956601
>[2,] -30.04880 1.39123405 33.036675
>[3,] -29.94070 -1.15741513 19.337497
>[4,] -29.92360 0.74300731 38.741367
>[5,] -29.81723 0.84565813 22.261605
>
>, , 2
>
> [,1] [,2] [,3]
>[1,] -30.01407 -0.5198845 20.85942
>[2,] -29.77586 -0.1705062 22.06274
>[3,] -29.96901 -0.4412471 21.42849
>[4,] -30.04079 0.4230790 28.35480
>[5,] -30.04794 0.3000821 50.09012
我想我需要把 ...
包裹起来;我曾随意尝试过 alist()
、list()
、bquote()
、expression()
、as.call()
和其他几种失败的方法。
1) 为什么
编辑:这是一个辅助问题,答案是"dots don't work with replicate"。好的,进入核心问题...f()
中的 mapply
似乎完全忽略了 ...
?
# =====================================
# = OK, Function for Refined Question =
# =====================================
f.edit <- function(n=3, ...){
l <- list(...)
replicate(2,mapply(rnorm, MoreArgs=list(n=n), l))
}
f.edit(mean=mus, sd=sds) # Doesn't work :(
2) 如何将对象的元素拆分为 length(object) 元素,并将它们作为单独的项目传递给 mapply
的 ...
?
我需要改进我的省略号。我真的被这个难住了。
以您正在尝试的方式 f.edit
...
您可以构建参数列表并执行
do.call
和mapply
function(n=3, ...){ dots<-list(...); replicate(2,do.call(mapply,c(list(FUN=rnorm),dots,list(MoreArgs=list(n=n))))) }
或者,提取参数直接调用
function(n=3, ...){ dots<-list(...); replicate(2,mapply(rnorm,dots$mean,dots$sd,MoreArgs=list(n=n))) }
在任何一种情况下,您都在正确的环境中强制计算 ...
,就像在 f2
和
f.edit
的问题在于 rnorm
被称为 rnorm(n,list(means=mus,sd=sds))
而不是 rnorm(n,means=mu,sd=sds)
。
但是,在这种情况下,您可能会考虑使用 rnorm
已经在其 mean 和 sd 参数上矢量化的事实
function(n,reps,mean,sd) {
m<-length(mean);
aperm(array(rnorm(n*m*reps,mean,sd),c(m,n,reps)),c(2,1,3))
}
问题是您不能将 ...
与 replicate
一起使用,如前所述 here。
我认为您已经找到了最易读的替代方法,其中包含帮助函数 f.inside()
以获取 ...
术语。
另一种选择是避免 replicate
并使用 lapply
(您最终忽略的索引值)
f3 <- function(n=3, ...){
lapply(1:2, function(i, ...) mapply(rnorm, MoreArgs=list(n=n), ...), ...)
}
f3(mean=mus, sd=sds)