我正在尝试建立一个标准的translog需求函数,即:
lnY = lnP + lnZ + lnY*lnZ + lnY^2 + lnZ^2
其中Y
= 需求、P
= 价格和Z
= 收入。
但是,当我在 nlme 或 lme4 中包含平方项时,它们会忽略它们。我试过了:
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 的等价物。
平方项只是没有出现在摘要(模型)中,而且我知道它们被忽略了,因为我已经用平方项创建了单独的向量并将它们传递进去并且估计值不同。
有人有什么建议吗?