私は、次のような標準のトランスログ要求機能を構築しようとしています。
lnY = lnP + lnZ + lnY*lnZ + lnY^2 + lnZ^2
ここでY、= 需要、P= 価格、およびZ= 収入です。
ただし、nlme または lme4 に 2 乗項を含めると、無視されます。私はもう試した：
model <- lme(asinh(gallons) ~ asinh(rprc) + asinh(rexp) + asinh(rexp)*asinh(rexp) + asinh(rprc)*asinh(rprc) + asinh(rprc)*asinh(rexp), random=~1|cuid, data = data)
そして私は試しました：
model <- lme(asinh(gallons) ~ asinh(rprc) + asinh(rexp) + asinh(rexp)^2 + asinh(rprc)^2 + asinh(rprc)*asinh(rexp), random=~1|cuid, data = data)
そして、lmerに相当するものを試しました。
二乗項は要約（モデル）に表示されません。二乗項を使用して別のベクトルを作成し、それらを渡したところ、推定値が異なるため、それらが無視されていることがわかりました。
誰かアドバイスはありますか？