r - 无法将参数传递给自己的函数中的函数

Question

我正在编写一个函数，我想在其中将一些参数传递给 crrstep 函数('crrstep' 包)，但我遇到了一个问题:不知何故，当我输入函数时无法识别函数中的参数“事件”在 crrstep.我猜 crrstep 看起来与我想要它看起来的环境不同，但即使在网上搜索解决方案数小时后，我似乎也无法弄清楚如何解决这个问题(我在编程方面非常缺乏经验..) .任何帮助将不胜感激!

这是一些模拟数据(来自 crrstep-documentation 的调整示例)和我的代码示例:

n <- 500
ftime <- rexp(n)
fstatus <- sample(0:2,n,replace=TRUE)
testdata <- matrix(runif(8*n),nrow=n)
testdata <- cbind(ftime,fstatus,testdata)
dimnames(testdata)[[2]] <- c('ftime','fstatus','x1','x2','x3','x4','x5','x6','x7','x8')
testdata <- as.data.frame(testdata)
formula1 <- ftime ~ 1 + x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8
rm(fstatus,ftime,n)

test.fun <- function(x,data,event){
require(crrstep)
select.mod<- crrstep(formula=x,,etype=event, failcode=1, cencode=0,data=data, 
                   direction = "backward", criterion = "AIC", crr.object = TRUE, 
                   trace = FALSE)
#Rest of function omitted for now
print(select.mod)
}

#Test
test.fun(x=formula1,data=testdata,event=fstatus) 
#I get: Error in eval(expr, envir, enclos) : object 'event' not found"

非常感谢!罗布

Answer 1

在依赖于评估数据框中的名称的函数中调用函数时，我使用 do.call ，它在传递给函数之前评估它的参数，因此使调试和编写代码都更简单，而且我觉得我可以更确定它在做什么。 (对于调试，只需使用 call 而不是 do.call ，这将显示该函数将尝试运行的内容；语法也略有不同，因此在执行此操作时也会删除调用中的列表结构。)

(感谢 Josh O'Brien 对这个想法的回答:https://stackoverflow.com/a/7668846/210673)

在这种情况下，它看起来像这样:

test.fun <- function(x, data, event){
  require(crrstep)
  select.mod <- do.call("crrstep", 
         list(formula=x, etype=substitute(event), failcode=1, cencode=0,
              data=as.name("data"), direction = "backward", criterion = "AIC", 
              crr.object = TRUE, trace = FALSE))
  print(select.mod)
}

test.fun(x=formula1, data=testdata, event=fstatus) 

substitute(event)告诉它使用赋予函数的名称，而不是名称 event . as.name("data")告诉它寻找 data在函数内而不是传递实际的数据框。另一个选项是 substitute(data)它将查找您拥有的实际数据框。

一个使用 lm 的例子

这是使用 lm 的非常相似行为的示例和 weights参数:

这是一个示例数据集和对 lm 的调用，不在另一个函数中。我打印 call响应的元素以查看它实际做了什么。

> set.seed(5)
> dd <- data.frame(x=1:10,y=round(rnorm(10,mean=10),1), z=round(runif(10,1,4),1))
> lm(y~x, weights=z, data=dd)$call
lm(formula = y ~ x, data = dd, weights = z)

自然的方式，它不起作用，因为它正在寻找 w在数据框中:

> f1 <- function(f,w,d){
+   lm(formula=f,weights=w, data=d)
+ }
> f1(y~x, z, dd)
Error in eval(expr, envir, enclos) : object 'w' not found

可以用字符串构建调用；这更简单一些:

> f2 <- function(f,w,d){
+   do.call("lm", list(formula=as.formula(f), weights=as.name(w), data=as.name(d)))
+ }
> f2("y~x", "z", "dd")$call
lm(formula = y ~ x, data = dd, weights = z)

或者可以使用 substitute ;在这里，我在我的实际数据集上调用函数 dd , 而不是 d函数内。如果我想使用 update，这以后可能会有用。 .

> f3 <- function(f,w,d){
+   do.call("lm", list(formula=f, weights=substitute(w), data=substitute(d)))
+ }
> f3(y~x, z, dd)$call
lm(formula = y ~ x, data = dd, weights = z)

但我也可以使用 d在函数内；这次注意data = d在通话中而不是 data = dd .

> f4 <- function(f,w,d){
+   do.call("lm", list(formula=f, weights=substitute(w), data=as.name("d")))
+ }
> f4(y~x, z, dd)$call
lm(formula = y ~ x, data = d, weights = z)

它也可以放入实际数据框，但调用时间更长。但是，如果您在每次调用之前以编程方式更改数据框并希望记录该数据框是什么，则可能需要这样做。 (不过，我更喜欢以更明确的方式保存该数据框，如果您以后确实需要的话。)

> f5 <- function(f,w,d){
+   do.call("lm", list(formula=f, weights=substitute(w), data=d))
+ }
> f5(y~x, z, dd)$call
lm(formula = y ~ x, data = list(x = 1:10, y = c(9.2, 11.4, 8.7, 
10.1, 11.7, 9.4, 9.5, 9.4, 9.7, 10.1), z = c(3.7, 3.2, 1.6, 1.7, 
1.4, 2.4, 2.3, 3.9, 1.4, 3.9)), weights = z)

还有一个表明你不能只使用 substitute没有do.call作为substitute在对 lm 的调用中执行.

> f6 <- function(f,w,d){
+   lm(formula=f,weights=substitute(w), data=d)
+ }
> f6(y~x, z, dd)
Error in model.frame.default(formula = f, data = d, weights = substitute(w),  : 
  invalid type (symbol) for variable '(weights)'