1

この質問に取り組む際に: 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() にラップしてもうまくいかない理由を誰かが説明できますか? 私には（少なくとも直感的には）そうすべきだと思われます。

4

1 に答える 1

0

どんな答えが欲しいのかわかりません。「理由」は、それがlmListが機能するように設計されているためです。あなたはすでに回避策を知っています。このような状況で失敗しないようにlmListを書き直す方法を知りたいですか？単純な関数ではありません。ドキュメントを注意深く読み、関数を調べてこれが意図的な動作であるかどうかを確認した後、パッケージの作成者を作成することをお勧めします。

于 2011-12-13T15:01:09.680 に答える