R 带sjplot的曲线图

我有一个混合模型（使用nlme），带有时间和时间二次型（time2）。如何让sjp.int（sjplot的）要求将其绘制成曲线

   model = lme(MR ~ Age* time + sex*time + ed* time+ time* HMD+time2*HMD, random=~1|subject, na.action="na.omit",data=mydata);summary(model)

我把sjplot称为：

  sjp.int(model, type="eff", swap.pred=T, mdrt.values="quart", show.ci=T)

---

二次项应具有相同的名称，即应使用

poly（）

或

i（）^2

。请参见以下示例：

library(sjPlot)
library(sjmisc)
data(efc)
fit <- lm(tot_sc_e ~ c161sex * e17age + c161sex * I(e17age^2), data = efc)
sjp.int(fit)

库（sjPlot）
图书馆（sjmisc）
数据（efc）
拟合你的二次项应该有相同的名称，即你应该使用
poly（）
或
i（）^2
。请参见以下示例：

library(sjPlot)
library(sjmisc)
data(efc)
fit <- lm(tot_sc_e ~ c161sex * e17age + c161sex * I(e17age^2), data = efc)
sjp.int(fit)

库（sjPlot）
图书馆（sjmisc）
数据（efc）
嗨，丹尼尔，谢谢你的帮助，这真是个好办法！还有一个问题，是否有办法绘制三方交互（例如，HMD*时间*组），然后绘制一个效果HMD*时间在我的两个组中分割的图？我尝试了facet.grid=TRUE，但它给了我奇怪的图（我尝试将因子转换为数字，但没有帮助），三方交互尚未实现，请参阅。还需要一些概念上的思考。我大牛，谢谢你的帮助，这确实是个窍门！还有一个问题，是否有办法绘制三方交互（例如，HMD*时间*组），然后绘制一个效果HMD*时间在我的两个组中分割的图？我尝试了facet.grid=TRUE，但它给了我奇怪的图（我尝试将因子转换为数字，但没有帮助），三方交互尚未实现，请参阅。还需要一些概念性的思考