为什么log（）语句导致lmList失败？

Question

在处理这个问题时：Interference between lme4 and nlme packages OP和我确定将log()中的响应向量包含在lme4::lmList的调用中会导致lmList失败。

如，

my.lmList.model <- lmList( log(response) ~ log(predictor) | group, mydata)

给出回复

Error in eval(expr, envir, enclos) : object 'response' not found

但如果您首先记录响应和预测器，它就可以正常工作，即

mydata$log.response <- log(mydata$response)
mydata$log.predictor <- log(mydata$predictor)
my.new.lmList.model <- lmList( log.response. ~ log.predictor | group, mydata)

有人可以解释为什么在第一个例子中将响应包装在log（）中并不起作用吗？在我看来（至少是直觉）它应该。

Answer 1

我不确定你想要什么答案。 “为什么”是因为lmList的设计是如何运作的。您已经知道了解决方法。你想知道如何重写lmList以使它在这种情况下不会失败吗？这不是一个简单的功能。在仔细阅读文档并查看函数以查看这是否是故意行为后，最好编写包作者。