R语言 调节变量回归分析的交互作用图

xyhw6mcr  于 2023-01-18  发布在  其他
关注(0)|答案(1)|浏览(316)

我正在为我的回归分析获取交互图。我已经测试过回归了。很好。代码如下:

plot(data$ocs, data$efs, pch=19)
abline(lm(data$efs~data$ocs))

model1 <- lm(efs ~ ocs * rws, data = data)
summary(model1)
interact_plot(model=model1, pred = ocs, modx = rws, data = data)

现在,除了interact_plot之外,一切正常。我返回以下错误消息:

Error in element_line(colour = text_and_line_color, linewidth = 0.5, linetype = 1,  : 
  unused argument (linewidth = 0.5)

这到底是什么意思?
interaction_plot是交互作用的函数!

rws = c(3.8, 3.8, 4.4, 3.2, 4, 3.8, 4.2, 4.6, 4.4, 3, 3.6, 3.4, 3.8, 4, 3.4, 
4, 2.2, 4, 3.4, 4, 2.8, 3, 3, 4.2, 3, 3.6, 3.8, 3.8, 3, 3.8, 
3.4, 3.2, 3.4, 3.4, 3.8, 4, 3.2, 3, 4.6, 4, 3.4, 3, 2.8, 3.8, 
3.8, 3.6, 2.8, 3.8, 4.6, 3, 3.8, 3.8, 3.4, 3.8, 2.4, 2.4, 4.2, 
3.6, 2.6, 3.6, 2.4, 3.2, 4, 4.2, 3.4, 3.2, 3.8, 4.8, 3.6, 4, 
3, 3.4, 4.2, 4.4, 3.8, 4.2, 4, 3.6, 4.4, 4, 3, 4.6, 4.4, 4, 3.4, 
4, 3.8, 4.8, 3.4, 3.4, 3.8, 3.8, 3.6, 3.6, 4.6, 3.6, 3.6, 2.6, 
2.8, 4.6, 4, 3.4, 3.8, 3.8, 3.6, 3, 3.2, 2.8, 3.6, 3.6, 3.4, 
3.6, 2.2, 4.2, 3.6, 3.8, 3.4, 3.8, 4.4, 2.4, 2, 3.4, 4.2, 3.2, 
3.2, 3.2, 3, 3.6, 3.6, 2.8, 3.6, 2.8, 3.6, 4.4, 3, 3.6, 3.2, 
4.8, 3.4, 3.8, 4.2, 3.2, 3.6, 3.8, 3.6, 3, 4, 3.8, 3.8, 4.2, 
3.6, 4.4, 3.6, 3.4, 3.8, 4.4, 3.4, 3.4, 4, 4.2, 3.8, 3.6, 4, 
4.6, 3.2, 4, 3.4, 3, 3.2, 3.8, 3.4, 3.4, 3, 4.8, 2.8, 3.6, 3, 
4, 3.4, 2.8, 3.2, 3.4, 3.8, 3.2, 4.2, 3.2, 3.6, 3.4, 4.2, 3.8, 
3.6, 4, 3.6, 3.4, 2.4, 4, 3.4, 3.6, 4, 3.4, 3.6, 3, 3.6, 3.6, 
2.4, 3.6, 1, 3.6, 4.2, 3.2, 2.6, 3.8, 3.6, 3.6, 3.2, 3.2, 3.2, 
3.6, 3, 4, 4, 3.2, 4.4, 3.4, 4.4, 3, 3.8, 3.6, 3.2, 3, 4.6, 4.8, 
3, 4.2, 3, 3.4, 4, 3.2, 3.2, 4.2, 4.2, 3, 2.6, 2.8, 2.2, 3.8, 
2.8, 3.8, 1.4, 3.8, 4.6, 3.8, 3.6, 4.2, 3.6, 3.4, 2.6, 3.8, 3.8, 
3.8, 2.8, 3.8, 4.8, 3.6, 3.6, 3.6, 4, 3.4, 3.6, 3.2, 2.6, 3.2, 
4, 3.4, 3.2, 4, 3, 4, 3.6, 3.4, 2.8, 3.2, 3, 3.4, 4.6, 3.4, 3.4, 
3.4, 4.2, 2.8, 2.4)

efs = c(3.5, 3.75, 4, 3.75, 4.25, 3.5, 4.5, 2.5, 4.25, 2, 2.5, 3, 4.5, 
2.25, 3.5, 3.5, 1.75, 4, 3, 2.75, 3.75, 3.5, 3.5, 2.25, 3, 2.25, 
3.25, 2.5, 3.75, 2, 3.5, 2.75, 4.25, 3.5, 4, 2.5, 3, 4.5, 2.5, 
3.75, 1.25, 3.75, 2.25, 3.25, 3, 4.75, 3.25, 5, 3.75, 3.5, 2.75, 
3.75, 3.75, 3.25, 4, 2.75, 4.75, 2, 3.5, 3.5, 2.75, 1.75, 1.75, 
3, 3, 3.5, 3.5, 2.25, 2.5, 2.25, 3.25, 3.75, 3, 4.25, 3.75, 4, 
4.5, 3.75, 3.75, 4.25, 3.5, 3, 3, 2.75, 4, 3, 4.25, 3.25, 3.25, 
3, 3.25, 2, 3.75, 2.5, 3.5, 4, 4, 3.25, 4.25, 4.75, 4, 4, 4.75, 
3.75, 4, 4, 4, 3.5, 3.25, 4.25, 4.5, 3.75, 4, 4.5, 3.5, 2.75, 
4, 3, 4, 3, 4, 4.5, 2.75, 4.25, 3.5, 4, 3.75, 3.75, 3.75, 4.25, 
5, 4.75, 3.25, 2.5, 2.75, 3.5, 4, 2.75, 3.25, 3, 3, 3.75, 4, 
3.75, 2.5, 3.75, 3.5, 4, 3.5, 2.5, 3.5, 2.75, 3.25, 2.75, 3, 
3.5, 3, 4.25, 3, 3.25, 2.75, 4.75, 2.25, 4, 3.75, 3.75, 4.75, 
3.25, 4, 3.75, 3.5, 3.75, 4.25, 4.25, 4.25, 3.75, 4, 3, 2.5, 
4.25, 3.25, 4.25, 3.75, 3.75, 3.75, 3.5, 3.75, 4, 3.25, 2.75, 
3.25, 4.25, 3, 3.25, 3, 3.75, 2.5, 2.75, 3.5, 3.5, 3.75, 4.75, 
3.5, 3.5, 4, 3.5, 3.5, 3, 4, 3.5, 3.75, 4, 3.25, 3.5, 3.25, 4.25, 
4.5, 3, 3, 4.75, 3, 3, 2.75, 3, 3.25, 3.25, 4.25, 4.25, 3.5, 
4, 2.25, 4.5, 3.75, 4, 4, 3, 3.75, 2.5, 3.5, 3.5, 3.5, 4.25, 
3.75, 4.25, 3.5, 2.75, 4, 3.25, 2.5, 4, 2.75, 2.5, 3, 2.5, 2.75, 
4.5, 3.25, 3.25, 3.75, 3.75, 3.5, 2.5, 3.5, 2.75, 3.5, 3.5, 3.75, 
4.5, 3, 3.25, 3.75, 4.5, 3.5, 3, 3.75, 4, 3.75, 4, 3.75, 3.75, 
3, 3.75, 3.25, 4.5, 2.25, 3.25, 4, 3.75, 3.5, 2.75, 3.75)

ocs = c(2.8, 3.2, 3.4, 2, 2.6, 2.4, 2.6, 2.2, 3.6, 2.4, 2.2, 3.4, 2.8, 
1.4, 2.8, 2, 1.8, 2.8, 2.6, 2.6, 2.8, 2.6, 3.6, 1.4, 2.8, 1.8, 
2.6, 1.8, 3, 2, 2.6, 2.2, 3.6, 2.8, 3, 1.2, 2.8, 2.8, 2.8, 2.8, 
2.6, 2.6, 1.8, 2, 2.2, 4.4, 3.2, 4.4, 3.8, 3.4, 2.8, 2.6, 2.6, 
2.6, 2.8, 1.8, 4, 3.4, 2.6, 3.2, 2.8, 1.8, 1.2, 2, 2.4, 3.4, 
2.2, 1.6, 2.6, 1.8, 2.4, 3.2, 2.6, 2.8, 2.8, 2.4, 3.2, 2.6, 2.8, 
3.2, 3.4, 1.4, 1.6, 2, 3.2, 2.8, 3.6, 4, 2, 1.6, 2.8, 1.6, 2.6, 
2.6, 3.4, 2.8, 3.4, 2.6, 3.4, 3, 3.8, 3, 4, 3.2, 3.4, 3.2, 4.4, 
2.8, 3.2, 3.4, 4, 2, 3.8, 2.8, 2.8, 2.2, 3.4, 2.2, 3, 2.4, 3, 
3.2, 2, 1.8, 3, 2.2, 3.4, 3.6, 3.6, 3.2, 4.2, 3.6, 3, 2.6, 2.4, 
3.6, 4.2, 1.6, 3.4, 2.6, 3, 2.6, 4, 2.2, 2.8, 3.4, 1.8, 3.4, 
2, 2.6, 2.6, 2.6, 3, 1.8, 2.2, 2.6, 3, 3.4, 3.6, 1.6, 1.8, 4, 
2.4, 2.6, 2.6, 2.6, 3.6, 3, 3.8, 2.6, 3, 3.2, 4.6, 3.8, 2.8, 
2.6, 3.8, 2.2, 2.8, 3.2, 3.4, 3.4, 2.6, 1.4, 2.4, 2.8, 3.4, 2.6, 
2.8, 2.4, 2.6, 2.2, 3.2, 2.4, 2.6, 3, 2.4, 2.8, 3.2, 3, 3.6, 
4, 2.6, 3.2, 4, 2.8, 0, 3, 3.2, 3, 3.4, 3.2, 2.4, 2.2, 3, 3, 
3.6, 2.6, 2.8, 3.6, 2.4, 2.2, 2.4, 2.6, 2.8, 1.8, 3.2, 2.6, 3.4, 
3, 2.2, 4, 3, 3.6, 4.2, 1.8, 2.8, 2.4, 1.8, 3.4, 3.4, 3.8, 2.6, 
3.6, 2.4, 2, 3.4, 3.4, 0, 2.8, 2.4, 1.4, 2.8, 2.2, 2.4, 3.4, 
3.2, 1.8, 2, 3.8, 3.4, 1.8, 2, 2.6, 2.6, 3.4, 3.8, 3.6, 2.8, 
3.2, 3.4, 3.2, 2.4, 2.6, 3, 3.2, 3, 2.4, 3.4, 3.6, 3.6, 3, 3.6, 
4, 2.2, 2.8, 2.8, 3.2, 2.2, 1.6, 3.4)
3qpi33ja

3qpi33ja1#

首先,将数据放入 Dataframe 并运行模型:

data<-as.data.frame(cbind(rws,efs,ocs))
model1 <- lm(efs ~ ocs * rws, data = data)

接下来,sjPlot包在这种情况下非常有用。

library(sjPlot)
sjPlot::plot_model(model1, type='int')

您还可以使用相同的包创建整洁的表格:

sjPlot::tab_model(model1)

相关问题