R语言|4. 轻松绘制临床基线表Table 1 临床三线表绘制
R语言|4. 轻松绘制临床基线表Table 1 regular table regular
欢迎关注公众号
第四次考核 Jimmy 学徒考核 Linux安装软件 rnaseq上游分析-2 ascp kingfisher数据下载ena (qq.com)
https://mp.weixin.qq.com/s?__biz=Mzg2NDcxMzYwNg==&mid=2247483760&idx=1&sn=dccdb31eba73a0f1330f5e78d1de1bf7&chksm=ce646f99f913e68feb1d69c137c8ebbac24d8b79f5f4c172556af5f7b34e2f36189aaffbd2b5&token=2005279620&lang=zh_CN#rd 2021-01-29 11:00
Table 1是总结研究样本中变量组成及变量间关系的一张表。
R语言可以通过几行代码快速生成Table 1。
本期用到的是“tableone”，是我制作Table 1时最喜欢用的包，今天分享给大家，觉得有用的小伙伴希望多加点赞转发支持。
目录
一.文献中的Tabl 1
二. tableone包制作Table 1（2列）
三. tableone包制作Table 1（4列）
四. tableone包制作Table 1（5列）
五. tableone其他代码展示
一
文献中的Table 1
图片来自：https://www.jianshu.com/p/3328fd0689f2
连续和分类变量相关性检验详见之前文章：
R语言|2. 连续变量相关性比较 和 R语言|3. 分类变量检验
1-1. 变量名+总人群（两列）
研究目的为“xx疾病的预后因素分析”时常用，即不涉及比较。
这种表只需把握好连续变量为正态还是非正态分布即可。正态分布写均值±标准差，非正态分布写中位数+上下四分位数。
分类变量写数目和百分数。
例：研究某类型乳腺癌的预后，所以Table 1两列
参考文献：Lewis GD，et al. Prognosis of lymphotropic invasive micropapillary breast carcinoma analyzed by using data from the National Cancer Database. Cancer Commun (Lond). 2019 Oct 21;39(1):60.
1-2：变量+某变量分列+p值
研究目的为 “xx变量对预后的比较分析” 时常用。
因涉及比较，所以在Table 1常将用该变量将全部人群分为几类，观察其他变量在这几个分类中是否有差异（p值）。
这种表就涉及到了连续变量和分类变量的相关性检验，这两种检验在R语言|2. 连续变量相关性比较 和 R语言|3. 分类变量检验已有介绍。
下面这篇文章的目的是研究两种病理类型的乳腺癌临床病理特点的比较。Table 1中作者总结了出这两个病理类型在肿瘤位置、组织分级等变量中有何差异
参考文献：Yu JI,et al. Differences in prognostic factors and patterns of failure between invasive micropapillary carcinoma and invasive ductal carcinoma of the breast: matched case-control study. Breast. 2010 Jun;19(3):231-7.
1-3. 变量名+总数+某变量分组+p值
1-1和1-2的综合。可以向读者展示更多细节。通常5列。
文章目的是两种病理类型的乳腺癌临床病理特点的比较。
参考文献：Hashmi AA, et al. Clinicopathologic features of invasive metaplastic and micropapillary breast carcinoma: comparison with invasive ductal carcinoma of breast. BMC Res Notes. 2018 Jul 31;11(1):531.
二
tableone包制作Table 1（2列）
以下数据来自SEER数据库的5000例乳腺癌病例，公众号回复：20210126获得
总结：过程其实就是输入3个条件【指定需要汇报的变量;指定这些变量中的分类变量；指定哪些连续变量是非正态分布】-----构建函数-----输出表格
1
安装R包、数据放入工作目录
#1.加载R包，电脑里若无tableone，去掉下面的#install.packages("tableone")library(tableone)#2.清理运行环境rm(list = ls())#3.读入数据aa<- read.csv('20210126.csv')
2
查看数据构成
#4.查看数据前6行head(aa)
源数据尽量用英文单词或缩写而非0.1.2.3.4代替。因软件默认识别为连续变量，需要转换，很麻烦。用英文在整理表格时也方便，省时省力。
#5.查看数据数据性质str(aa)
数据有5000病例；16变量，2连续变量（年龄和时间）；14分类变量及它们亚变量的个数和名字。
#6.提取变量的名字names(aa)
age=年龄，age70=以70岁为界分组的年龄，race=种族，marry=婚姻，t=t分期，n=n分期，tnm=tnm分期，er=雌激素受体，pr=孕激素受体，her2，g=组织分级，sur=手术，rt=放疗，che=化疗，status=死亡否，time=时间
3
连续变量正态性检验
shapiro.test(aa$age)#p＞0.05才符合正态分布shapiro.test(aa$time)#p＞0.05才符合正态分布
age和time SW检验p<0.05, 不符合正态性分布，Table1应汇报中位数+四分位数
前面的代码是只是查看数据性质，从第4步开始正式写Table 1代码。
4
输入Table 1的条件
#条件1# myVars的（）中输入想要在Table 1出现的变量【英文引号以英文逗号隔开】
myVars <- c("age","age70"，"race", "marry","t", "n", "tnm", "er", "pr","g","her2", "sur","rt","che","time","status")#16个变量
#条件2# catVars的（）内指明上述中哪些是分类变量
catVars <- c("age70"，"race", "marry","t", "n", "tnm", "er", "pr","g","her2", "sur","rt","che",,"status")#14个分类变量
#条件3# tableone包默认输出均数+标准差，所以nonvar的（）输入那些连续变量是非正态分布的（输出中位数+四分位数）
nonvar <- c("time","age") # 指定哪些变量是非正态分布变量5
构建Table 函数
table<- CreateTableOne(vars = myVars, #条件1factorVars = catVars, #条件2data = aa, #源数据addOverall = TRUE) #增加overall列
6
输出结果
table1 <- print(table, #构建的table函数（包括条件1.2）showAllLevels=TRUE, #显示所有变量nonnormal = nonvar) #条件3
age和time均汇报为中位数+四分位数，分类变量均汇报为数目加百分数
7
保存至Excel
write.csv(table1, file = "table1.csv")三
tableone包制作Table 1（4列）
4列的Table 1其实就是在构建table 1时加入了一个变量分列条件而已。
需要注意的是此时应先验证连续变量和分类是哪种检验方法。
详见：R语言|2. 连续变量相关性比较 和R语言|3. 分类变量检验
前面操作1-4不变
library(tableone) #1.加载R包rm(list = ls()) #2.清理运行环境aa<- read.csv('20210126.csv') #3.读入数据head(aa) #4.查看数据前6行str(aa) #5.查看数据数据性质names(aa) #6.提取变量的名字
#条件1不变# myVars的（）中输入想要在Table 1出现的变量
myVars <- c("age","age70"，"race","marry","t","n", "tnm","er","pr","g","her2", "sur","rt","che","time","status") #16变量
#条件2不变# catVars的（）内指明上述中哪些是分类变量
catVars <- c("age70"，"race", "marry","t", "n", "tnm", "er", "pr","g","her2", "sur","rt","che",,"status")
5
构建Table函数
#构建table 函数，加入条件4，strata = " " 。英文引号内填入需要分列的变量，例如本研究想探索放疗对预后的影响，则为strata = "rt"
table <- CreateTableOne(vars = myVars, #条件1factorVars = catVars,#条件2strata = "rt", #条件4data = aa, #原始数据);table
#条件3不变#. 指定非正态分布连续变量变量
nonvar <- c("time","age") #条件5新加入# 假如有T<5变量应使用Fisher精确检验，本文数量大，无需Fisher精确检验
exactvars <- c("a", "b")#附加细节条件#：
catDigits = 2, contDigits = 3, pDigits = 4,修改连续变量小数位数为2位,分类变量百分比位数为3位,调整小数位数为4位；
6
输出结果
把构建的table+条件1-5+附加细节条件放入print（）函数
table1<- print(table, #构建的table函数（带条件1.2.3）nonnormal = nonvar,#条件4#exact = exactvars,#条件5catDigits = 2,contDigits = 3,pDigits = 4, #附加条件showAllLevels=TRUE, #显示所有变量quote = FALSE, # 不显示引号noSpaces = TRUE, # #删除用于对齐的空格printToggle = TRUE) #展示输出结果
7
保存为Excel
write.csv(table1, file = "table1.csv")四
tableone包制作Table 1（5列）
制作4列的所有条件不变（1-5+附加条件），在构建table函数时的代码增加一个addOverall = TRUE
5
构建Table函数
table <- CreateTableOne(vars = myVars, #条件1factorVars = catVars,#条件2strata = "rt", #条件4data = aa, #原始数据addOverall = TRUE);table#条件6加入overall
6
输出Table 1（代码不变）
table1<- print(table, #构建的table函数（带条件1.2.3）nonnormal = nonvar,#条件4#exact = exactvars,#条件5catDigits = 2,contDigits = 3,pDigits = 4, #附加条件showAllLevels=TRUE, #显示所有变量quote = FALSE, # 不显示引号noSpaces = TRUE, # #删除用于对齐的空格printToggle = TRUE) #展示输出结果
函数CreateTableOne()所以细节条件汇总，即第5步。
5
构建Table函数
table <- CreateTableOne(vars, #指定哪些变量是Table 1需要汇总的变量strata,#指定进行分类的变量，不写则只出Overall列data,# 变量的数据集名称factorVars,#指定哪些变量为分类变量，指定的变量应是vars参数中的变量includeNA = FALSE,#为TRUE则将缺失值作为因子处理，仅对分类变量有效test = TRUE,#默认为TRUE，当有2个或多个组时，自动进行组间比较testApprox = chisq.test,#默认卡方检验argsApprox = list(correct = TRUE),# 进行连续校正的chisq.testtestExact = fisher.test,#进行fisher精确检验argsExact = list(workspace = 2 * 10^5),# 指定fisher.test分配的内存空间testNormal = oneway.test, # 连续变量为正态分布进行的检验，默认为oneway.test，两组时相当于t检验argsNormal = list(var.equal = TRUE), # 假设为等方差分析testNonNormal = kruskal.test,# 默认为Kruskal-Wallis秩和检验argsNonNormal = list(NULL),#传递给testNonNormal中指定的函数的参数的命名列表smd = TRUE,#如果为TRUE（如默认值）并且有两个以上的组，则将计算所有成对比较的标准化均值差。addOverall = FALSE#仅在分组时使用）将整个列添加到表中。Smd和p值计算仅使用分层的列阵进行。);table
#print（）函数细节汇总，即第6步
6
输出Table 1
table1<- print(x, #CreateTableOne()的 <- 前的名字(x=table)catDigits = 1, #连续变量小数位1位contDigits = 2,#分类变量保留2位pDigits = 3,#p值保留3位quote = FALSE,#默认值为FALSE。如果为TRUE#，则包括行名和列名在内的所有内容都用引号引起来，以便您可以轻松地将其复制到Excel。missing = FALSE,#是否显示丢失的数据信息explain = TRUE,#显示百分比时是否在变量名称中添加解释，即（％）添加到变量名称中。printToggle = TRUE,#如果为FALSE，则不输出test = TRUE,#是否显示p值。默认为TRUE。smd = FALSE,#是否显示标准化均值差异。默认为FALSE。如果存在多个对比，则显示所有可能的标准化均值差的平均值。noSpaces = FALSE,#是否删除为对齐而添加的空格。padColnames = FALSE,#是否用空格填充列名以居中对齐。默认值为FALSE。如果noSpaces = TRUE，则不进行varLabels = FALSE,#是否用从labelled :: var_label（）函数获得的变量标签替换变量名。format = c("fp", "f", "p", "pf")[1],#默认值为“ fp”频率（百分比）。您也可以选择仅“ f”频率，“仅p”百分比和“ pf”百分比（频率）。showAllLevels = FALSE,#是否显示所有级别。cramVars = NULL,#字符向量，用于指定两个级别的分类变量，对于这两个级别的变量，应在一行中显示两个级别。dropEqual = FALSE,#是否删除“ =第二级名称”描述，指示为两级分类变量显示哪个级别。exact = NULL,#字符向量，用于指定p值应为精确测试值的变量。nonnormal = NULL,#字符向量，用于指定p值应为非参数检验的变量的变量。minMax = FALSE#对于非正态变量，是否使用[min，max]而不是[p25，p75]。默认值为FALSE。)
R语言|4. 轻松绘制临床基线表Table 1 临床三线表绘制相关推荐
- R语言使用table1包绘制（生成）三线表、使用单变量分列构建三线表、编写自定义三线表结构（将因子变量细粒度化重新构建三线图）、为不同的变量显示不同的统计信息
R语言使用table1包绘制(生成)三线表.使用单变量分列构建三线表.编写自定义三线表结构(将因子变量细粒度化重新构建三线图).为不同的变量显示不同的统计信息 目录
- 三线表是什么？R语言使用table1包绘制（生成）三线表、使用单变量分列构建三线表、编写自定义三线表结构（将因子变量细粒度化重新构建三线图）、编写自定义函数在三线表中添加p值
三线表是什么?R语言使用table1包绘制(生成)三线表.使用单变量分列构建三线表.编写自定义三线表结构(将因子变量细粒度化重新构建三线图).编写自定义函数在三线表中添加p值 目录
- R语言广义线性模型Logistic回归模型亚组分析及森林图绘制
R语言广义线性模型Logistic回归模型亚组分析及森林图绘制 #Logistic回归案例 6 亚组分析森林图 library(forestplot) rs_forest <- read.csv ...
- R语言生存分析寿命表（life table）实战案例：比较两种药物治疗感染患者的生存时间
R语言生存分析寿命表(life table)实战案例:比较两种药物治疗感染患者的生存时间 目录
- R语言使用epiDisplay包的roc.from.table函数可视化临床诊断表格数据对应的ROC曲线并输出新的诊断表、设置cex参数指定AUC值及其文本标签字体的大小
R语言使用epiDisplay包的roc.from.table函数可视化临床诊断表格数据对应的ROC曲线并输出新的诊断表(diagnostic table).输出灵敏度.1-特异度.AUC值等.设置c ...
- R语言广义加性模型GAMs分析温度、臭氧环境数据绘制偏回归图与偏残差图
最近我们被客户要求撰写关于广义加性模型的研究报告,包括一些图形和统计输出. 视频:R语言广义相加模型(GAM)在电力负荷预测中的应用 拓端tecdat:R语言广义相加模型(GAM)在电力负荷预测中的应 ...
- R语言学习笔记之——数据处理神器data.table
前言 数据处理在数据分析流程中的地位相信大家都有目共睹,也是每一个数据从业者面临的最为繁重的工作任务. 在实际应用场景下,虽然SQL(SQL类专业的etl语言)是数据处理的首选明星语言,性能佳.效率高 ...
- R语言计算医学统计学中RR、OR和HR三个关于比值
计算队列研究的相对风险度(Calculating Relative Risks for a Cohort Study) 生物医学统计数据中一种非常常见的数据集是一项队列研究,在这项研究中,您可以了解接 ...
- r语言plot函数x轴y轴名字_Matplotlib入门-1-plt.plot( )绘制折线图
在Python学堂1-8中,我们通过一个实例(温湿度变化曲线对比分析图)对Matplotlib模块中的常用折线图进行了初步的认知学习.本章及后续几章内容是对前面内容的总结以及细化.系统性的带领大家共同 ...
最新文章
- Intel® oneAPI Toolkit编译VASP
- ABAP类的继承、多态、重载
- 火焰效果材质实现_「游戏开发」使用Unity实现魔法火焰效果
- git删除本地分支、删除远程分支 复制分支
- 光缆弹性模量计算_光缆的制造、种类、施工、选用方法（超全）
- 信息学奥赛C++语言： 数字方阵
- 疯狂java的第四章答案_疯狂java讲义第四章第一题
- 62. 不同路径（JavaScript）
- opencv源码解析之(5)：CommandLineParser类的简单理解
- 科创板第二批受理名单公布，为何AI独角兽 “全军覆没”？
- ios 推送通知服务证书不受信任（Apple Push Service certificate is not trusted）
- H5公众号-canvas海报分享图+生成二维码
- 【@NotNull和@NotBlank的区别】
- TypeError: xxx() missing 1 required positional argument xxx
- c语言中 hr是什么字符,Heat Treatment 的缩写（HR）C
- 人脸识别门禁有哪些功能呢？
- 个人训练赛第二十场：对撞
- 使用ADK（windows评估和部署工具包）制作PE
- Android MQTT协议
- idea下载安装教程(2020 | 2019)