我想用时间来模拟我的响应变量，我用了一个二次项。
**通常只有一个线性项时，我使用emmeans包的函数emtrends来获得我的因子'Side'的两个水平的两个趋势之间的差异，如何找到趋势的差异时，有一个二次项“I（Age^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
1条答案
按热度按时间zbwhf8kr1#
在没有可重复数据的情况下，很难正确回答，但emtrends可以处理多项式项。对于您的模型，不要使用Age + I（Age）^2，而是尝试：
然后
会给予你线性和二次趋势的成对对比。
你的方差分析结果表明，由于Side，趋势没有差异，所以我可能不会在这个特定情况下做事后检验，但这取决于你。