计算两组患者生存率（考虑截尾数据） R语言代码

在R语言中，计算两组患者的生存率并考虑截尾数据通常涉及到生存分析中的 Kaplan-Meier (K-M) 曲线估计和log-rank检验。下面是一个简单的示例，假设你有两个数据框`data_group1`和`data_group2`，每个数据框都有"Survival_Time"列（生存时间）和"group"列（组别）：

# 首先确保已安装survival包（如果未安装）
if (!require("survival")) install.packages("survival")

# 假设生存时间和分组信息存储在data_group1和data_group2中
group1 <- data_group1$group == "Group1"
group2 <- data_group2$group == "Group2"

# 使用 survfit 函数创建生存函数
km_fit_group1 <- survfit(Surv(data_group1$Survival_Time, event = data_group1$Event_column) ~ group1)
km_fit_group2 <- survfit(Surv(data_group2$Survival_Time, event = data_group2$Event_column) ~ group2)

# 绘制生存曲线
plot(km_fit_group1, main = "Survival Curves for Group 1", xlab = "Time", ylab = "Survival Probability")
lines(km_fit_group2, col = "red", lty = 2) # 第二组的红色虚线

# 计算Log-Rank统计量
lr_test <- survdiff(formula = Surv(Time, Status) ~ group, data = rbind(data_group1, data_group2))
summary(lr_test)

# 输出结果
cat("Log-Rank Test Results:\n", lr_test$p.value, "\n