从另一个函数内部调用mlogit（），在使用attach时确定变量的问题

Question

我需要从另一个函数内部调用mlogit() R function。

这是用于演示目的的功能：

#-------------------------
# DEMO FUNCTION
#-------------------------
# f = formula (string)
# fData = data.frame
# cVar = choice variable (string)
# optVar = alternative variable (string)
##########################
mlogitSum <- function(f, fData,  cVar="choice", optVar="option"){
  library(mlogit)
  r2 <- mlogit(as.formula(f), shape = "long", data = fData, alt.var=optVar, choice = cVar)
  return(summary(r2))
}

显然存在环境问题，因此mlogit（）函数找不到全局声明的变量作为参数。

此示例不起作用：

mydata <- read.csv(url("http://www.ats.ucla.edu/stat/r/dae/mlogit.csv"))
attach(mydata)
library(mlogit)
mydata$brand<-as.factor(mydata$brand)
mlData<-mlogit.data(mydata, varying=NULL, choice="brand", shape="wide")

myFormula <-"brand~1|female+age"
var1 <- "brand"
var2 <- "alt"
mlogitSum(myFormula, fData = mlData, var1, var2)

如果变量是在主环境中分配的，它可以工作：

mydata <- read.csv(url("http://www.ats.ucla.edu/stat/r/dae/mlogit.csv"))
attach(mydata)
library(mlogit)
mydata$brand<-as.factor(mydata$brand)
fData<-mlogit.data(mydata, varying=NULL, choice="brand", shape="wide")

myFormula <-"brand~1|female+age"
cVar <- "brand"
optVar <- "alt"
mlogitSum(myFormula, fData, cVar, optVar)

或者，如果我从函数内部全局分配变量，它也可以工作

#-------------------------
# DEMO FUNCTION
#-------------------------
# f = formula (string)
# fData = data.frame
# cVar = choice variable (string)
# optVar = alternative variable (string)
##########################
mlogitSum_rev <- function(f, fData,  cVar="choice", optVar="option"){
  fData<<-fData
  cVar<<-cVar
  optVar<<-optVar
  #return(head(lcmData))
  library(mlogit)
  #mi serve per poi estrarre model.matrix(r2), per il resto sarebbe ridondante
  r2 <- mlogit(as.formula(f), shape = "long", data = fData, alt.var=optVar, choice = cVar)
  return(summary(r2))
}

mydata <- read.csv(url("http://www.ats.ucla.edu/stat/r/dae/mlogit.csv"))
attach(mydata)
library(mlogit)
mydata$brand<-as.factor(mydata$brand)
mlData<-mlogit.data(mydata, varying=NULL, choice="brand", shape="wide")

myFormula <-"brand~1|female+age"
var1 <- "brand"
var2 <- "alt"
mlogitSum_rev(myFormula, mlData, var1, var2)

有关如何避免全局分配变量的任何想法？

Answer 1

tl; dr这似乎是mlogit中的一个错误，您可以自行修复（见下文）或要求维护者修复。

在mlogit内部，该函数尝试评估数据如下：

nframe <- length(sys.calls()) ## line 11
... 
data <- eval(mldata, sys.frame(which = nframe))  ## line 44

对于R的范围结构，这是一个中等复杂的混乱 - 它试图在当前帧之上的帧中评估mldata，如果有人做了一些棘手的事情（但非常合理！），就像在函数内调用mlogit一样，它将会失败。

我通过运行fix(mlogit)解决了问题（有点！），它将把你转储到编辑器中并允许你修改函数。我把第44行换成了

data <- eval(mldata, parent.frame())

之后代码似乎有效。

如果这适合你，你可以（1）fix() mlogit每次你需要使用它：（2）下载源（.tar.gz）包的副本，修改它，并安装它;或者（3）[最好！]联系软件包维护者，让他们知道问题，并要求他们发布修补版本......

PS取决于您的一般数据分析协议，您可能想要摆脱使用attach的习惯：Why is it not advisable to use attach() in R, and what should I use instead?