审稿人:你的回归模型靠谱吗?请提供回归诊断图!谢谢!

学术   科学   2024-03-04 11:31   浙江  
如果在返修的过程中,审稿人提出一个问题:你的回归模型是否靠谱?请画出回归诊断图。 

那该怎么办?

这篇文章提供两种代码供大家选择。

首先查看数据集,代码如下:

summary(mtcars)

下一步,建立一个包含两个自变量的线性回归模型,并且查看模型情况,代码如下: 

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)
R包get!

一行代码就可以画出6张回归诊断图,代码如下:

check_model(mymodel)

同样,也可以使用下方代码将上述图片一张一张的画出来:

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

公众号的线上课程
1. 《R语言和统计新手课程》
2. 《回归:从入门到进阶》

统计咨询
《服务介绍和经典合作案例》

公众号核心成员的成果发表
《SCI医学1区影响因子9分论文》

公众号核心成员担任SCI杂志Associate Editor!
《JAD杂志Associate editor》
《Frontiers in Neuroscience, Frontiers in Neurology and Frontiers in Psychiatry杂志的神经退行性病变板块》


▌本文由R语言和统计首发
▌课程相关咨询可添加R师妹微信: kefu_rstats
▌编辑:June
▌邮箱:contact@rstats.cn
▌网站:www.rstats.cn
我们致力于让R语言和统计变得简单!




R语言和统计
我们定期更新与R有关的内容,比如R编程基础,作图,实用R包的解读,统计学基础知识,前沿的统计方法,机器学习等等。
 最新文章