3

logistf::logistf我正在使用and函数编写一个logistf::forward函数。我将在这里给出一个使用包中的最小工作sex2示例logistf

data(sex2)
fwSel <- function(datamod) {
  fitnull <- logistf(case ~ 1, data = datamod, pl = FALSE) 
  fw <- forward(fitnull)
  return(fw)
}
fwSel(sex2)

我得到以下输出:

步骤 0:启动模型

is.data.frame(data) 中的错误:找不到对象“datamod”。

有人对此有解释吗?

4

1 回答 1

6

这是您可以在 R 中遇到的典型错误。它已被再次询问,不幸的是,它根据不同函数在不同环境中的工作方式以及函数如何尝试根据使用来查找数据parent.envor parent.frame。这可能是以下两个问题之一:

  • 延迟评估问题:
    尝试force(datamod)在您的logistf函数之前使用,因为您的 datamod 当前未在您的自定义函数中评估。如果存在以下问题,这可能不起作用:
  • 您的 datamod 数据集存在于函数的执行环境中。如果函数链中的一个函数使用对 aparent.frame()或对 的调用parent_env(),这将导致问题,因为 R 在不同环境中查找数据的方式不同。解决此问题的唯一方法是在全局环境中启动 datamod,即:

data(sex2)
datamod <- sex2
fwSel <- function(datamod) {
    fitnull <- logistf(case ~ 1, data = datamod, pl = FALSE) 
    fw <- forward(fitnull)
    return(fw)
}
fwSel(sex2)

这肯定会起作用,因为无论如何都会搜索全局环境。

我发现这个链接是一种很好的方式来了解它与parent.env它的不同之处parent.frame以及使用这两个内部函数会如何导致你所面临的问题。

我根据链接中的函数制作了一个新示例,可以准确地演示您的问题:

f = function() {
  print(c(f=environment(), defined_in=parent.env(environment()),  
    called_from=parent.frame()))

  #search for mydata in calling environment
  try(print(get('mydata',env=parent.frame())))  

  #search for mydata in parent evnironment
  try(print(get('mydata',env=parent.env(environment())))) 
  }

g = function() {
  mydata <- c(1,2,3)
  print(c(g=environment(), f()))
  }   

> g()
$f
<environment: 0x0000000030868df8>

$defined_in
<environment: R_GlobalEnv>

$called_from
<environment: 0x000000003086a360>

#the first get works perfect
[1] 1 2 3       

#the second produces an error
Error in get("mydata", env = parent.env(environment())) : 
  object 'mydata' not found

$g
<environment: 0x000000003086a360>

正如您在上面看到的,get与调用环境一起使用是有效的,而get与父环境一起使用会失败并产生错误。这也是(可能)在您的函数中发生的事情。

于 2015-02-19T11:49:26.200 回答