A Future for R: Common Issues with Solutions

A Future for R: Common Issues with Solutions

In the ideal case, all it takes to start using futures in R is to replace select standard assignments (<-) in your R code with future assignments (%<-%) and make sure the right-hand side (RHS) expressions are within curly brackets ({ ... }). Also, if you assign these to lists (e.g. in a for loop), you need to use a list environment (listenv) instead of a plain list.

However, as show below, there are few cases where you might run into some hurdles, but, as also shown, they are often easy to overcome. These are often related to global variables.

If you identify other cases, please consider reporting them so they can be documented here and possibly even be fixed.

Issues with globals and packages

Missing globals (false negatives)

If a global variable is used in a future expression that conditionally overrides this global variable with a local one, the future framework fails to identify the global variable and therefore fails to export it, resulting in a run-time error. For example, although this works:

plan(multisession)

reset <- TRUE
x <- 1
y %<-% { if (reset) x <- 0; x + 1 }
y
## [1] 1

the following does not work:

reset <- FALSE
x <- 1
y %<-% { if (reset) x <- 0; x + 1 }
y
## Error: object 'x' not found

It is recommended to avoid above constructs where it is ambiguous whether a variable is global or local. To force variable x to always be global, insert it at the very being of the future expression, e.g.

reset <- FALSE
x <- 1
y %<-% { x; if (reset) x <- 0; x + 1 }
y
## [1] 2

Comment: The goal is to in a future version of the package detect globals also in expression where the local-global state of a variable is only known at run time.

do.call() - function not found

When calling a function using do.call() make sure to specify the function as the object itself and not by name. This will help identify the function as a global object in the future expression. For instance, use

do.call(file_ext, list("foo.txt"))

instead of

do.call("file_ext", list("foo.txt"))

so that file_ext() is properly located and exported. Although you may not notice a difference when evaluating futures in the same R session, it may become a problem if you use a character string instead of a function object when futures are evaluated in external R sessions, such as on a cluster. It may also become a problem with futures evaluated with lazy evaluation if the intended function is redefined after the future is resolved. For example,

> library(future)
> library(listenv)
> library(tools)
> plan(sequential)
> pathnames <- c("foo.txt", "bar.png", "yoo.md")
> res <- listenv()
> for (ii in seq_along(pathnames)) {
+   res[[ii]] %<-% do.call("file_ext", list(pathnames[ii])) %lazy% TRUE
+ }
> file_ext <- function(...) "haha!"
> unlist(res)
[1] "haha!" "haha!" "haha!"

get() - object not found

The base R function get() can be used to get the value of a object by its name. For example, get("pi", envir = baseenv()) will get the value of object pi in the ‘base’ environment, i.e. it corresponds to base::pi. If no other objects named pi exists on the search path, we could have used get("pi") and pi, respectively. It is not unusual to see code snippets such as:

> a <- 1:3
> b <- 4:6
> c <- 3:5
> my_sum <- function(var) { sum(get(var)) }
> y <- my_sum("a")
> y
[1] 6

If we attempt to call my_sum() via a future, we will get an error (if the future is resolved in an external R process);

> library(future)
> plan(multisession)
> f <- future(my_sum("a"))
> y <- value(f)
Error in get(var) : object 'a' not found

