脚本之家 服务器常用软件
快捷导航
您的位置：首页 → 软件编程 → R语言 → R语言生存分析

R语言关于生存分析知识点总结

 更新时间：2021年05月17日 08:10:31   作者：w3cschool  
在本篇文章里，小编给大家整理的是一篇关于R语言生存分析的相关知识点及实例内容，有兴趣的朋友们跟着学习下吧。

生存分析处理预测特定事件将要发生的时间。 它也被称为故障时间分析或分析死亡时间。 例如，预测患有癌症的人将存活的天数或预测机械系统将失败的时间。

命名为survival的R语言包用于进行生存分析。 此包包含函数Surv（），它将输入数据作为R语言公式，并在选择的变量中创建一个生存对象用于分析。 然后我们使用函数survfit()创建一个分析图。

安装软件包

install.packages("survival")

语法

在R语言中创建生存分析的基本语法是

Surv(time,event)
survfit(formula)

以下是所使用的参数的描述 

  • time是直到事件发生的跟踪时间。
  • event指示预期事件的发生的状态。
  • formula是预测变量之间的关系。

例

我们将考虑在上面安装的生存包中存在的名为“pbc”的数据集。 它描述了关于受肝原发性胆汁性肝硬化（PBC）影响的人的生存数据点。 在数据集中存在的许多列中，我们主要关注字段“time”和“status”。 时间表示在接受肝移植或患者死亡的患者的登记和事件的较早之间的天数。

# Load the library.
library("survival")

# Print first few rows.
print(head(pbc))

当我们执行上面的代码，它产生以下结果及图表

  id time status trt      age sex ascites hepato spiders edema bili chol
1  1  400      2   1 58.76523   f       1      1       1   1.0 14.5  261
2  2 4500      0   1 56.44627   f       0      1       1   0.0  1.1  302
3  3 1012      2   1 70.07255   m       0      0       0   0.5  1.4  176
4  4 1925      2   1 54.74059   f       0      1       1   0.5  1.8  244
5  5 1504      1   2 38.10541   f       0      1       1   0.0  3.4  279
6  6 2503      2   2 66.25873   f       0      1       0   0.0  0.8  248
  albumin copper alk.phos    ast trig platelet protime stage
1    2.60    156   1718.0 137.95  172      190    12.2     4
2    4.14     54   7394.8 113.52   88      221    10.6     3
3    3.48    210    516.0  96.10   55      151    12.0     4
4    2.54     64   6121.8  60.63   92      183    10.3     4
5    3.53    143    671.0 113.15   72      136    10.9     3
6    3.98     50    944.0  93.00   63       NA    11.0     3

从上述数据，我们正在考虑分析的时间和状态。

应用Surv（）和survfit（）函数

现在我们继续应用Surv()函数到上面的数据集，并创建一个将显示趋势图。

# Load the library.
library("survival")

# Create the survival object. 
survfit(Surv(pbc$time,pbc$status == 2)~1)

# Give the chart file a name.
png(file = "survival.png")

# Plot the graph. 
plot(survfit(Surv(pbc$time,pbc$status == 2)~1))

# Save the file.
dev.off()

当我们执行上面的代码，它产生以下结果及图表

Call: survfit(formula = Surv(pbc$time, pbc$status == 2) ~ 1)

      n  events  median 0.95LCL 0.95UCL 
    418     161    3395    3090    3853

，使用R生存分析

上图中的趋势有助于我们预测在特定天数结束时的生存概率。

到此这篇关于R语言关于生存分析知识点总结的文章就介绍到这了,更多相关R语言生存分析内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！

您可能感兴趣的文章:

相关文章

  • R语言实现二进制文件读写操作

    R语言实现二进制文件读写操作

    这篇文章主要介绍了R语言实现二进制文件读写操作，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧
    2021-03-03
  • R语言 数据表匹配和拼接 merge函数的使用

    R语言 数据表匹配和拼接 merge函数的使用

    这篇文章主要介绍了R语言 数据表匹配和拼接 merge函数的使用说明，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
    2021-03-03
  • R语言对CSV文件操作实例讲解

    R语言对CSV文件操作实例讲解

    在本篇文章里小编给大家整理了一篇关于R语言对CSV文件操作实例讲解内容，有兴趣的朋友们可以学习下。
    2021-04-04
  • 利用R语言绘制时间序列图的操作

    利用R语言绘制时间序列图的操作

    这篇文章主要介绍了利用R语言绘制时间序列图的操作，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
    2021-03-03
  • R语言ggplot2实现将多个照片拼接到一起

    R语言ggplot2实现将多个照片拼接到一起

    本文主要介绍了R语言ggplot2实现将多个照片拼接到一起，文中通过示例代码介绍的非常详细，需要的朋友们下面随着小编来一起学习学习吧
    2021-07-07
  • R语言seq()函数的调用方法

    R语言seq()函数的调用方法

    这篇文章主要介绍了R语言seq()函数的调用方法,本文给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下
    2021-03-03
  • R语言对数据库进行操作的实例详解

    R语言对数据库进行操作的实例详解

    在本篇文章里小编给大家分享的是一篇关于R语言对数据库进行操作的实例详解内容，有兴趣的朋友们跟着学习下吧。
    2021-05-05
  • R语言中corrplot标题居中及eps格式输出

    R语言中corrplot标题居中及eps格式输出

    由于论文的需要，很多R语言绘图需要eps格式矢量图的方式进行输出，同时最近也需要绘制相关性图，因此决定写个博客作为记录，有需要的朋友可以借鉴参考下
    2021-11-11
  • 如何用R语言绘制散点图

    如何用R语言绘制散点图

    这篇文章主要介绍了如何用R语言绘制散点图，帮助大家更好的理解和学习使用R语言，感兴趣的朋友可以了解下
    2021-03-03
  • 关于R语言lubridate包处理时间数据的问题

    关于R语言lubridate包处理时间数据的问题

    这篇文章主要介绍了关于R语言lubridate包处理时间数据的问题,本文给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下
    2021-05-05

最新评论