Как оценить вызов в 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)
