平时一遍又一遍的分析数据,还要一遍一遍的誊分析结果然后画表格。烦不烦人?

有没有想过让R自动生成所有的结果表格?可以直接发表的那种?今天就给大家写写如何做,提醒一下:关注我,你不会后悔的。

回归结果表格自动生成

先给大家做一个很简单回归分析,体验一下R自动输出结果的强大,首先载入相关包和数据

library(sjPlot) library(sjmisc) library(sjlabelled) data("efc") efc <- as_factor(efc, c161sex, c172code)

上面的代码应该不用解释吧,我用的是自带的欧洲家庭调查数据,你也可以试着用自己的数据,接下来拟合模型

m1 <- lm(barthtot ~ c160age + c12hour + c161sex + c172code, data = efc)

我用c160age + c12hour + c161sex + c172code这些变量来预测生活自理能力barthtot的回归模型。

自动结果输出只需要tab_model方法就可以

tab_model(m1)


可以看到,输出直接是基本上可以发表的结果表格,不用费心编辑拷贝啦,简直完美。输出是一个html的文件大家可以直接转为doc做相关编辑,很方便的。

多个回归结果自动展示

在上面的例子中我们只拟合了一个模型,我们也可以拟合2个模型,模型结果同时在一个表格中展示,比如我们再拟合一个m2

m2 <- lm(neg_c_7 ~ c160age + c12hour + c161sex + e17age, data = efc)

同时展示2个模型结果:

tab_model(m1, m2)


可以看到在输出的结果表格中2个模型并排放置,对同一个变量,每个模型都输出了系数置信区间和p值

自定义输出

有时候我们根据杂志要求,有可能并不需要提供某个参数,此时可以用tab_model()自定义哪些参数输出,哪些隐藏。比如说我们要输出系数的标准差,就可以将show.std设置为True。其他的可以设置的参数见下图:


小结

今天给大家写了如何用R生成了一个可以发表的结果表格,这些表格都是可以存为word文档再手工编辑的,真的非常好用,同时像我今天介绍的tab_model方法还有很多功能,比如输出的合并,表格顺序的调整,自动标记p值等等,关注我,以后慢慢给大家写。感谢大家耐心看完。发表这些东西的主要目的就是督促自己,希望大家关注评论指出不足,一起进步。内容我都会写的很细,用到的数据集也会在原文中给出链接,你只要按照文章中的代码自己也可以做出一样的结果,一个目的就是零基础也能懂,因为自己就是什么基础没有从零学Python和R的,加油。

(站外链接发不了,请关注后私信回复“数据链接”获取本头条号所有使用数据)

往期精彩:

R数据分析:如何用R做验证性因子分析及画图,实例操练