Question

我构建了以下mwe，它显示了一系列嵌套函数，这些函数在foreach语句中调用，MPI作为后端。所有这些都包装在R包中。问题是，当我加载包，然后调用fun3时，我无法将fun3on的值从TRUE更改为FALSE。有什么建议？完整的代码是 https://github.com/jamaas/toymod3.git

期待调用

fun1（fun3on = FALSE）

更改输出值，但它没有？

#' Test function level 1
#' @param var11 first variable for function 1
#' @param var12 second variable for function 1
#' @param var13 third variable for function 1
#' @export fun1

fun1 <- function (var11=10, var12=8, var13=5,
                  var21=0.05, var22=9.876,
                  fun3on=TRUE, ...) {

    results <- data.frame (foreach::`%dopar%`(
            foreach::`%:%`(foreach::foreach(j = 1:var11,
                                            .combine = cbind,
                                            .packages="toymod2",
                                            .export = "fun3on"),
            foreach::foreach (i = 1:var12, .combine=rbind)),
                   out3 <- replicate(var13, do.call(fun2, list(var21, var22)))
        )
)
## save outputs for subsequent analyses if required
saveRDS(results, file = paste("./outputs/", var13 ,"_", var12, "_", var11, "_",
                              format(Sys.time(), "%d_%m_%Y"), ".rds", sep=""))
return(results)
}

#' Test function level 2
#' @param var21 first variable for function 2
#' @param var22 second variable for function 2
#' @export fun2

fun2 <- function (var21=0.05, var22=9.876, ...) {
    out2 <- `if` (rpois(1, var21) > 0,
                    var22 * fun3(...),
                  0)
}

#' Test function level 3
#' @param var31 first variable for function 3
#' @param var32 second variable for function 3
#' @param var33 third variable for function 3
#' @param fun3on turn the formula on or off
#' @export fun3

fun3 <- function (var31=1.396, var32=14.387,
                  var33=3.219, fun3on = TRUE, ...) {
    out3 <- `if` (fun3on,
                  var31 * rnorm(1, mean=var32, sd= var33),
                  500)
}

Answer 1

问题在于您没有将fun3on传递给fun2，因此无法将其传递给fun3。

在do.call(fun2, list(var21, var22))中，将其更改为do.call(fun2, list(var21, var22, fun3on = fun3on, ...))。确保将它们作为命名参数传递（我添加了点，因为您似乎在fun3中有一些。）

PS：除非你使用令人回味的名字，否则不要给你的回报命名。

Answer 2

看起来很复杂。我发表评论但由于这个可能的答案较长并且代表了回答您问题的尝试，所以

我猜你的问题就在这里

在fun1中，您拨打fun2

fun2 <- function (var21=0.05, var22=9.876, ...) {
    out2 <- `if` (rpois(1, var21) > 0,
                    var22 * fun3(...),
                  0)
}

fun2调用fun3(...)而不传递参数

fun3 <- function (var31=1.396, var32=14.387,
                  var33=3.219, fun3on = TRUE, ...) {
    out3 <- `if` (fun3on,
                  var31 * rnorm(1, mean=var32, sd= var33),
                  500)
}

但fun3on

中的TRUE是默认fun3

也许您将(...)更改为(fun3on=fun3on)

fun2 <- function (var21=0.05, var22=9.876, ...) {
    out2 <- `if` (rpois(1, var21) > 0,
                    var22 * fun3(fun3on=fun3on),
                  0)
}

它可能有用???

无法使用doMPI更改foreach函数内的布尔变量值

2 个答案: