存在二次项的趋势：二次趋势的差异

Question

我想按时间为我的响应变量建模。我用了一个二次项。

通常当我只有一个线性项时，我使用 emmeans 包的函数 emtrends 来获得我的因子'Side'的两个水平的两个趋势之间的差异，当有二次项时如何找到趋势的差异“我（年龄^2） emtrends 能做到吗？

> mod.quadra<-lmer(RV~Side*(Age + I(Age^2))+(1|ID),data)
> Anova(mod.quadra)
                      Chisq Df Pr(>Chisq)    
Side                2.0791  1  0.1493267    
Age                 12.0248  1  0.0005250 ***
I(Age^2)            10.9917  1  0.0009152 ***
Side:Age            1.8857  1  0.1696876    
Side:I(Age^2)       2.0099  1  0.1562804   

我通常是怎么做的（当我只有一个线性项时）

>mtrend<-emtrends(LM.fit, "laterality", var = "Age")

 Side       Age.trend    SE  df lower.CL upper.CL
 contra          1.03 0.648 159   -0.251     2.31
 ipsi            2.03 0.648 159    0.755     3.31

>pairs(mtrend)
 contrast      estimate    SE  df t.ratio p.value
 contra - ipsi    -1.01 0.891 149  -1.129  0.2606