问题描述
我想用置信区间绘制斜率,描述两个斜率与预测值(或边际效应)图的差异。作为说明,我使用下面的“mtcars”数据并进行以下交互模型:
data("mtcars")
main_model <- lm(mpg ~ am + wt + am*wt,mtcars)
summary(main_model)
其中,我的因变量“mpg”是一个连续变量,我的自变量是“am”(一个二元变量)和“wt”(一个连续变量)。
在建立模型后,我使用“ggpredict”来估计因变量的预测值,基于“am”的两个可能值和“wt”的条件:
predict <- ggpredict(main_model,c("wt[all]","am"),ci.lvl = 0.95)
然后,我使用“plot”绘制预测值:
plot_pred <- plot(predict,colors="bw") + labs(
x="Weight (1000 lbs)",y="Miles/(US) gallon") +
theme(legend.position="bottom")
plot_pred
图中有两个斜率:1.当“am”为0时(实线的斜率); 2.当“am”为1时(虚线的斜率)。现在,我想绘制一个类似的图形,但只有一个斜率是这两个斜率(1 和 2)之间的差值,置信区间为 95%。
是否有任何软件包可以直接生成我想要的情节,或者在调用“ggpredict”或“情节”时?
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)