1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
|
source("incl/start.R")
library("datasets") ## cars data set
library("stats") ## lm(), poly(), xtabs()
plan(batchtools_local)
## CRAN processing times:
## On Windows 32-bit, don't run these tests on batchtools
if (!fullTest && isWin32) plan(sequential)
message("*** Globals - formulas ...")
message("*** Globals - lm(<formula>) ...")
## From example("lm", package = "stats")
ctl <- c(4.17, 5.58, 5.18, 6.11, 4.50, 4.61, 5.17, 4.53, 5.33, 5.14)
trt <- c(4.81, 4.17, 4.41, 3.59, 5.87, 3.83, 6.03, 4.89, 4.32, 4.69)
group <- gl(2, 10, 20, labels = c("Ctl", "Trt"))
weight <- c(ctl, trt)
## Truth:
fit0 <- lm(weight ~ group - 1)
print(fit0)
## Explicit future
f <- future({ lm(weight ~ group - 1) })
fit <- value(f)
print(fit)
stopifnot(all.equal(fit, fit0))
## Future assignment
fit %<-% { lm(weight ~ group - 1) }
print(fit)
stopifnot(all.equal(fit, fit0))
message("*** Globals - lm(<formula>) ... DONE")
message("*** Globals - one-side formulas, e.g. xtabs(~ x) ...")
x <- c(1, 1, 2, 2, 2)
## Truth:
tbl0 <- xtabs(~ x)
print(tbl0)
## Explicit future
f <- future({ xtabs(~ x) })
tbl <- value(f)
print(tbl)
stopifnot(all.equal(tbl, tbl0))
## Future assignment
tbl %<-% { xtabs(~ x) }
print(tbl)
stopifnot(all.equal(tbl, tbl0))
message("*** Globals - one-side formulas, e.g. xtabs(~ x) ... DONE")
message("*** Globals - lm(<formula>, data = cars) ...")
exprs <- list(
# "remove-intercept-term" form of no-intercept
a = substitute({ lm(dist ~ . - 1, data = cars) }),
# "make-intercept-zero" form of no-intercept
b = substitute({ lm(dist ~ . + 0, data = cars) }),
# doesn't do what we want here
c = substitute({ lm(dist ~ speed + speed ^ 2, data = cars) }),
# gets us a quadratic term
d = substitute({ lm(dist ~ speed + I(speed ^ 2), data = cars) }),
# avoid potential multicollinearity
e = substitute({ lm(dist ~ poly(speed, 2), data = cars) })
)
for (kk in seq_along(exprs)) {
expr <- exprs[[kk]]
name <- names(exprs)[kk]
mprintf("- Globals - lm(<formula #%d (%s)>, data = cars) ...\n",
kk, sQuote(name))
fit0 <- eval(expr)
print(fit0)
f <- future(expr, substitute = FALSE)
fit <- value(f)
print(fit)
stopifnot(all.equal(fit, fit0))
} ## for (kk ...)
message("*** Globals - lm(<formula>, data = cars) ... DONE")
message("*** Globals - map(x, ~ expr) ...")
## A fake purrr::map() function with limited functionality
map <- function(.x, .f, ...) {
if (inherits(.f, "formula")) {
expr <- .f[[-1]]
.f <- eval(bquote(function(...) {
.(expr)
}))
}
eval(lapply(.x, FUN = .f, ...))
}
inner_function <- function(x) { x + 1 }
outer_function <- function(x) {
map(1:2, ~ inner_function(.x))
}
y0 <- outer_function(1L)
str(y0)
f <- future({ outer_function(1L) })
y <- value(f)
str(y)
stopifnot(all.equal(y, y0))
y %<-% { outer_function(1L) }
str(y)
stopifnot(all.equal(y, y0))
message("*** Globals - map(x, ~ expr) ... DONE")
message("*** Globals - formulas ... DONE")
source("incl/end.R")
|
