学习笔记

Posted 2022-12-30 统计学小王子

tags:

篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了学习笔记相关的知识，希望对你有一定的参考价值。

0引言

开始学习生存分析了，选择了Statistical Analysis of Failure Time Data 2nd¹这本书来学习。为保证学习效果。会把一些知识点以及总结发出来共同学习。水平有限，发现错误希望评论区及时指正。

1、失效时间数据（failure time data）

定义：关心的事件发生在已经出现的事件的之前，称关心的事为失效事件数据。
主要的应用领域：生物医药研究（biomedical studies）或者产业寿命（industrial life）
两个主要的问题：
– 解释变量对失效时间的依赖性。(dependence of failure time on explanatory variables)
– 潜在失效时间分布模型的估计和具体形式。(the estimation and specification of models for the underlying failure time distribution)

1.1 癌症案例

附录

右截断数据展示

plot.right_censored <- function(data)
  group = data[[1]]
  group_censored = data[[2]]
  n <- length(group);Max = max(group) + max(group)*0.05
  plot(x = c(1,group[1]),y = c(1,1),xlim = c(0,Max),
  ylim = c(0,n+1), type = "l", xlab = "x", ylab = "")
  for(i in 1:n)
    lines(c(1,group[i]), c(i,i), type = "l")
    if(group_censored[i] == 1) 
      points(group[i], i, pch = 8)
      text(group[i]+group[i]*0.04, (i-0.1),"右截断")
    
  
  abline(v = 0, lty = 2)
  text(20,n/2+0.5,"个体感染时间")

主函数

# 导入自编函数
source("right_censored.r")
# 载入包
library(lattice)
library(ggplot2)
# 右截断数据可视化
group1 <- c(143,146,265,304,216,244)  # 原数据
group1_censored = c(rep(0,4),rep(1,2))  # 是否右截断
group2 <- c(142,156,296,323,204,344)
group2_censored = c(rep(0,4),rep(1,2))
censored_Data1 <- data.frame(
  group1 = group1,
  group1_censored = group1_censored
)
censored_Data2 <- data.frame(
  group2 = group2,
  group2_censored = group2_censored
)
class(censored_Data1) <- "right_censored"
class(censored_Data2) <- "right_censored"
par(mfrow = c(2, 1))
plot(censored_Data1)
plot(censored_Data2)