Как оценить вызов в data.frame?

В одном пакете, над которым я работаю, я сталкиваюсь с проблемой, которая выглядит просто - но я не могу понять:

Подфункция задается двумя аргументами:

obj, a data.frame
foo, a call

Например:

> head(obj)
  cadmium copper lead zinc  elev
1    11.7     85  299 1022 7.909
2     8.6     81  277 1141 6.983
3     6.5     68  199  640 7.800
4     2.6     81  116  257 7.655
5     2.8     48  117  269 7.480
6     3.0     61  137  281 7.791
> foo
log(cadmium)
> class(foo)
[1] "call"

В этом примере я хочу создать вектор x <- log(obj$cadmium). Как мне это сделать? Я попытался использовать with(), но я не получил ожидаемого результата:

> with(obj, foo)
log(cadmium)

foo - это вызов, созданный пользователем путем задания преобразования в столбце data.frame obj:

my_function(obj, foo = log(cadmium)) { ... }

dput() фрагмента данных:

obj <- structure(list(cadmium = c(11.7, 8.6, 6.5, 2.6, 2.8, 3), copper = c(85L, 
81L, 68L, 81L, 48L, 61L), lead = c(299L, 277L, 199L, 116L, 117L, 
137L), zinc = c(1022L, 1141L, 640L, 257L, 269L, 281L), elev = c(7.909, 
6.983, 7.8, 7.655, 7.48, 7.791)), .Names = c("cadmium", "copper", 
"lead", "zinc", "elev"), class = "data.frame", row.names = c("1", 
"2", "3", "4", "5", "6"))

Ответ 1

Вам нужно оценить вызов, например, используя eval():

foo <- call("log", quote(cadmium))
with(obj, eval(foo))

который дает:

> with(obj, eval(foo))
[1] 2.4595888 2.1517622 1.8718022 0.9555114 1.0296194 1.0986123

где obj - фрагмент данных, которые вы показали.

eval() также имеет аргумент envir, указывающий среду, в которой вычисляется выражение. Таким образом, вы можете делать то, что хотите, без with() напрямую используя eval():

> eval(foo, envir = obj)
[1] 2.4595888 2.1517622 1.8718022 0.9555114 1.0296194 1.0986123

Ответ 2

Что-то вроде

z <- data.frame(a1=1:5,b1=LETTERS[1:5],c1=letters[1:5])
foo <- quote(log(a1))
eval(foo,envir=z)

но будьте осторожны: как только вы начнете использовать eval, вы опускаетесь на более низкие глубины R. Я все еще не полностью понимаю различия между вызывающими кадрами, охватывая фреймы и т.д. и т.д....

Ответ 3

Не похоже, что вы правильно определили функцию:

> foo <- function(x) log(x)
> class(foo)
[1] "function"

Тогда do:

x <- foo(obj$cadmium)

Или:

x <- foo(obj[["cadmium"]])

Или:

x <- with(obj, foo(cadmium)