Невозможно изменить значение булевой переменной внутри функции foreach с помощью doMPI

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)
}

это может сработать???