NLME和lme4忽略平方項

Question

我想建立一個標準的超越對需求的功能，那就是：

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) 

而且我試過11聚物爲等同。

的平方項就是不以摘要（模型）顯示了，我知道，因爲我已經創建與平方項不同的載體，並通過那些和估計是不同的他們被忽略。

有人有什麼建議嗎？

Answer 1

在式中，術語^2用於創建的變量的相互作用。例如，(x+y+z)^2創建主效應和與兩個變量的所有可能交互，即x + y + z + x:y + x:z + y:z。因此，公式中的x^2與x相同。

此外，*也用於創建交互，例如，x*y創建x + y + x:y。因此，公式中的x*x與x相同。

要在公式中創建平方值，必須使用函數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)