零基础怎么办?没时间学怎么办?
这份R语言集合让你1分钟轻松掌握R语言,再也不用被导师嫌弃了~!
PS:小编按不同需求分类,请各位看官按需自提!
干货来袭
今日(2019-8-29)----生存曲线绘制
项目基础版中阶版高级版懒人版生存曲线样图
获取方式原代码操作文档+原代码视频+操作文档+原代码点击详解
FUNs首日福利
中阶版原代码 (10行代码就可以完成)
#数据读入及检查------------------------------------------------------------
HCC
HCC [1:5,]
S
class(S)
S
#普通K-M生存图------------------------------------------------------------
survfit(S ~ HCC$Group )
survfit(Surv(HCC$OS, HCC$COD) ~ HCC$Group,data = HCC)
sfit
summary (sfit)
plot (sfit)
#彩色K-M生存图------------------------------------------------------------
#单纯彩色K-M生存图1------------------------------------------------------
ggsurvplot(sfit)
中阶版教程详解
第一步:数据读取。对于新生来说,如何读取数据是R语言的第一大障碍!
1. 新建excel表格---按下图所示录入数据(Group为分组[1=A;2=B],COD表示生存情况[1=死亡;0=存活],OS表示生存时间 [以月份为单位]。
2. 文件----另存为---地址选择桌面(也可以存在自己特定的文件中---文件名:OS---保持类型 CSV(逗号分隔)。
3. 打开RStudio----新建R script
#数据读入及检查------------------------------------------------------------
#数据赋值给HCC就是我们的表格数据
#特别提醒:地址栏中是“/”,而不是“".
HCC
#显示前5行数据
HCC [1:5,]
如下表所示
Group COD OS
1 1 1 15
2 1 0 1
3 1 1 11
4 1 0 10
5 1 1 5
第二部分,生存率的计算和曲线绘制
#准备各自R包
library (survival)
library (dplyr)
library (ggplot2)
library (survminer)
#计算生存率------------------------------------------------------------
#HCC$OS和HCC$COD分别是我们表格中OS和COD那一列数据
#特别提醒:R语言中是严格区分大小的
S
class(S)
S
如下方所示
> S
[1] 15 1+ 11 10+ 5 23 14 11 16+ 4 24+ 27 11+ 9 26+ 18+ 3 13 6 25 13+ 4 24+ 16+
[25] 20+ 1 21 1 19 24 10+ 23 18+ 2 17 6 19 27 21+ 26 26+ 1+ 1 26 14+ 5 21+ 6
[49] 8+ 39 9+ 58 59 17+ 8+ 19+ 30+ 56 14+ 48+ 11+ 10 29 5+ 10+ 29+ 10 40+ 58 22 58 13+
[73] 8+ 6 53+ 25+ 21+ 40 58 27 19 32 8+ 22 10+ 12 7+ 29 8+ 43+ 21+ 23 59 54 45 44
[97] 18+ 18 14
#生存曲线绘制------------------------------------------------------------
survfit(S ~ HCC$Group )
survfit(Surv(HCC$OS, HCC$COD) ~ HCC$Group,data = HCC)
sfit
summary (sfit)
plot (sfit)
生存曲线就出来,黑白不好看,想要彩色怎么办?
第三部分,图片美化,1秒变彩色
#立马变彩图
ggsurvplot(sfit)
生存曲线就出来了,是不是还没看够!
目前越来越多的TOP的杂志都要求95%CI, 和各对应点数据?
读者:小编,小编,我想要图片更美,和你一样美,像高阶版的那种
小编:持续关注,定期送福利,送高阶版教程