c#读取excel两列数据并绘制xy曲线_R语言生存曲线绘制

零基础怎么办?没时间学怎么办?

这份R语言集合让你1分钟轻松掌握R语言,再也不用被导师嫌弃了~!

PS:小编按不同需求分类,请各位看官按需自提!

干货来袭

今日(2019-8-29)----生存曲线绘制

项目基础版中阶版高级版懒人版生存曲线样图

d3e2be8ddd924cbe8a66203369f8bd08
0b1a1659947547c6a2ce9c373d467942
90ad6166acc34948a90ab01b5f29766e

获取方式原代码操作文档+原代码视频+操作文档+原代码点击详解

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表示生存时间 [以月份为单位]。

7e20fe30cbca4469a20934135c860fbf

2. 文件----另存为---地址选择桌面(也可以存在自己特定的文件中---文件名:OS---保持类型 CSV(逗号分隔)。

b931307179294e0c8940a0292b5d93ba

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)

d3e2be8ddd924cbe8a66203369f8bd08

生存曲线就出来,黑白不好看,想要彩色怎么办?

第三部分,图片美化,1秒变彩色

#立马变彩图

ggsurvplot(sfit)

0b1a1659947547c6a2ce9c373d467942

生存曲线就出来了,是不是还没看够!

目前越来越多的TOP的杂志都要求95%CI, 和各对应点数据?

读者:小编,小编,我想要图片更美,和你一样美,像高阶版的那种

小编:持续关注,定期送福利,送高阶版教程

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值