nlme 和 lme4 忽略平方项答案

【问题标题】：nlme and lme4 Ignoring squared termsnlme 和 lme4 忽略平方项
【发布时间】：2013-09-02 21:53:14
【问题描述】：

我正在尝试建立一个标准的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 的等价物。

平方项只是没有出现在摘要（模型）中，而且我知道它们被忽略了，因为我用平方项创建了单独的向量并将它们传递进去，并且估计值不同。

有人有什么建议吗？

【问题讨论】：

标签： r statistics regression lme4 mixed-models

【解决方案1】：

在公式中，术语^2用于创建变量的交互。例如，(x+y+z)^2创建主要效果和与两个变量的所有可能的交互，即987654323 @。因此，公式内的987654324 x。

此外，还使用* x*y 987654328 @。因此，公式内的987654329 x。要在公式内创建平方值，您必须使用I，i.，I(x^2)或I(x*x)。

lme(asinh(gallons) ~ asinh(rprc) + asinh(rexp) +
      I(asinh(rexp)^2) + I(asinh(rprc)^2) + asinh(rprc)*asinh(rexp), 
    random=~1|cuid, data = data)

【讨论】：

谢谢！但我确实努力创造一个互动。例如，我希望价格与自己互动，因为价格变化的边际影响取决于价格。 span>
@ User2611609您不能使用Squared术语吗？ span>
是的！谢谢你。我误解了。我没有意识到它的工作原理。我现在可以实际简化我的整体式代码。谢谢！