使用 sjplot 绘制曲线图
Curvilinear plots with sjplot
我有一个时间和时间二次方 (time2) 的混合模型(使用 nlme)。我怎样才能让 sjp.int (of 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)
请注意,交互项的预测变量的顺序很重要,并且在公式中应该相同(即在上述情况下,c161sex
首先出现,然后是 e17age
)。
我有一个时间和时间二次方 (time2) 的混合模型(使用 nlme)。我怎样才能让 sjp.int (of 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)
请注意,交互项的预测变量的顺序很重要,并且在公式中应该相同(即在上述情况下,c161sex
首先出现,然后是 e17age
)。