如何使用R语言做生存分析?

首先,我们这里是介绍使用survival这个包来做。

软件的安装:install.packages(survival)

加载包:library(survival)

我们这里利用R自带的数据进行分析。

attach(melanom)

names(melanom)

Kaplan-Meier分析

survfit(Surv(days,status==1))

surv.all<-survfit(Surv(days,status==1))

summary(surv.all)

画生存曲线

plot(surv.all)

接着根据性别进行分组画生存曲线

surv_sex<-survfit(Surv(days,status==1))~sex)

plot(surv_sex,conf.int=T,col=c("red","green"))

接着进行对数秩检验

survfit(Surv(days,status==1)~sex)

Cox分析

summary(coxph(Surv(days,status==1)~sex))

本次分享就到这里,希望本文对大家有所帮助。

(0)

相关推荐

  • R语言生存分析

    欢迎来到医科研,这里是白介素2的读书笔记,跟我一起聊临床与科研的故事, 生物医学数据挖掘,R语言,TCGA.GEO数据挖掘. R语言生存分析  生存分析是医学数据挖掘中的重要内容 R语言中用于生存分析 ...

  • 生信编程16.对有临床信息的表达矩阵批量进行生存分析

    有一些五六年前的学生们都成长为了各个生物信息学相关公司的小领导,而且他们都有了自己的公众号,知乎号,也算是一番人物.最近他们跟我反馈面试找不到或者说很难直接考核筛选到认真干活的生信工程师,挺有意思的. ...

  • 单因素和多因素cox回归分析

    参考 生信交流平台 https://zhuanlan.zhihu.com/p/339754573欢迎关注微信公众号"生信交流平台"6 人赞同了该文章前面我们讲过一个R函数搞定风险评 ...

  • estimate或者CIBERSORT结果真的是很好的临床预后指标吗

    肿瘤免疫微环境我们讲了很多内容了,目录是: estimate的两个打分值本质上就是两个基因集的ssGSEA分析 针对TCGA数据库全部的癌症的表达量矩阵批量运行estimate 不同癌症内部按照est ...

  • R数据分析:生存分析的做法和结果解释

    今天给大家写写生存分析: Survival analysis corresponds to a set of statistical approaches used to investigate th ...

  • 生存分析就是一个任人打扮的小姑凉

    最近接到粉丝的TCGA分析需求,想看看指定基因在指定癌症是否具有临床意义(也就是生存分析是否有统计学显著效果咯!)其实很早以前我在生信技能树就号召粉丝讨论过这个问题:集思广益-生存分析可以随心所欲根据 ...

  • 用R语言做单因素方差分析及多重比较

    SPSS方差分析的应用已经做得非常好了,绝大多数的方差分析问题均可通过SPSS"点菜单"的方式得以解决,R语言在统计和可视化方面有自己的特色, 我们不妨来对比着学习. 选用R语言自 ...

  • 博士师兄手把手教你用R语言做PCA分析,不存在学不会!

    手把手教你用R语言做PCA主成分分析,不存在学不会 作者简介 本文作者Trigo Hoang,作者目前在香港攻读博士学位,硕士期间发表了多篇生信相关的SCI,累计影响因子35+,公众号简书会记录作者学 ...

  • 精心整理(含图PLUS版)|R语言生信分析,可视化

    为了能更方便的查看,检索,对文章进行了精心的整理(PLUS).建议收藏,各取所需,当前没用也许以后就用到了呢! R资料+计划 R语言精品资料年中无套路赠送 R-plotly|交互式甘特图(Gantt ...

  • R语言配对检验分析案例

    原文链接 :http://tecdat.cn/?p=3424 什么是检验对? 检验对的形式 (x1,x2) 出现在两种情况中: 对同一实体执行两次测量.例如,一项评估新型胰岛素疗效的临床研究将为每位患 ...

  • 技术贴 | R语言——肠型分析:介绍、方法

    导读 2011年,肠型(Enterotypes)的概念首次在<自然>杂志上由Arumugam等[1]提出,该研究发现可以将人类肠道微生物组分成稳定的3种类型,因为这3种类型不受年龄.性别. ...

  • 技术贴 | R语言物种组成分析和绘图

    本文由阿童木根据实践经验而整理,希望对大家有帮助. 原创微文,欢迎转发转载. 导读 宏基因组分析分为物种分析和功能分析两大块.物种组成分析是物种分析中最基本最常见的分析方法.利用R语言堆叠图,我们可以 ...

  • R语言数据可视化分析案例:探索BRFSS数据

    原文链接:http://tecdat.cn/?p=9284 加载包 library(tidyr) library(knitr) opts_chunk$set(echo = TRUE, fig.alig ...

  • BMJ子刊:手把手教你用R语言实现meta分析~

    Meta分析是循证医学中非常重要的工具,如果配上一个好的科学问题,它可以整合当前已有的证据,以量化的方式给出相对公正的结论,临床价值也可以非常巨大! 那如何使用R做Meta分析? 是不是觉得无从下手? ...