如果在返修的过程中,审稿人提出一个问题:你的回归模型是否靠谱?请画出回归诊断图。 下一步,建立一个包含两个自变量的线性回归模型,并且查看模型情况,代码如下: mymodel <- lm(mpg ~ hp + wt, data = mtcars)
summary(mymodel)
par(mfrow = c(2, 2)) # 改变图片输出格局
plot(mymodel)
par(mfrow = c(1, 1)) # 将图片格局恢复为默认形式
如果希望只画出其中的一张图,而并不是全部画出,可以这么做,代码如下:
plot(mymodel, which = 1)
plot(mymodel, which = 2)
plot(mymodel, which = 3)
plot(mymodel, which = 5)
也可以选择其他R包来解决上面的问题,安装并且载入R包: install.packages("performance")
library(performance)
同样,也可以使用下方代码将上述图片一张一张的画出来:
plot(check_model(mymodel))[[1]]
plot(check_model(mymodel))[[2]]
plot(check_model(mymodel))[[3]]
plot(check_model(mymodel))[[4]]
plot(check_model(mymodel))[[5]]
plot(check_model(mymodel))[[6]]
好啦,今天的内容就到这里。如果有帮助,记得分享给需要的人!https://cran.r-project.org/web/packages/performance/index.html
公众号核心成员担任SCI杂志Associate Editor!▌课程相关咨询可添加R师妹微信: kefu_rstats