描述统计是统计学重要的一部分内容,尤其是在医学中的应用更广。很多医学、或者做统计的SCI文章,开头就是统计的内容,一般是由三线表的方式呈现的。之前小编也只会用excel统计然后自己制作表格,但这样效率很慢。这里我们介绍一个R包---table1,可以非常简单快速的完成统计工作,并制出三线表!
安装包和示例数据
install.packages("table1")
install.packages("boot")
library(table1)
library(boot)
data <- melanoma
数据就是一般的描述数据,性别、status都是用数字表示的
图片对数据进行转化,将sex等换成字符
data$status <- factor(data$status,
levels=c(2,1,3),
labels=c("Alive", "Melanoma death", "Non-melanoma death"))
data$sex <- factor(data$sex, levels=c(1,0), labels=c("Male", "Female"))
data$ulcer <- factor(data$ulcer, levels=c(0,1),labels=c("Absent", "Present"))
图片
之后作图,分组用status,table1函数的用法具体可参考帮助函数!一般格式为~不同变量+变量+......|分类变量,data。
table1(~ factor(sex) + age + factor(ulcer) + thickness | status, data=data)
图片
如果用性别做分类变量,则效果如下。
table1(~ status + age + factor(ulcer) + thickness | factor(sex), data = data)
图片
还可以为变量加上单位或者改变变量名
units(data$age) <- "years"
units(data$thickness) <- "mm"
table1(~ sex + age + ulcer + thickness | status, data=data, overall="Total")
图片
最后得到的表格可以复制,直接粘贴到PPT或者AI中进行编辑整理!
下节我们将继续说说table1制作统计表并添加统计分析,以及对表格的各种样式调节!
网友评论