Почему оператор log () приводит к сбою lmList?

Question

При работе над этим вопросом: Помехи между пакетами lme4 и nlme мы с 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, чтобы он не работал в таких условиях?Это не простая функция.Лучше всего написать авторам пакета после тщательного прочтения документации и просмотра функции, чтобы выяснить, является ли это преднамеренным поведением.