科研有捷径,输入代码,一键获取科研成果!就是这么省事,来具体看下有多方便!
搜索http://985.so/a9kb查看全部代码(目前共计50+持续新增中),也可以点击右侧【目录】,可以看到更多有趣的代码~
接下来学习下table1包生成该表格的具体用法:
table1包使用时,其中要包含在表中的变量由“+”分隔,“分层”变量(创建列)出现在条件符号“|”的右侧。
#安装并加载table1包进行绘图
install.packages(“table1”)
library(table1)
library(boot)
melanoma2 <- melanoma
# Factor the basic variables that
# we’re interested in
melanoma2$status <-
factor(melanoma2$status,
levels=c(2,1,3),
labels=c(“Alive”, # Reference
“Melanoma death”,
“Non-melanoma death”))
table1(~ factor(sex) + age + factor(ulcer) + thickness | status, data=melanoma2)
这样一个简单的统计表格就出来了,但是因为我们没有设置变量名和类别标签,所以看起来表格表达的信息可能不是很直观。接下来对这里进行改进,我们可以为分类变量(sex和ulcer)创建带有描述性标签的因子,按照我们想要的方式标记每个变量,并为连续变量(age和thickness)指定单位,如下所示:
melanoma2$sex <-
factor(melanoma2$sex, levels=c(1,0),
labels=c(“Male”,
“Female”))
melanoma2$ulcer <-
factor(melanoma2$ulcer, levels=c(0,1),
labels=c(“Absent”,
“Present”))
label(melanoma2$sex) <- “sex”
label(melanoma2$age) <- “Age”
label(melanoma2$ulcer) <- “Ulceration”
label(melanoma2$thickness) <- “Thickness”
units(melanoma2$age) <- “years”
units(melanoma2$thickness) <- “mm”
table1(~ sex + age + ulcer + thickness | status, data=melanoma2, overall=”Total”)
考文献及网址:
[1]https://cran.r-project.org/web/packages/table1/vignettes/table1-examples.html