This is because the static code inspection done on the future expression my_sum("a") does not reveal object a as a global object. In that expression alone, there are only three objects: the function my_sum(), the primitive function (, and the string "a", and none of those are object a. The future framework will also scan these three objects for globals, which in this example means that it scans also my_sum(). This recursive search for globals will identify three additional globals, namely, the primitive function {, the function sum(), and the function get(), but, as before, none of these source will identify a as a global object. In order for a to be identified, the future framework would need to have a built-in understanding on how get(var) works, which would be a daunting task, especially if it need to know how it acts for different data types of var and various choices on arguments envir and enclos. In fact, this can often not be inferred until run time, that is, it is not possible to identify what objects are needed without actually running the code. In short, it is not possible to automatically identify global variables specified via a character string.

The workaround is to tell the future framework what additional globals are needed. This can be done via argument globals using:

> f <- future(my_sum("a"), globals = structure(TRUE, add = "a"))
> y <- value(f)
> y
[1] 6

or by injecting variable a at the beginning of the future expression, e.g.

> f <- future({ a; my_sum("a") })
> y <- value(f)
> y
[1] 6

Note that, independently of the future framework, it is often a bad idea to use get(), and related functions mget() and assign(), in R code. Searching the archives of R forums, such as the R-help and R-devel mailing lists, will reveal numerous suggestions against using them. A good rule of thumb is:

If you find yourself using get() in your code, take a step back, and reconsider your implementation. There is most likely a better solution available.

For example, consider this, slightly more complex, example:

> a <- 1:3
> b <- 4:7
> c <- 3:5
> my_sum <- function(var) { sum(get(var)) }
> y <- sapply(c("a", "b", "c"), FUN = my_sum)
> y
 a  b  c
 6 22 12

Instead of using “free roaming” objects a, b, and c, it’s better to put those values in a list (or a data frame of of the same length);

> data <- list(a = 1:3, b = 4:7, c = 3:5)

This will in turn allow us to perform the same calculations without having to use get();

> my_sum <- function(x) { sum(x) }
> y <- sapply(data, FUN = my_sum)
> y
 a  b  c
 6 22 12

glue::glue() - object not found

The future framework will fail to identify globals that are declared via character strings. The above section gives an example of this where get() is used and explains why it is not feasible to automatically identify string-embedded globals from such code. Another example, is when using glue() from the glue package to generate strings dynamically, e.g.

> library(glue)
> a <- 42
> s <- glue("The value of a is {a}.")
> s
The value of a is 42.

Attempt to perform the same via a future that is resolved in another R session will produce an “object not found” error;

> library(glue)
> library(future)
> plan(multisession)
> a <- 42
> s %<-% glue("The value of a is {a}.")
> s
Error in eval(parse(text = text, keep.source = FALSE), envir) : 
  object 'a' not found

As explained in the previous section, the workaround is to specify what additional global variables there are, which can be done as:

> s %<-% glue("The value of a is {a}.") %globals% structure(TRUE, add = "a")
> s
The value of a is 42.

An alternative solution is to guide the future framework by adding the missing globals as “dummy” variables, e.g.

> s %<-% { a; glue("The value of a is {a}.") }
> s
The value of a is 42.

Missing packages (false negatives)

Occasionally, the static-code inspection of the future expression fails to identify packages needed to evaluated the expression. This may occur when an expression uses S3 generic functions part of one package whereas the required S3 method is in another package. For example, in the below future generic function [ is used on data.table object DT, which requires S3 method [.data.table from the data.table package. However, the future and globals packages fail to identify data.table as a required package, which results in an evaluation error:

> library(future)
> plan(multisession)

> library(data.table)
> DT <- data.table(a = LETTERS[1:3], b = 1:3)
> y %<-% DT[, sum(b)]
> y
Error: object 'b' not found

The above error occurs because, contrary to the master R process, the R worker that evaluated the future expression does not have data.table loaded. Instead the evaluation falls back to the [.data.frame method, which is not what we want.

Until the future framework manages to identify data.table as a required package (which is the goal), we can guide future by specifying additional packages needed:

> y %<-% DT[, sum(b)] %packages% "data.table"
> y
[1] 6

or equivalently

> f <- future(DT[, sum(b)], packages = "data.table")
> value(f)
[1] 6

Note, do not use library() or loadNamespace() to resolve these problems. It is always better to use the above packages approach.

Non-exportable objects

Certain types of objects are tied to a given R session and cannot be passed along to another R process (a “worker”). An example of a non-exportable object is is XML objects of the xml2 package. If we attempt to use those in parallel processing, we may get a error when the future is evaluated (or just invalid results depending on how they are used), e.g.

> library(future)
> plan(multisession)
> library(xml2)
> xml <- read_xml("<body></body>")
> f <- future(xml_children(xml))
> value(f)
Error: external pointer is not valid

The future framework can help detect this before sending off the future to the worker;

> options(future.globals.onReference = "error")
> f <- future(xml_children(xml))
Error in FALSE : 
  Detected a non-exportable reference ('externalptr') in one of the globals
('xml' of class 'xml_document') used in the future expression

For additional details on non-exportable objects and examples of other R packages that use objects that may cause problems in parallel processing, see Vignette ‘Non-Exportable Objects’.

Trying to pass an unresolved future to another future

It is not possible for a future to resolve another one unless it was created by the future trying to resolve it. For instance, the following gives an error:

> library(future)
> plan(multisession)
> f1 <- future({ Sys.getpid() })
> f2 <- future({ value(f1) })
> v1 <- value(f1)
[1] 7464
> v2 <- value(f2)
Error: Invalid usage of futures: A future whose value has not yet been collected
 can only be queried by the R process (cdd013cb-e045-f4a5-3977-9f064c31f188; pid
 1276 on MyMachine) that created it, not by any other R processes (5579f789-e7b6
 -bace-c50d-6c7a23ddb5a3; pid 2352 on MyMachine): {; Sys.getpid(); }

This is because the main R process creates two futures, but then the second future tries to retrieve the value of the first one. This is an invalid request because the second future has no channel to communicate with the first future; it is only the process that created a future who can communicate with it(*).

Note that it is only unresolved futures that cannot be queried this way. Thus, the solution to the above problem is to make sure all futures are resolved before they are passed to other futures, e.g.

> f1 <- future({ Sys.getpid() })
> v1 <- value(f1)
> v1
[1] 7464
> f2 <- future({ value(f1) })
> v2 <- value(f2)
> v2
[1] 7464

This works because the value has already been collected and stored inside future f1 before future f2 is created. Since the value is already stored internally, value(f1) is readily available everywhere. Of course, instead of using value(f1) for the second future, it would be more readable and cleaner to simply use v1.

The above is typically not a problem when future assignments are used. For example:

> v1 %<-% { Sys.getpid() })
> v2 %<-% { v1 }
> v1
[1] 2352
> v2
[1] 2352

The reason that this approach works out of the box is because in the second future assignment v1 is identified as a global variable, which is retrieved. Up to this point, v1 is a promise (“delayed assignment” in R), but when it is retrieved as a global variable its value is resolved and v1 becomes a regular variable.

However, there are cases where future assignments can be passed via global variables without being resolved. This can happen if the future assignment is done to an element of an environment (including list environments). For instance,

> library(listenv)
> x <- listenv()
> x$a %<-% { Sys.getpid() }
> x$b %<-% { x$a }
> x$a
[1] 2352
> x$b
Error: Invalid usage of futures: A future whose value has not yet been collected
 can only be queried by the R process (cdd013cb-e045-f4a5-3977-9f064c31f188; pid
 1276 on localhost) that created it, not by any other R processes (2ce86ccd-5854
 -7a05-1373-e1b20022e4d8; pid 7464 on localhost): {; Sys.getpid(); }

As previously, this can be avoided by making sure x$a is resolved first, which can be one in various ways, e.g. dummy <- x$a, resolve(x$a) and force(x$a).

Footnote: (*) Although sequential futures could be passed on to other futures part of the same R process and be resolved there because they share the same evaluation process, by definition of the Future API it is invalid to do so regardless of future type. This conservative approach is taken in order to make future expressions behave consistently regardless of the type of future used.

Miscellaneous

Capturing errors, outputting their messages, and returning a default value

Sometimes a function call produce an error for a particular input. In such cases, we might want to return a default value, say, a missing value, instead of signaling an error. This can be done using:

res <- tryCatch({
  unstable_calc(x)
}, error = function(e) {
  NA_real_
})

Here, res takes the value of unstable_calc(x), unless it produces an error, in case it takes value NA_real_.

In addition to the above, we could produce a warning whenever we get an error and replace it with a missing value. We can do this as:

res <- tryCatch({
  unstable_calc(x)
}, error = function(e) {
  warning(conditionMessage(e))
  NA_real_
})

This will turn the error into a warning with the same message. If we want to just output the message without producing a warning, we can use message(conditionMessage(e)).

Importantly, we must not use just warning(e) or message(e), although it appears to work at a first glance. If we do, we will end up re-signaling the error but without interruption. It is an important distinction that will reveal itself if used within futures. The above example with warning(conditionMessage(e)) will work as expected, but if we use warning(e) the future framework will produce an error, not a warning.

Using source() in a future

Avoid using source() inside futures. It is always better to source external R scripts at the top of your main R script, e.g.

library(future)
source("./my-script.R")

f <- future({
  ...
})

However, if you find yourself having to source a script inside a future, or inside a function, make sure to specify argument local = TRUE, e.g.

f <- future({
  source("./my-script.R", local = TRUE)
  ...
})

This is because source() defaults to local = FALSE, which has side effects. When using local = FALSE, any functions or variables defined by the R script are assigned to the global environment - not the calling environment as we might expect. This may make little different when calling source() from the R prompt, or from another script. However, when called from inside a function, inside local(), or inside a future, it might result in unexpected behavior. It is similar to using assign("a", 42, envir = globalenv()), which is known be a bad practice. To be on the safe side, it is almost always better call source() with local = TRUE.

Clashes with other packages

Sometimes other packages have functions or operators with the same name as the future package, and if those packages are attached after the future package, their objects will mask the ones of the future package. For instance, the igraph package also defines a %<-% operator which clashes with the one in future if used at the prompt or in a script (it is not a problem inside package because there we explicitly import objects in a known order). Here is what we might get:

> library(future)
> library(igraph)

Attaching package: 'igraph'

The following objects are masked from 'package:future':

    %<-%, %->%

The following objects are masked from 'package:stats':

    decompose, spectrum

The following object is masked from 'package:base':

    union

> y %<-% { 42 }
Error in get(".igraph.from", parent.frame()) : 
  object '.igraph.from' not found

Here we get an error because %<-% is from igraph and not the future assignment operator as we wanted. This can be confirmed as:

> environment(`%<-%`)
<environment: namespace:igraph>

To avoid this problem, attach the two packages in opposite order such that future comes last and thereby overrides igraph, i.e.

> library(igraph)
> library(future)

Attaching package: 'future'

The following objects are masked from 'package:igraph':

%<-%, %->%

> y %<-% { 42 }
> y
[1] 42

An alternative is to detach the future package and re-attach it, which will achieve the same thing:

> detach("package:future")
> library(future)

Yet another alternative is to explicitly override the object by importing it to the global environment, e.g.

> `%<-%` <- future::`%<-%`
> y %<-% { 42 }
> y
[1] 42

In this case, it does not matter in what order the packages are attached because we will always use the copy of future::`%<-%`.

Syntax error: “non-numeric argument to binary operator”

The future assignment operator %<-% is a binary infix operator, which means it has higher precedence than most other binary operators but also higher than some of the unary operators in R. For instance, this explains why we get the following error:

> x %<-% 2 * runif(1)
Error in x %<-% 2 * runif(1) : non-numeric argument to binary operator

What effectively is happening here is that because of the higher priority of %<-%, we first create a future x %<-% 2 and then we try to multiply the future (not its value) with the value of runif(1) - which makes no sense. In order to properly assign the future variable, we need to put the future expression within curly brackets;

> x %<-% { 2 * runif(1) }
> x
[1] 1.030209

Parentheses will also do. For details on precedence on operators in R, see Section ‘Infix and prefix operators’ in the ‘R Language Definition’ document.

Error: invalid ‘type’ (environment) of argument with magrittr

Another example where the future assignment operator %<-% requires curly brackets is when using the magrittr infix operator %>%, e.g.

> library(magrittr)
> x %<-% 1:100 %>% sum
Error in sum(.) : invalid 'type' (environment) of argument

The reason for this error is that x %<-% 1:100 is passed to sum() by %>%. To fix this, use:

> x %<-% { 1:100 %>% sum }
> x
[1] 5050

R CMD check NOTEs

The code inspection run by R CMD check will not recognize the future assignment operator %<-% as an assignment operator, which is not surprising because %<-% is technically an infix operator. This means that if you for instance use the following code in your package:

foo <- function() {
  b <- 3.14
  a %<-% { b + 1 }
  a
}

then R CMD check will produce a NOTE saying:

* checking R code for possible problems ... NOTE
foo: no visible binding for global variable 'a'
Undefined global functions or variables:
  a

In order to avoid this, we can add a dummy assignment of the missing global at the top of the function, i.e.

foo <- function() {
  a <- NULL ## To please R CMD check
  b <- 3.14
  a %<-% { b + 1 }
  a
}